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SE  CTION  I 


INTRODUCTION 


This  report  covers  the  primary  technical  work  completed  during 
the  first  six  months  of  the  subject  contract  to  design  and  develop  a 
prototype  moving  base  gravity  gradiometer.  The  scope  of  work  during 
this  time  period,  and  for  the  remaining  time  during  Phase  One,  is 
limited  to  establishing  the  instrument  and  system  design  parameters. 

As  such  the  nature  of  this  work  is  restricted  primarily  to  analytical  and 
theoretical  studies.  The  Phase  One  conclusions  will  be  documented  in 
the  R&D  Design  Evaluation  Report  and  presented  at  a  Design  Review. 

This  first  six  months  has  seen  many  positive  accomplishments. 
The  most  significant  findings  are  documented  on  the  following  pages 
of  this  report.  Because  of  the  complex  and  interactive  nature  of  the 
many  tasks  required  under  this  contract,  and  because  this  report  covers 
a  specific  time  span,  it  is  necessary  to  regard  some  of  the  findings  as 
being  less  than  absolute,  thus  open  to  revision  as  the  Phase  One  study 
proceeds.  The  progress  has  been  in  general  accordance  with  the 
original  management  and  technical  plan.  The  only  significant  deviation 
from  the  September  1971  Proposal  has  been  in  the  mix  between  analy¬ 
tical  and  laboratory  work;  there  have  been  almost  no  laboratory  exper¬ 
iments  conducted  with  the  original  hard- bearing  sensor  during  this 
reporting  period.  This  decision  results  from  desires  expressed  by  the 
AFCRL  Contract  Monitor  and  from  a  reestablishment  of  task  priorities 
due  to  funding,  personnel  and  time  limitations. 

Significant  technical  progress  and  accomplishments  toward 
finalization  of  the  design  of  a  prototype  moving  base  rotating  gravity 
gradiometer  (RGG)  have  been  achieved,  most  of  which  are  documented 
on  the  following  pages. 
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SECTION  II 


SUMMARY 


During  the  first  six  months  of  work  under  the  contract,  HRL 
has  conducted  analyses,  studies,  and  originated  several  preliminary 
designs.  There  has  been  an  evolution  of  baseline  sensor  configura¬ 
tions  beginning  with  the  spherical  floated  sensor  concept  presented  in 
our  September  1971  proposal.  We  have  conducted  individual  error 
analyses  and  have  started  on  a  comprehensive  error  budget.  We  have 
selected  many  of  the  basic  design  parameters  for  the  subsystems  of 
the  gradient  tensor  sensing  system.  The  principal  technical  findings 
and  accomplishments  are  described  briefly  in  this  summary.  Detailed 
discussions  of  each  topic  follow  in  the  body  of  this  report  (Section  III), 
entitled  "Technical  Findings  and  Accomplishments.  " 

Section  III-A  presents  a  complete,  consistent  and  organized 
system  error  budget,  including  a  listing  of  all  known  error  excitation 
sources  and  error  mechanisms.  This  detail  is  necessary  for  a  system 
as  complex  and  difficult  as  the  moving  base  gravity  gradiometer. 

This  section  provides  a  tabular  format  which  indicates  errors  reflected 
at  the  individual  sensor  output  as  well  as  at  the  individual  gradient 
tensor  element.  Due  to  the  comprehensive  nature  of  this  task,  the 
table  is  not  yet  completed,  however  many  errors  of  individual  sub¬ 
systems  have  been  estimated  and  appear  in  other  sections  of  this 
report.  In  addition.  Section  III-A2  presents  RGG  coordinate  defini¬ 
tion  ri  which  will  be  adopted  for  all  system  error  analyses,  budgeting 
and  system  design. 

Derivations  of  particular  sensor  error  mechanisms  are  presen¬ 
ted  in  Section  III-B.  The  thermad  noise  de’-ivation  (III-Bl)  corrects  an 
error  in  our  previously  reported  work.  The  new  thermal  noise  analy¬ 
sis  has  a  direct  impact  on  the  basic  sensor  arm  size  requirement. 

Our  previous  analysis  indicated  for  a  sufficiently  small  thermal  noise 
contribution  (1/3  E.  U.  ),  that  a  sensor  arm  inertia  of  18,  600  gm-cm^ 
W’as  required.  The  corrected  result  indicates  a  requirement  of 
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35,  610  gm-cm  .  Assuming  similar  arm  geometry  and  material 
densities,  this  reflects  in  an  arm  size  increase  of  13.  8  percent. 

A  derivation  of  the  errors  doe  to  rotational  field  effects  (111-52) 
Illustrates  the  importance  of  considering  correlations  between  angular 
rate  about  orthogonal  axes.  This  results  in  adding  considerable  com¬ 
plexity  in  computation  of  estimates  of  errors  caused  by  sensor  case 
angular  disturbances. 

Errors  due  to  arm  mass -unbalance  are  derived  (III-B3)  as  well 
as  the  sensor  error  sensitivity  caused  by  interaction  of  the  Earth's 
magnetic  field  with  paramagnetic  materials  in  the  sensor  arms  (III-B4). 

Design  concepts  and  parameters  for  the  various  subsystems  are' 

presented  in  Section  TTT_r  •  j 

'  ^  description  of  the  rationale  used 

o  .select  the  nominal  design  values  for  the  sensor  resonant  frequency 
y  and  mechanical  amplification  at  resonance  (O)  (111-Cl).  Selection 
o  these  basic  sensor  parameters  along  with  a  thermal  noise  design 

of  1/3  E.  U.  also  defines  the  minimum  sensor  arm  inertia,  hence 
us  minimum  size.  With  the  above  selection  of  basic  sensor  design 
parameters  completed,  many  parameters  of  the  other  subsystems  are 
then  bounded.  For  example,  the  design  requirements  of  the  RGG  elec¬ 
tronic  circuits  presented  in  the  next  section  (111-C2)  were  heavily 
influenced  by  the  basic  sensor  parameter  selection. 

a  d  •  ion  of  the  sensor  arm  mass 

nd  inertia  balance  probler.  .  Values  of  arm  mass  and  moment  of 

inertia  balance  which  are  felt  to  be  practically  achievable  in  the  labora- 
ory  and  maintainable  in  operational  usages  are  shown.-  The  change  of 
balance  resulting  from  the  sensor  being  spun-up  are  estimated  and  it  is 
shown  that  these  changes  are  within  the  balance  adjustment  range  of 
our  proposed  mercury  balance  tubes.  The  resultant  sensor  mass 
unbalance  error  sensitivity  to  vibrations  at  one  and  three-times  spin 
requency  is  shown  to  produce  errors  sufficiently  below  the  required 
error  budget  limit  for  the  anticipated  sensor  vibration  environment. 

Estimates  for  the  sensor  translational  and  angular  vibration 
environment  are  discussed  in  the  next  section  (I11-C4I.  These  esti¬ 
mates  assume  the  carrying  vehicle  to  be  a  KC-135  aircraft  operated 
during  times  when  the  flight  condition  is  no  more  severe  than  a  "mild 


turbulence"  condition.  Vibration  power  spectral  aensities  are  estimated 
for  the  aircraft,  the  base  of  the  inertial  platform,  and  the  stable  element 
of  the  inertial  platform.  To  generate  these  estimates,  preliminary 
performance  parameters  of  both  the  vibration  isolation  mount  and  the 
inertial  platform  were  needed.  To  determine  the  performance  param¬ 
eters  of  each  of  these  subsystems,  some  preliminary  design  concepts 
for  each  subsystem  were  generated.  Results  of  the  vibration  estimates 
reveal  that  the  translational  vibration  environment  of  the  stable  ele¬ 
ment  of  the  inertial  platform  is  sufficiently  small  such  that,  with  the  arm 
mass  balance  capability  we  expect  to  achieve,  the  system  performance 
goal  can  easily  be  met.  The  results  also  indicate  that  if  the  sensor 
were  hard-mounted  to  the  platform  stable  element,  the  rotational  field 
error  would  be  approximately  2-1/2  EU,  one  sigma.  Since  2-1/2 
EU  exceeds  our  budgeted  allocation  for  this  error,  two  alternatives 
can  be  considered.  One  alternative  is  to  add  an  additional  stage  of 
angular  isolation  between  the  platform's  stable  element  and  the  sensor. 
Several  sensor  angular  isolation  support  systems  have  been  considered 
and  some  preliminary  designs  are  presented  in  Section  III-C5.  The 
second  alternative  is  to  reduce  the  angular  vibration  levels  of  the 
platform  stable  element.  The  primary  cause  of  the  2-1/2  EU 
error  is  due  to  the  assumed  platform  gimbal  bearing  friction.  The 
friction  level  was  based  on  the  use  of  conventional  ball-type  platform 
gimbal  bearings.  One  method  of  reducing  this  friction  level  would  be  to 
use  hydrostatic  air  bearings  instead  of  ball  bearings. 

In  the  next  section  (III-C5)  preliminary  designs  are  shown  for 
each  of  the  following  sensor  angular  vibration  isolation  concepts: 

1.  Rotating  floated  spherical  sensor. 

2.  Two  axis  mechanical  gimbal  ring  sensor 
mount. 

3.  Two  axis  air  pad  gimbal  mount. 

The  basic  sensor  size  shown  for  each  of  these  preliminary 
design  sketches  was  approximately  the  same.  The  comparison  table, 
contained  in  this  section,  summarizes  the  advantages  and  disadvantages 
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of  each  concept.  System  No.  3  above,  appears  the  most  favorable 
of  the  three  at  this  time.  Section  III-C6  presents  a  proof  that  only 
four  air  gimbal  pads,  located  at  the  intersection  of  the  corners  of  a 
circumscribed  equilateral  tetrahedron,  provide  an  isoelastic 
suspension  system. 

Section  III-C7  deals  with  the  characteristics,  requirements, 
and  error  mechanisms  which  affect  the  sensor  spin-bearing  perfor¬ 
mance,  Included  is  the  current  spin-bearing  procurement  specifica¬ 
tion,  Most  of  these  analyses  are  equally  applicable  in  bounding  the 
rctor  spin  motor  and  its  servo  loop  performance. 

The  last  topic  of  Section  C  {III-C8)  is  a  summary  of  the  pres¬ 
ently  conceived  RGG  baseline  sensor  parameters.  The  values  shown 
are  not  necessarily  final  firm  values  but  large  deviations  are  not 
expected  during  the  remaining  Phase  I  design  effort. 

Section  D  summarizes  our  efforts  in  finding  and  selecting 
suitable  materials  for  constructing  the  gradiometer.  Because  of  the 
dual  requirements  for  use  of  materials  which  are  both  highly  mechani¬ 
cally  stable  and  highly  non -magnetic,  considerable  effort  was  required 
to  discover  and  locate  sources  of  supply  for  some  of  these  materials. 

It  is  important  to  note  that  many  materials  usually  considered  to  be 
non-magnetic  have  completely  intolerable  magnetic  permeabilities 
for  this  application. 

These  technical  findings  and  accomplishments  cover  the  wcrk 
that  has  been  completed  on  the  moving  base  RGG  to  this  date.  A 
number  of  similar  studies  and  designs  are  in  progress  and  their 
results  will  be  included  in  the  final  report  on  Phase  I. 
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SECTION  III 


TECHNICAL  FINDINGS  AND  ACCOMPLISHMENTS 


The  following  sections  treat  various  aspects  of  the  moving  base 
rotating  gravity  gradiometer  design  problem.  Appropriate  sections 
are  generally  grouped  together,  but  due  to  the  complexity  and  inter¬ 
relations  of  the  problems,  s'ome  sections  discuss  several  subjects. 

A.  ERROR  BUDGET 

1.  Outline  for  RGG  System  Error  Budget 

2.  RGG  Coordinate  Definitions 

These  sections  provide  a  coordinated  and  systematic  approach 
to  the  RGG  error  analysis  and  error  budget. 

1.  Outline  for  RGG  System  Error  Budget 

This  section  prefients  an  outline  and  format  for  the  RGG  System 
Error  Budget.  It  was  prepared  both  to  aid  in  organizing  the  error 
analysis  and  to  illustrate  an  orderly  method  of  presenting  the  results. 

The  error  budget  is  organized  by  hardware  element;  i.  e.  ,  it 
groups  error  contributions  by  hardware  element.  This  grouping  was 
chosen  —  as  opposed  to  gr  ouping  by  error  type,  error  excitation  source, 
or  error  mechanism  —  so  that  performance  requirements  for  any  given 
hardware  element  can  easily  be  gleaned  from  the  summary  table. 

Figure  1  illustrates  the  system  breakdown  structure  for  a 
typical  operational  gradient  tensor  measuring  system.  It  consists  of 
four  major  systems: 

•  Vibrat^oi?  Isolation,  Alignment  and  Leveling  System 

(VIALS).  This  system,  consisting  of  two  major 
subsystems,  provides  a  certain  amount  of  transla¬ 
tional  vibration  isolation,  serves  as  an  attitude 
Reference  and  attenuates  angular  rate  inputs  to  the 
gradiometers. 
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System  Breakdown  Structure. 
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•  Gradiometer(s).  This  system  is  comprised  of  ttiree 
identical  rotating  gravity  gradiometers.  The  hardware 
elements  comprisuig  each  gradiometer  are  detailed  in 
the  figure  and  the  term  "gradiometer”  atr  used  here  is 
intended  to  include  all  of  the  hardware  elements  to  be 
developed  under  the  AFCRL  RGG  contr."  jc.  The  output 
of  the  gradiometer  will  consist  of  the  demodulated  and 
filtered,  in-phase  (I),  and  quadrature-phase  (Q)  sensor 
signals  in  digital  form. 

•  Data  Processing  and  Recording  System.  This  system 
can  only  be  defined  after  a  specific  application  is 
selected.  Typically,  it  would  consist  of  a  digital  cen¬ 
tral  processor,  I/O,  and  recording  system.  Its  func¬ 
tions  would  include  combining  the  outputs  of  the  three 
gradi;)meters  to  form  the  gradient  tensor  elements, 
compensating  for  mass  proximity  effects,  performing 
coordinate  rotations  of  the  gradient  tensor,  computing 
deflection  of  the  vertical,  recording  gradient  data,  etc. 
We  are  planning  to  study  the  data  smoothing  and  com¬ 
bining  of  the  three  gradiometer  outputs  during  Phase  II. 

•  Ancillary  Support  System.  This  system  will  consist 
of  electrical  and  pneumatic  power  supplies,  special 
air  conditioning  equipment  and  the  like  necessary  to 
support  the  rest  of  the  systems. 

a.  Summary  Table  Explanation 

The  proposed  method  of  presenting  the  error  analysis 
is  illustrated  in  T«\ble  I,  "RGG  System  Error  Summary.  "  Explanation 
of  the  table  is  as  follows.  A  single  error  process  is  represented  by  one 
horizontal  line  entry  in  the  table.  As  an  example  of  one  error  process, 
for  the  specified  contributing  subsystem  (e.  g.  ,  spin  bearing)  the 
error  excitation  source  (e.  g. ,  translational  vibration  at  spin  frequency) 
acts  thru  the  error  mechanism  (e.  g.  ,  spin  bearing  variation  in  torque 
with  load  acts  thru  sum-mode  mismatch)  to  produce  an  error.  This 
error  is  then  classified  as  to  error  type  (e.  g. ,  a  bias  shift).  Each 
error  process  is  evaluated  for  either  a  vertical  or  horizontal  spin-axis 
sensor  orientation  and  the  error  in  amplitude  and/or  phase  at  the 
gradiometer  output  is  evaluc  ted.  Further,  the  reflection  of  this  error 
process  at  the  appropriate  ten  lor  element(s)  is  also  evaluated. 

The  following  sections  list  in  detail  the  specific  error  types, 
error  excitation  sources,  and  error  mechanisms. 


9 


10 


T723 


b.  Error  Types 

We  are  concerned  with  the  following  types  of  error  as 
they  affect  the  accuracy  of  measurement  of  each  of  the  elements  of 
the  mass  attraction  gravity  gradient  tensor: 


Initial  Bias  (compensatable) 

Shift  of  Bias  from  Initial  Value  (uncompensatable) 


Noise  (i.  e.  ,  random  variation  about  bias) 
Scale  Factor  Stability 
Scale  Factor  Linearity 


(  Both  of  these  to  include  result¬ 
ant  effect  in  output  in  the 
presence  of  the  Earth's  nominal 
eravitv  gradients. 


c.  Error  Excitation  Sources 

The  following  fifteen  excitation  sources  hopefully 
represent  all  of  those  sources  which  can  possibly  cause  errors  in 
gradient  tensor  output; 


(1)  Translational  Acceleration 
•  Steady  DC 


(2) 


(3) 


•  Transient  changes  in  DC 

•  Vibration  (initial  or  changes  from  iniiial 

levels) 

•  Self-induced  vibration 
Angular  Rates  and  Accelerations 

•  Steady  DC 

•  Transient  changes  in  DC 

•  Vibration  and  changes  from  those  existing 

at  time  of  calibration 

•  Self-induced  vibration 

Temperature  -  (Nominal  operating  temperature  results 
in  thermal  noise  effects) 
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(4) 


Temperature  Variation 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 


d. 

(1) 


Ambient  Pressure 
Ambient  Humidity 
Magnetic  Fields 
Electric  Fields 


Acoustic  Di8*:urbance 
Angular  Orientation 
Mags  Proximity  Noise 
Prime  Power  Variations 
Time  Standard  Variations 


Component  Inherent  Characteristics 

Material  Stability  —  This  will  include  stability  of  dimen¬ 
sional  properties  as  well  as  other  parameter  changes 
(e.g.,  transistor  P's,  Youngs  modulus,  damping  coeffi¬ 
cient,  etc.  )  resulting  from  creep,  aging,  crystal  growth, 
temperature  cycling,  etc. 

Error  Mechanisms 


Translational  Acceleration  Sf.nsitivity 
•  Differential  arm  mass  unbalance 


•  Arm  CG  offset  along  spin  axis 


Transducer  sensitivity  to  a  ,  a  ,  a 

X  y  z 

•  Gage  compression-tension  sensitivity 

•  Compression  induced  torsion 

•  etc. 

•  etc. 

•  Pivot  anisoelasticity 


Arm  anisoelasticity 
Rotor  Mass  Unbalance 


(2) 


(3) 


(4) 


•  Spin  bearing  torque  variations 


• 

Spin  motor  torque  variations 

Speed  control  servo  pickoff 

Sum-mcde  excited  phase  or  bias  shift  error 

Power  input/output  capacitor(s)  torque  variation 

Signal  processing  and  data  recording 

Speed  control  servo 

Angula 

r  Acceleration  Sensitivity 

Rotational  field  error 

Arm  center  of  mass  offset 

Spin  axis  torque  variation 

• 

Signal  processing  and  data  recording 

Thermal  Noise  -  (Sensitivity  to  nominsv!  operating 

temperature) 

• 

Sensor  arm/damping 

• 

Transducer 

• 

Transducer  amplifier 

Temperature  Variation  —  (Sensitivity  to  variations  in 

operating  temperature) 

• 

Arm  differential  mass  unbalance  shift 

• 

Arm  inertia  balance  shift  ' 

, 

sum-mode 
'  mismatch 

• 

Spring  rate  balance  shift  , 

factor 

• 

Arm  inertia  shift 

1  resonant 
►  frequency 

• 

Spring  rate  shift 

^  shift 

• 

Arm  anisoelasticity  shift 

• 

Transducer  gain  or  phase  shift 

•  Difference  mode  damping  shift 
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Electronic  subsystemfl  bias,  gain,  and 
phase  shifts 


(5) 


(6) 

(7) 


(8) 


Frequency  reference  oscillator 
Sensor  amplifier 
FM  modulator  and  transmitter 
FM  receiver  and  demodulator 
Digital  data  processors 
DC  data  processors 
Speed  control  servo 

•  servo  electronics 

•  phase/speed  reference  subsystem 
Ambient  Pressure  Variations 

•  Effect  on  air  bearing  run^iing  torque 

•  Effect  on  sensor  operating  temperature 
Humidity  Variations 

Electromagnetic  Effects 

•  Platform  torquer  interactions 

•  Earth's  fieild  induced  torque  variation  due  to: 

changing  orientation  in  field 

field  magnitude  variation  over  earth 

•  Interactions  from  other  vehicle  sources 
Electrostatic  Effects 

•  Torque  variations  in  power  input/output  capacitors 

•  Charge  build  up  on  rotor 
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(9)  Acoustic  Effects 


•  Acoustically  transmitted  vibration  input  to 
sensor  case 

•  Acoustic  excitation  of  differential  arm  torques 

(10)  Angular  Orientation  Errors 


Sensor 

• 

Phase  changes  after  initial  phase 
calibration 

• 

Spin  axis  alignment 

VIALS 

• 

Platform  initial  alignment  (vertical 
and  azimuth) 

• 

Platform  drift  uncertainty  (vertical 
and  azimuth) 

(11)  Mass  Proximity  Noise 

•  Gimbal  motion 

•  Vehicle  consumables 

•  Vehicle  orientation 

•  Vehicle  borne  mass  motion 

•  Crew 

•  Fuel  slosh 

•  External  masses 

•  Displaced  air/water  density  variation 

•  Sea  wave  motion 

•  Other  vehicle  proximity  (in  van  testing) 
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(12)  Prime  Power  Variatinna 


(13) 


(14) 


Prime  Electric  Power 

•  Speed  control  servo 

•  Digital  data  processors 

•  DC  data  processors 
Air  Bearing  Supply  (if  used) 

•  Sensor  temperatuie 

•  Spin  bearing  torque 

•  Speed  control  servo 
Time  Standard  Variations 

•  Over-all  sensor  phase  and  amplitude 
Component  Inherent  Characteristics 

•  Sensor  Q  variation  with  amplitude 

•  Sensor  resonant  frequency  variation  with 
amplitude 

•  Transducer  nonlinearities 

•  Transducer  amplifier  nonlinearities, 
^^plitudc  and  phase  distortion 

•  Transducer  amplifier  frequency  response 

•  FM  converter  nonlinearity,  distortion, 
and  stability 

•  FM  converter  resolution  limit 

FM  transmitter  nonlinearity,  distortion 
and  stability 

•  FM  receiver  nonlinearity,  distortion 
and  stability 

•  FM  to  digital  conversion 
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•  FM  to  analog  conversion 

•  Digital  data  pi  ^cessing  accuracy 

•  Analog  data  processing  accuracy 

•  Digital  data  recording  accuracy 

•  Analog  data  recording  accuracy 
Materials  Stability 

Same  considerations  as  (4). 

RGG  Coordinate  DefimHo»c 

IZTXTLT  - 

phases  of  the  RGG  d  testing,  evaluation,  and  operational 

t'  es  oi  me  RGG  development  program  it  is 

standard  definitions  for  many  of  the  RG^l  ..r.rneuTVi 
establishes  such  standards  for  the  d.f  • 

for  the  RGG  rotor  and  stator  assemblilT^Til'  coordinates 

general  subject  of  standardization  be  considered^'T""”'''" 

Assembly  Referenr» 

wieK  ehe  notin';  pHrjp':::!'!";'''  ■''"”'''  '■> 

axis  ®  eaken  poaiHve  in  ihe  nominal  dineclion  o,  IhTa" 

Further,  when  viewing  ti,c  RCf  rn  .n 
designate  the  "upper'  o  outsH 

inside  arm  as  No  2  n  ‘  ” 

Of  arm  No.  1  as  T  and  tbar^f  ^rm  NoTar'"”'-'"''”^^^  Principal  inertia 

on  k.  Define  the  principal  inertias  of  each  arm 

as  follows:  C, 

A  =  Minimum  Transverse  arm  inertia 
B  =  Maximum  Transverse  arm  inertia 

G  =  Polar  Inertia 
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Based  on  these  definitions,  the  inertial  tensors  of  each  arm 
are  stated  as  eqs.  (1)  and  (2). 

A  liAj  +  jjBj  +l^Cj  (1 

$2  ^  ii  +  jj  (2 

Because  of  the  symmetry  of  the  arms,  these  definitions  do  not  specify 
a  unique  direction  for  the  transverse  principal  axes  of  each  arm. 

This  choice  must  be  made  at  the  assembly  of  the  rotor;  however,  it  is 
recommended  that  some  scheme  be  devised  to  allow  disassembly  and 
reassembly  of  the  rotor  to  the  identical  reference  in  order  to  provide 
as  much  continuity  as  possible  for  the  instrument  parameters  that 
have  been  determined  on  the  basis  of  the  initially  selected  reference. 

Stator  Assembly  Reference  Frame 

It  is  proposed  that  the  SAR  frame  be  defined  such  that 
its  polar  axis  (z)  coincides  with  the  nominal  spin  axis  of  the  rotor  and 
IS  positive  in  the  direction  the  spin  vector.  Further,  it  proposed 
that  the  transverse  axes  (xy)  of  the  SAR  frame  be  right-handed  on  z  and 
that  coincidence  of  the  RAR  frame  (ijk)  with  the  SAR  frame  (j^)  be 
designated  as  the  mechanical  phase  reference  of  the  rotor  with  respect 
to  the  stator.  This  la^definition  provides  the  basis  for  describing 

the  RAR  unit  vectors  (ijk)  in  terms  of  the  SAR  unit  vectors  (1^)  as 
shown  in  eqs.  (3),  (4)  and  (5). 


X  cos  w  t  +  y  sin  w  t 

S''  s 

(3) 

-X  sin  u>  t  +  y  cos  w  t 
s  s 

(4) 

z 

(5/ 
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C.  Definition  Consequences 

As  a  consequence  of  these  definitions,  the  "body  torques" 
of  the  arms  which  drive  the  RGG  differential  mode  will  have  unique 
time  histories  in  terms  of  stator-referenced  parameters.  For 
example,  let's  consider  the  differential  torque  due  to  gravity  gradient 
and  rotational  field  inputs  to  the  RGG.  T’  ^  gravity  gradient  torque  on 
arm  No.  1  about  its  center  of  mass  is  given  approximately  as  eq.  (6). 


X  r  •  r  dm 


The  integral  is  taken  over  the  mass  of  the  arm,  and  when  it  is  expanded 
in  principal  axes  of  inertia  yields  eq.  (7). 

Lgj  a  r(Cj-Bj)r.^  +  T(A,-c,)r;^  +  k(B,-Aj)r;.  (7) 

The  rotational  field  input  is  just  the  negative  of  the  anisoinertial 
moment  defined  by  eq.  (8). 


L  ^  -0)  X  $,  •  w 

wi  ~  1 


Expansion  of  eq.  (8)  in  principal  axes  of  inertia  yields  eq.  (9). 

=  r  (Bj-Cj)  +  j  (Cj-Aj)  +  k  (Aj-Bj)  w.Wj  (9) 

Finally,  the  gravity  gradient  and  rotational  field  torques  acting  on  each 
arm  about  their  polar  axes  are  expressed  from  eqs.  (7)  and  (9)  as 
eqs.  (10)  and  (11). 

Tj  S  (Bj-Aj)  (r..  -  w.w.)  (10) 

a  (A^-Bj)  (F;.  -  w.»j) 
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The  differential,  torque  may  be  enpreased  as  eq.  (12)  when  the  arm 
trarxsverse  inertias  are  "matched". 


T,  ^  2{B-A) 


r..  -  co.w. 
1  j 


(12) 


stator  eo„7  !  differential  torque  in  terms  of 

s  ator  coordmates  rather  than  rotor  coordinates.  This  may  be  aceom- 

P  IS  e  using  the  coordinate  transformation  defined  by  eqs.  (3)  and  (4) 

he  results  are  stated  as  eqs.  (13),  (14),  and  (15). 


''ij  ■  ^jcy  +  T  (r 


"  2"yy-r^)  Sin2u  t 


w.tii  ~  cj  oj  cos  2a)  t  +  —  /t  1^  •  •> 

-  J  X  y  “  +  2  '  V  " 


(13) 

(14) 


■d  =  (B-A)  -  20)^0)^  j  cos  2a)^t 

+  (B-A)  [(Fyy  -  r^)  +  .  „^)jsin  20)^ 


(15) 


The  differential  mode  torque  of  eq.  (15)  is  acted  on  by  the  RGG  carrier 
xlter  process  to  produce  Ihe  normalized  output  of  eq.  (16). 


2 

p 

0) 

r\ 

2r..  = 

(J 

Q 

Td 

ij 

c4  I  o  c  j.  2 

S  +  -5^  0+0) 

W  O 

B-A 

(16) 

where 


0)  &  2f  - 

O  -  S 
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The  filter  action  of  eq.  (16)  produces  a  90°phase  lag  at  2  w  such  that 
when  the  signals  are  constant,  the  resultant  output  may  be  stated  as 
eq.  (17). 


2a)  (jj 

X  y 


sin  qw  t 
^  s 


( r. 


XX 


-  r  )  +  (w^  -  co^) 
yy  '  y 


cos  2a)  t 
s 


(17) 


Observe  the  sign  reversal  that  has  occurred  between  eqs.  (15)  and  (17) 
in  the  second  part  of  the  equation.  This  result  obviously  provides 
the  cross. gradient  output  signal  in  the  "sine"  demodulation  channel  as 
before;  however,  the  output  of  the  •’cosine"  demodulation  channel 
differs  in  sign  from  previous  analyses. 

Summary  and  Conclusions 

A  set  of  RGG  coordinate  definitions  has  been  proposed 
to  provide  continuity  between  analytical  and  experimental  data.  It  is 
the  intent  that  future  error  studies  adhere  rigorously  to  these  defini¬ 
tions  in  order  that  consistent  statistical  correlation  is  maintained  in 
these  studies  - 


B.  ERROR  MECHANISM  DERIVATIONS 


1.  Thermal  Noise  Fluctuations  in  the  Rotating 
Gravity  Gradiometer 

2.  Covariance  Functions  of  the  Rotational  Field 
Error  of  the  RGG 

3.  RGG  Errors  Due  to  Transverse  Differential 
Arm  Mass  Unbalance 

4.  RGG  Torques  Due  to  Paramagnetic 
Materials 

These  sections  discuss  some  of  the  mechanisms  by  which 
errors  can  be  introduced  into  the  gravity  gradient  tensor  elements. 
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1. 


Thermal  Noise  Fluctuations  in  the  Rotatinfe  Gravity  Gradiometpr 


Intrc'duction 

Physical  devices  whose  function  is  the  measurement 
and/or  processing  of  very  weak  signals  are  fundamentally  performance 
-  limited  by  the  thermal  fluctuations  of  their  dissipative  elements.  The 
specific  force  "signals"  which  the  second  order  gravity  gradiometer 
must  measure  are  extremely  small  (of  the  order  of  lO’^^g  per  cm 
with  an  integration  time  of  the  order  of  ten  seconds).  For  this  reason, 
thermal  noise  must  be  treated  as  a  quantitative  parameter  in  the 
d^jign  of  such  instruments. 

Brief  Conceptual  Background  of  Thermal  Noise 

Over  forty  years  ago,  Nyquist  (Ref.  6)  presented  a 
quantitative  treatment  of  the  thermal  fluctuations  in  linear  electrical 
systems  based  on  the  principles  of  thermodynamics  and  statistical 
mechanics.  More  recently,  a  general  quantitative  treatment  of 
"the  thermal  fluctuations  of  generalized  forces  in  linear  dissipative 
systems,"  embodying  the  Nyquist  result  as  a  special  case,  was  given 
by  Callen  and  Welton  (Ref.  7).  An  important  result  of  this  work  is 
the  "generalized  Nyquist  relation"  which  asserts  that  "a  system  with 
a  generalized  resistance  Rw  exhibits,  in  equilibrium,  a  fluctuaing 
force  V"  whose  mean  square  value  is  given  by  (1). 


<  >  s  2kT  J ^ 


The  right  side  of  eq.  (1)  may  be  viewed  in  terms  of  the  integral  of  a 
model  power  spectrum  of  the  generalized  force  and  re-expressed 

as  a  "two-sided"  integral  in  terms  of  the  circular  rather  than  the 
angular  frequency  as  eq.  (2). 
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(2) 


where 


<  > 


(f)  df 


S^(f)  4  2  kT  R(f) 


(3) 


This  form  of  the  "generalized  Nyquist  relation"  is  based  on  the  equi- 
partition  law,  and  itb  limits  of  validity  are  the  same  as  those  of  the 
equipartition  law,  i.e.  it  is  valid  when  (kT  >>  Kv*). 

Equation  (3)  may  be  employed  to  define  the  model  noise  spectra 
of  the  basic  dissipative  elements  of  electrical  and  mechanical  systems, 
i.e.  ,  electrical  resistance  and  viscous  damping,  as  in  eqs.  (4)  and  (5). 


S  (f)  =  2  kT  R 
e 

(4) 

Sp(f)  =  2  kT  D 

(5) 

Thus,  the  basic  electrical  and  mechanical  dissipative  elements  may  be 
modeled  as  noise-free  components  with  associated  white  noise  genera¬ 
tors  of  voltage  and  force  (or  torque),  respectively,  having  the  power 
spectra  of  eqs.  (4)  and  (5).  These  concepts  form  the  basis  for  the 
thermal  noise  analysis  of  the  rotating  gravity  gradiometer. 

c.  Analytical  Approach 

The  objective  of  this  analysis  is  to  develop  analytical 
expressions  for  the  thermal  noise  induced  errors  of  the  measured 


*k  1.  38  X  10'^^  erg/°K 
H  6.  6  X  10  erg-see 


T  =  Absolute  Temperature,  °K 
V  =  Frequency,  Hz. 
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gravity  gradient  tensor  elements  in  terms  of  the  design  parameters 
of  an  individual  rotating  gravity  gradiometer. 

Measurement  of  all  of  the  gradient  tensor  elements  requires  a 
set  of  three  rotating  gravity  gradiometers  whose  individual  output 
measurements  are  uniquely  combined  to  yield  the  prescribed  tensor 
elements.  The  fact  that  six  measurements  are  available  for  the 
estimation  of  five  independent  variables  suggests  the  possibility  that 
this  apparent  redundancy  might  be  exploited  to  minimize  the  total 
measurement  errors  of  the  gradient  tensor  elements.  Although  worthy 
of  investigation,  this  subject  is  considered  beyond  the  scope  of  the 
thermal  noise  analysis.  For  the  purposes  of  this  analysis,  it  will  be 
assumed  that  the  three  gradiometers  are  mutually  orthogonal  and  that 
their  output  measurements  are  conventionally  combined  to  obtain  the 
tensor  elements.  This  process  provides  the  necessary  "link"  in  this 
analysis  between  the  thermal  noise  errors  of  the  individual  instruments 
and  those  of  the  gradient  tensor  elements. 

The  analytical  approach  consists  of  developing  a  signal  process¬ 
ing  model  for  a  single  instrument,  deriving  the  thermal  noise  errors 
of  this  model,  and  relating  these  errors  to  those  of  the  gravity  gradient 
tensor  elements. 

d.  Signal  Processing  Model 

In  order  to  consider  two  independent  noise  sources  with 
different  points  of  entry  into  the  signal  process,  it  is  convenient  to 
investigate  this  process  in  two  sections:  (1)  the  electromechanical 
portion  associated  with  generating  and  transducing  the  basic  signals, 
and  (2)  the  electronic  portion  associated  with  amplification,  prefilter- 
ingi  phase  sensitive  demodulation,  and  postfiltering  of  the  basic  signals. 
The  two  independent  noise  sources  of  the  model  are  equivalent  torque 
noise  with  a  white  power  spectrum  of  the  form  of  eq.  (5)  and  equivalent 
voltage  noise  with  a  white  power  spectrum  of  the  form  of  eq.  (4).  The 
reason  that  these  sources  are  treated  separately  is  that  they  are 
filtered  by  different  processes.  The  equivalent  torque  noise  is  filtered 
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by  the  electromechanical  filter  of  the  sensor  while  the  equivalent 
voltage  noise  enters  the  signal  process  beyond  the  electromechanical 
filter  and  does  not  receive  the  benefit  of  this  filter  process.  A  simpli- 
fied  functional  block  diagram  of  the  signal  processing  model  is  shown 

as  Fig.  2,  wherein  the  points  of  entry  of  these  independent  noise  sources 

are  indicated  as  n_,  and  n  . 

T  e 


lt4S-l 


Figure  2.  Signal  Processing  Model  Block  Diagram 

The  signal  input  to  the  process  is  the  spatially  modulated  gravity 
gradient  signal  (2r..).  The  filters  H^(s)  and  H^Cs)  are  narrow~band 
filters  centered  at  the  carrier  frequency  of  signal  modulation.  The 
first,  Hj(s),  represents  the  equivalent  electromechanical  filter  of  the 
sensor,  and  the  second,  H^is),  represents  an  electronic  pre-filter. 

The  filtereo  carrier  signal,  2r.j,  represents  the  resultant  signal  plus 
noise  entering  the  linear  demodulators;  these  are  defined  by  the  multi¬ 
pliers,  (2  cos  w^t)  and  (2  sin  u^t).  In  this  analysis  the  demodulation 
frequency,  co^,  is  assumed  to  coincide  with  the  center  frequency  of  the 
carrier  filters.  The  post  filters,  H^(s),  in  each  channel  are  low-pass, 

and  the  fiMl  outputs  of  signal  plus  noise  are  identified  symbolically  as 
r  (t)  and  r  (t). 

C  S 
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In  order  to  make  use  of  the  functional  signal  processing  model 
of  Fig.  2,  it  is  necessary  to  relate  its  parameters  to  those  of  the  physi¬ 
cal  devices.  For  this  purpose,  the  simplified  model  of  the  rotating 
gravity  gradiometer  shown  in  Fig.  3  will  be  employed. 


Figure  3.  Rotating  Gravity  Gradiometer  Simplified  Model. 

The  system  consists  of  the  inertias  of  the  two  "arms"  of  the 
sensor  and  the  three  supporting  springs  and  dampers.  In  addition, 
three  white  noise  torque  generators  (T^,  Tj,  T^)  are  shown  to  be 
associated  with  the  damping  elements  according  to  eq.  (5).  The 
equation  of  motion  of  each  arm  of  the  system  may  be  expressed  as 
eqs.  (6)  and  (7)  in  terms  of  the  principal  arm  inertias  (A,B,  C)  and  the 
spatially  modulated  gravity  gradient  signal 

ce,  +  D^(e,-§2'  +  Ko(6i-V  +  0,6,  +  KjOj 


=  (B-A)  r..  +  T  +  T, 
ij  o  1 

ce^  +  +  K^(92-9,)  +  +  K,92 

=  (A-B)  r..  -  T  +  T, 
ij  o  2 


(6) 


(7) 
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The  "differential  mode"  of  this  system  contains  the  desired  signal 
information,  and  its  differential  equation  is  obtained  by  subtracting 
eq.  (7)  from  (6).  The  result  is  shown  as  eq.  (8),  where  -G2. 

C6  +  (2D^+Dj)  d  +  (2K^+Kj)G 

=  r..  +  2T^  .  T,  .  T2  (8) 

Equation  (8)  may  be  solved  for  the  differential  angle  8  and  then  scaled 

by  the  factor  |c/B.aJ  |e>o/Qil  to  produce  the  normalized  form  of 
eq.  (9). 


"o^°l 

2T  +  T,  -  T  ' 

?  p  +  ®  ^  2 

S  "t  (w  /Q.  )S  +  oj  ^ 

•  O  1  O 

ij  qC 

(9) 


where 


n  ^  (B-A)/C 
Wq/Qj  #  (2D^  +Dj)/C 

"o  =  +K,)/C 

(2fijl*'  4  e(<s^/qQj 

Equation  (9)  is  the  functional  form  required  for  use  in  the  signal 

processing  model  of  Fig.  2.  The  appropriate  identities  are  given  by 
eqs.  (10)  and  (11). 
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(10) 


Hj(S)  4 


£,2  ,  Oca.  2 

S  +  jg-*  S  CO 


+  "^i  -  T,] 


(11) 


The  power  spectral  density  of  ^rp  may  be  found  from  eq.  (5)  using  the 
assumption  that  the  torque  noise  sources  are  statistically  independent. 
Each  generator  has  an  associated  spectrum  StAO  as  follows; 

(f)  =  2kT  (12) 


S._  (f)  =  2kT  D  (13) 

(f)  =  2kT  Dj  (14) 


Because  of  independence,  the  white  power  spectrum  of  may  be 
expressed  as  eq.  (15)  using  eq.  (11). 


N^(o)  = 


■M‘[ 


43„  (f)  +  S„  (f)  +  S„  (f) 


(15) 


Substitution  of  eqs.  (12),  (13),  (14)  into  eq.  (15)  yields  eq.  (16). 


N^(o) 


(16) 
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Using  the  definitions  following  eq.  (9),  a  normalized  form  of  eq.  (16)  is 
shown  as  eq.  (17). 


(17) 


N^(o)  = 


2kT 


2(0 


This  completes  the  modelling  of  the  electromechanical  part  of  the  signal 
process. 

Modelling  of  the  electrical  noise  which  enters  past  the  electro¬ 
mechanical  filter  process  requires  definition  of  an  electrical  signal  scale 
factor  and  an  equivalent  noise  resistance  for  the  specific  circuitry 
involved.  No  attempt  will  be  made  here  to  investigate  the  noise  charac¬ 
teristics  of  specific  circuits.  Rather,  the  form  of  the  power  spectrum 
of  the  electrical  noise  and  its  associated  carrier  filter  are  shown  as 
eqs.  (18)  and  (19). 


N^(c)  4  (2VT  rJ 


H2(S)  4 


(18) 


(19) 


The  parameter  is  the  gradient  scale  factor  associated  with  the  elec¬ 
trical  node  for  which  is  the  equivalent  noise  resistance. 

The  remaining  part  of  the  signal  processing  model  of  Fig.  2 
requiring  definition  is  the  low-pass  filter,  H^(s).  In  this  analysis  a 
simple  second  order  filter  with  equal  roots  will  be  assumed  as  defined 
by  eq.  (20). 


H  (S)  4 
o 


a 

o 

STir 


•i2 


(20) 
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e.  Instrument  Output  Noise 

The  noise  variances  of  the  signal  output  functions,  ^  ,(t) 
and  g(t)»  may  be  determined  from  the  noise  power  spectra  of  these 
functions.  In  turn  these  spectra  may  be  determined  on  the  basis  of  the 
signal  processing  model  of  Fig.  2.  In  effect,  these  spectra  are  the 
result  of  the  filter  process  operating  on  the  white  noise  input  spectra. 
Simply  stated,  these  variances  are  proportional  to  the  amplitudes  of 
the  white  noise  spectra  and  to  the  equivalent  noise  bandwidths  of  the 
signal  processes  which  act  on  the  spectra  as  shown  by  eq.  (21). 


=  2N_,{o)  Af^,  +  2N  (o)  Af 
1  1  e  c 


The  white  power  spectral  amplitudes  are  defined  by  eqs.  (17)  and  (18), 
and  the  equivalent  noise  bandwidths,  Af.^,  and  Af^,  must  be  evaluated  on 
the  basis  of  the  various  operations  in  the  signal  process  model  of  Fig.  2. 
Ecch  equivalent  noise  bandwidth  may  be  determined  independently  of  the 
ocher  by  similar  methods. 

To  illustrate  one  approach,  consider  one  channel  of  the  signal 
process  model  with  the  normalized  white  noise  torque  function,  n.j,(t), 
as  its  input.  This  is  shown  as  Fig.  4. 


It4>-S 


Figure  4.  Signal  Process  Model  with  n^(t)  as  Input. 

The  spectrum  at  the  output  of  the  carrier  filters  Hj(s)  H2(s),  is  simply 
stated  as  eq.  (22). 
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S^(f)  =  N^(o)  Hj(S)H  (S) 


'S  =  jZrrf 


The  spectrum  of  the  output  of  the  demodulation 
derived  as  follows; 


process,  y  (t),  may  be 


Y  (t)  ^  2x(t)  cos  0}  t 

o 


The  autocorrelation  function  of  the  demodulator  output  for  a  etationary 
process  is  just  the  expectation  of  eq.  (24). 

Ry^(r)  -  Ejy^lt)  Y^it  +  r))  ,24, 

From  eq.  (23)  the  desired  product  may  be  written  eq.  (25). 

^c^*^  ~  x(t  +  T)  pos  to^T  +  cos  w^(2t  +  T)j  (25) 

Therefore,  from  eqs.  (24)  and  (25)  the  autocorrelation  of  y  (t)  is 
eq.  (26).  ’ 


R  (T)  =  2R  (T) 


cos  (a)  T 


The  Fourier  traneform  of  eq.  (26)  yields  the  power  spectrum  of  the 
demodulator  output  as  eq.  (27). 

^  ■^®x('-fo)  (27) 

Note  that  the  operation  of  eq.  (27)  is  merely  a  frequency  shift  of  the 
spectrum  of  the  carrier  filter  output,  x(t).  This  process  is  illustrated 
in  Fig.  5.  Note  that  the  demodulator  output  spectrum  has  amplitude 
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^ig.  5(a).  Torque  Noise  Spectrum. 


_ 1 

"  “  7! 

1 

jr - N,(oi 

-f. 

1 

- 

0 

Fig.  5(b).  Carrier  Filter  Output  Spectrum. 


2N  W  at  .ero  fre,uency.  This  is  the  result  of  contributions  from  both 
x(  o)  zero  frequency  region. 

This  ““‘put  is  fUtered  by  the  low-pass  filter,  H  (s). 

Th.s  process  removes  the  second  harmonics  of  the  carrier  and  provides 

aermg  on  the  low  frequency  spectrum  of  the  demodulator  output.  The 

requency  part  of  the  demodulator  output  spectrum  may  be  expressed 

m  terms  of  an  equivalent  low  pass  filter  parametricaily  related  to  the 
carrier  filters,  Hj(s),  as  in  (28). 
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2 


(28) 


S  (f)  s  2N^(o) 


Heq(j2.f) 


After  low  pass  filtering,  the  power  spectrum  of  the  output  function, 
x^(t),  may  be  expressed  from  eqs.  (28)  as  (29). 


S  (f) 

X 


s  2N^(o) 


H^^(j2TTf)  H^(j2TTf) 


(29) 


The  output  variance  is  obtained  by  integration  of  eq.  (29)  over  all 
frequencies  as  in  eq.  (30). 


+  C0 , 


^  2Nrj,{o)J  H^q(j2TTf)  H^(j2TTf) 


df 


(30) 


The  integral  portion  of  eq.  (30)  is  defined  as  the  equivalent  torque  noise 
bandwidth 


+  00 

s  J 

-00 

Note  that  the  definition  of  eq.  (31)  when  substituted  into  eq.  (30)  yields 
one  part  of  the  assertion  of  eq.  (21).  The  other  channel  of  the  signal 
process  containing  the  sin  w^t  demodulation  function  may  be  analyzed 
in  the  same  way  to  yield  the  same  result  as  stated  by  eq.  (32). 

o-g  -  2N^(o)  Af^  (32) 

The  equivalent  noise  bandwidth  (Af^)  acting  on  the  electrical  noise  may 
be  computed  as  in  (31)  by  using  the  equivalent  of  H2(s)  rather  than  the 
equivalent  of  Hj(s)  H2(s)  in  the  integral. 


H^^(j2rf)  H^(j2,rf) 


df 


(31) 
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The  equivalent  noise  bandvridths,  A£^,  may  be  related  to 

the  various  filter  parameters  in  terms  of  the  "equivalent  integration 
times"  of  the  filter  processes  as  defined  by  eqs.  (33)  and  (34). 


Af^  # 

1/2t^ 

(33) 

Af  ^ 
e 

1/2t 

e 

(34) 

To  a  first  approximation,  these  equivalent  integration  times  may  be 
expressed  as  the  sum  of  the  equivalent  time  constants  of  the  various 
filters  in  the  signal  process  as  shown  in  eqs.  (35)  and  (36). 


Tt  =  Tj  +  +  2t^ 


ST  +  2t 
2  o 


where 


^i  ■  ^^l^“o  '  refer  to  eq.  (10) 
\  ^ Q  >  refer  to  eq.  (19) 

■■■q  =  to  eq.  (20) 


The  instrument  output  variance  defined  by  eq.  (21)  may  be 
expressed  explicitly  as  a  function  of  the  instrument  parameters  by 
using  the  power  spectral  amplitudes  of  eqs.  (17)  and  (18)  and  the 
equivalent  bandwidth  relations  of  eqs.  (33),  (34),  (35)  and  (36).  This 
result  is  stated  as  eq.  (37). 
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(37) 


The  standard  deviation  of  the  instrument  output  functions  due  to  torque 
noise  only  is  stated  as  eq.  (38). 


-  1  /  2kT 

■  h  Y 


(38) 


f .  Gravity  Gradient  Tensor  Noise 

To  relate  the  thermal  noise  errors  of  the  individual 
instruments  of  a  system  to  the  resultant  errors  of  the  individual 
gravity  gradient  tenso:  elements,  it  is  necessary  to  define  the  process 
to  be  employed  for  tensor  element  computation.  For  this  purpose, 
consider  a  set  of  three  mutually  orthogonal  instruments  whose  spin 
axes  coincide  with  a  measurement  reference  frame  (XYZ).  The  outputs 
of  each  instrument  contain  signal  plus  noise  as  follows: 


A 

r  =  r  -  r  +  ti 

cx  zz  yy  cx 


A 

r  =  2  r  +  n 

sx  yz  sx 


A 

r  =  r  .  r  +  ^ 

cy  XX  zz  cy 


A 

r  =  2  r  +  n 

sy  xz  sy 


A 

r  =  r  .  r  +  ti 

cz  yy  XX  cz 


A 

r  =  2  r  +  n 

s  z  xy  s  z 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 
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Equations  (40),  (42),  and  (44)  contain  the  cross-gradient  elements,  and 
these  instrument  outputs  when  divided  by  two  yield  the  measured  cross¬ 
gradient  elements  as  shown  in  eqs.  (45),  (46),  and  (47). 


A  1 

r  =  r  +4-0 

yz  yz  2  sx 


A  1 

r  =  r  +  4-  Ti 

xz  xz  2  sy 


(45) 

(46) 


1 


A 

r  =  r  +  T  0 

xy  xy  2  s  z 


(47) 


To  obtain  the  trace  elements  of  the  tensor,  one  must  employ  eqs.  (39), 
(41),  and  (43)  in  pairs  and  use  the  relation  of  eq.  (48). 


r  +  r  +  r  =  o 

XX  yy  zz 


(48) 


For  example,  substraction  of  eq.  (41)  from  (39)  yields  eq.  (49). 


A  A  ■ 

r  -  r  =  2  r 

cx  cy  zz 


(r  +  r  )  +  0  ~  r\ 

XX  yy  cx  'ey 


(49) 


Substitution  of  eq.  (48)  into  (49)  yields  eq.  (50). 


A  A 

r  -  r  =3r  +  r\ 

cx  cy  zz  cy 


(50) 


Division  of  eq.  (50)  by  three  yields  one  of  the  measured  trace  ele¬ 
ments  as  eq.  (51). 


A 

r  =  r  +  ^ 

zz  zz  3 


J  ["^cx  "  ’^cy] 


(51) 
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The  other  trace  elements  may  be  obtained  by  a  similar  process,  and 
the  results  ai-e  stated  as  eqs.  (52)  and  (53). 


^xx  "  ^xx  ^  T  [’^cy  “  ’^cz] 

~  ^yy  ^  ^  cz  “  ’^cx] 


On  the  assumption  that  the  thermal  noise  errors  of  each  instrument 
have  equal  variances  and  are  statistically  independent,  the  standard 
deviations  of  the  tensor  element  measurements  may  be  written  directly 
from  the  functional  relations  of  eqs.  (45),  (46),  (47)  for  the  cross - 
gradient  elements  and  eqs.  (51),  (52),  (53,  for  the  trace  elements. 

These  results  are  stated  as  (54)  for  the  standard  deviation  of  the  cross¬ 
gradient  elements  and  as  (5  5)  for  the  trace  elements. 


a 


ij 


(54) 


O' . . 
11 


(55) 


The  standard  deviations  of  the  instrument  outputs,  o  and  o  ,  are 

6  9 

both  equal  to  the  square  root  of  eq.  (37). 

g*  Summary  and  Conclusions 

The  thermal  noise  errors  of  the  rotating  gravity  gradi- 
ometer  have  been  analyzed  on  the  basis  of  the  "generalized  Nyquist 
relation"  derived  by  Callen  and  Welton.  A  signal  processing  model  was 
developed  for  an  individual  instrument,  and  its  spectral  response  to 
both  mechanical  and  electrical  thermal  noise  was  analyzed.  Expres¬ 
sions  for  the  instrument  output  variances  were  developed,  and  these 
variances  were  employed  in  a  model  gravity  gradient  system  to  obtain 
the  standard  deviations  of  the  resultant  gravity  gradient  tensor 
elements. 
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A  fundamental  result  of  this  study  is  that  the  instrument  output 
variance  due  to  torque  noise  is  inversely  proportional  to  the  product  of 
the  sr riser's  electro-mechanical  time  constant  (t^)  and  the  total  equiva¬ 
lent  time  constant  of  the  filter  process 

2.  Convariance  Functions  of  the  Rotational  Field  Errors 
of  the  Rotating  Gravity  GradiometeF 

a.  Introduction 

The  second  order  gradient  of  the  specific  force  field  at 
a  point  in  a  gravity  field  (as  viewed  l^y  an  observer  in  a  rotating  frame 
of  reference)  is  a  linear  combination  of  the  second  order  gradients  of 
the  gravity  and  rotational  fields.  All  so-called  gravity  gradiometers 
are  actually  "specific  force"  gradiometers;  and  because  Df  this,  the 
measurements  of  such  instruments  are  "contaminated"  by  the  inertial 
augular  velocity  of  their  measurement  frames  of  referenc:e.  In  order 
to  obtain  the  gravity  gradient  tensor  elements  from  these  measurements, 
it  is  necessary  to  correct  the  measurements  for  the  rotational  field 
effects.  The  differences  between  the  actual  rotational  field  gradients 
and  the  quantities  employed  for  compensation  are  defined  as  the  rota¬ 
tional  field  measurement  errors. 

The  intent  of  this  memorandum  is  to  develop  the  convariance 
iunctions  of  the  rotational  field  errors  at  the  input  to  a  single  rotating 
gravity  gradiometer  under  the  constraint  that  angular  velocities  of  the 
RGG  are  sample  functions  of  wide  sense  stationary  random  processes 
with  jointly  gaussian  probability  density  functions.  In  addition  a  method 
for  determining  the  variances  of  the  RGG  rotational  field  output  errors 
is  presented. 

b.  RGG  Rotational  Field  Input  Errors 

F rom  the  viewpoint  of  spectral  analysis,  the  rotational 
field  errors  of  the  RGG  are  most  conveniently  treated  in  the  non¬ 
spinning,  measurement  reference  frame  of  the  instrument.  This 
approach  allows  formulation  of  the  equivalent  input  errors  to  each  channel 
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of  the  RGG  signal  process  and  determination  of  the  resultant  output 
error  based  on  an  equivalent  RGG  signal  process  filter  operation  on 
each  input  error. 

The  RGG  rotational  field  error!"  are  functions  of  the  components 
of  the  RGG  case  referenced  inertial  angular  velocity  normal  to  the 
RGG  spin  reference  axis.  These  angular  velocity  components  may  be 
expressed  as  the  sums  of  a  zero-mean  random  variable  and  a  mean 
value  (in  the  statistical  sense)  as  eqs.  (1)  and  (2). 

cj  ^  X  +  M  (1) 

CJ  =  Y  +  M  (2) 

y  y 

The  "XY"  designations  refer  to  the  RGG  case-reference  axes  normal 
to  its  spin-reference  axis. 

In  this  analysis,  it  will  be  assumed  that  the  rotational  field 
effects  are  compensated  by  functions  of  the  statistical  means  of  the 
angular  velocity  components  such  that  the  rotational  field  input  errors 
will  be  defined  for  each  channel  as  eqs.  (3)  and  (4). 


E 

(cjJ  ■ 

.  +  (mJ  -  M^) 

(3) 

c 

y 

X  X  y 

E 

A 

2cj  CJ 

-  2M  M 

(4) 

s 

X 

y  X  y 

Expansion  of  eqs.  (3)  and  (4)  using  (1)  and  (2)  yields  (5)  and  (6)., 

E  =  Y^  -  X^  +  2  (M  Y  -  M  X)  (5) 

c  y  X  '  ' 

E  =  2X')  +  2  (M  X  +  M  Y) 

S  '  y  X 


(6) 


These  equations  are  formulated  completely  in  terms  of  zero-mean 
random  variables  with  the  mean  values  of  the  angular  velocities  serving 
only  as  coefficients.  The  mean  values  of  these  input  errors  aro  obtained 
directly  from  eqs.  (5)  and  (6)  as  (7)  and  (8). 

^  [^c]  ®  ®  <’> 

=  E  [eJ  =  2E  [xy]  (8) 

These  means  may  be  expressed  also  in  terms  of  the  covariance  functions 
of  the  random  variables  evaluated  at  zero  time  shift  as  in  eqs.  (9) 
and  (10). 


M  =  C  (o)  -  C  (o)  (9) 

c  y  X  '  ' 

M  =  2  C  (o)  (10) 

s  xy'  '  ' 

c.  Input  Error  Covariance  Functions 

The  input  error  covariance  functions  may  be  expressed 
in  terms  of  their  autocorrelation  functions  and  mean  values  as  in 
Eqs.  (11)  and  (12). 


Cc  <->  =  Rc  •  *^0 
C,  (->  =  R.  (rl  - 


(11) 

(12) 


The  input  error  autocorrelation  functions  may  be  expressed  in  general 
as  the  statistical  expectations  of  the  time-shifted  products  of  each  error 
function  as  in  eqs.  (13)  and  (14). 


Rc(ti,t2)  =E  [Ec(ti)Ec(y]  (13) 

RJtj.y  =.E  [e^(1,)E^(1^)]  (14) 
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An  expansion  of  eqs.  (13)  and  (14)  based  on  (5)  and  (6)  will 
produce  a  sum  of  expectations  of  the  products  of  zero  mean  random 
variables  in  various  combinations.  In  general,  the  joint  probability 
density  function  of  each  of  these  combinations  must  be  known  in  order 
to  evaluate  all  of  these  expectations.  In  this  analysis  all  combinations 
of  the  zero-mean  random  variables  are  assumed  to  have  jointly  gaussian 
probability  densities.  Within  this  constraint,  we  may  employ  a  general 
expansion  formula  (Ref.  8)  for  the  expectation  of  the  product  of  "n" 
z^ro-mean,  jointly  gaussian,  real  random  variables  as  eqs.  (15)  and 
(16). 


E  [x,  X, 
E  [xj  X^ 


-  -  ^^!<]  -  I 

All  pairs 


TT  E  fx.x. 

TT  L  ^  J 


(15) 


(16) 


Equation  (16)  merely  states  that  the  expectation  of  the  product  of  an 
odd  number  of  random  variables  is  zero.  In  this  analysis,  the  number 
of  random  variables  in  -ny  product  does  not  exceed  four.  In  this  case, 
eq.  (15)  (evaluated  for  k=2)  may  be  expressed  as  (17). 


E  [X,X3X3XJ  =  E  [X^XJ  E  [X3XJ  .  E  [X.X3]  E  [x,xj 

+  E  [xjXj  E  [X2X3]  ,17) 

Equation  (17)  will  be  employed  in  the  expansion  of  (13)  and  (14)  using 
(5)  and  (6)  in  the  wide  sense  stationary  case,  i.  e.  tj  =  t  and  t^  ^t  +  t. 
The  first  step  in  the  expansion  is  shown  as  (18)  and  (19). 


*e> 
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Rc(t)  =  E  y  +  2  (M^.Y  -  M^X) 


Y^  -  X^  +  2  (M  Y  -  M  X) 
y  X 


t+T  J 


R_  (t)  =  E  4  XY  +  M  X  +  M  Y 
®  y  X 


XY  +  M  X  +  M  Y| 

y  X 


The  truncated  expansiona  o£  eqs.  (18)  and  (19)  are  shown  as  eqe.  (20) 

and  (21)  where  the  expectations  of  all  odd-numbered  products  have  been 
eliminated  on  the  basis  of  (16). 


R,(rl  =  E  .  X^),(y2  . 


+  4  E  (MyY  -  M^X)^  (M  Y  -  M^X) 


X  't+T 


R  (t)  =  4E  +  M 


[(XV)t  (XY),^^]  t  [x,  X,^^] 


K  ’'ttr] 


+  M  M  L  v'  Y 

X  y  I  t  t+T 


"  ^t  ^ttr] 


Exp3:,aion  of  eqs.  (20)  and  (21),  using  (17)  where  applicable,  yield, 
eqs.  (22)  and  (23)  in  convariance  notation. 


42 


R^(t)  =  2  [c^(T)  +  C^(T)  .  C^y(T)  .  C‘^(T)j 

*  4  [MjCy(T)  +  -  M,,My  |c^y(T)  +  C^(T)|] 

+  rc^(o)  -  2  Cjj(o)  C^(o)  +  C^(o)]  (22) 

R^(t)  =  4  [c^(T)  C^(T)  4  C^^(r)  C^^(T)] 

4  4  [m^  C^(t)  4  Cy(T)  4  My  |c^y(T)  4  Cy^(T)|] 

+  4  (o)  (23) 

xy 

Observe  that  the  last  terms  of  eqs,  (2Z)  and  (23)  are  simply  the 
squares  of  the  means  given  by  (9)  and  (10).  Thus,  in  accordance  with 
(11)  and  (12,  the  covariance  functions  of  the  input  rotational  field  errors 
are  expressed  as  eqs.  (24)  and  (25). 

C^(4)  =  2  [cl(r)  4  0^(4)  .  C^(4)  .  C^(4)] 

4  4  [m^  Cy(T)  4  M^  C^(T)  -  My  (c^y(T)  4  Cy^(T))] 

(24) 

C^(T)  =  4  [c^(4)  Cy(4)  4  C^y(T)  Cy^(4)l 

4  4  [m^  C^(t)  4  mJ  Cy(4)  4  M^  My  |c^(T)  4  Cy^(T)|] 

(25) 
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d.  Output  Error  Variances 

The  variances  of  the  rotational  field  errors  at  the  output 
of  the  RGG  signal  process  rr.ay  be  obtained  by  integration  of  the  output 
error  spectra  of  each  ch-^nnel  of  the  RGG  as  defined  by  eqs.  (26)  and 
(27). 


+  0O 


CO 


■/ 


(f)  df 

CO 


(26) 


so 


+  CO 

■/ 

-  CO 


S  (f)  df 
so  ' 


(27) 


The  output  error  spectra  may  be  obtained  from  the  input  error  spectra 
and  the  equivalent  filter  of  the  signal  process  as  shown  in  eqs.  (28) 
and  (29). 


CO 


so 


1%  « 

|h  (j2,  f) 

I  *=4 


'  (f) 

"  S,  (£) 


(28) 

(29) 


The  input  error  spectra  may  be  derived  in  terms  of  the  angular 
rate  spectra  by  taking  the  Fourier  transforms  of  the  covariance  func¬ 
tions  of  (24)  and  (25).  The  general  results  of  this  process  are  expressed 
as  eqs.  (30)  and  (31)  wherein  convolution  operations  on  the  rate  spectra 
are  required  for  the  nonlinear  terms  in  the  covariance  functions. 
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-  2  I*  S  (f)  ’i'  S  (f)  +  S  (f)  S  (f)l 

L  xy'  xy'  yx'  '  yx'  'J 

"4[M^Sy(f)  tM^s^lt)] 

-  4  [®xy<''  ^  Sy^(o]  (30) 

S  (£)  =  4  fs  (f)  S  (f)  +  S  (f)  *  S  (f)l 

s  1_  X  y'  '  xy'  yx'  'J 

+  4  I  S  (f)  +  S  (f)l 
1_  y  X  X  y  J 

+  4MM  Fs  {f)+S  (f)l  (31) 

X  y  I  xy'  '  yx  J  '  ' 

An  approximate  method  of  determining  the  output  variances  is 
to  replace  the  output  spectra  in  eqs.  (26)  and  (27)  with  the  input  spectra 
and  to  restrict  the  integration  to  a  narrow  frequency  band  centered  at 
zero  frequency  of  width  equal  to  the  equivalent  bandwidth  of  the  RGG 

1  (33). 

S^(f)  df  (32) 

S^(f)  df  (33) 


signal  process  as  shown  in  eqs.  (32)  and 


L 


Af 

T 


+  - 


Af 
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Summary  and  Conclusions 

The  covarian-e  functions  of  the  rotational  field  errors 
were  derived  for  the  case  where  the  random  angular  rates  of  the  RGG 
measurement  frame  have  wide  sense  stationary,  jointly  gaussian 
probability  density  functions.  Also  a  method  for  determining  the  vari 

ances  of  the  rotational  field  output  errors  of  each  channel  of  the  RGG 
was  presented. 

Application  of  the  results  of  this  study  presents  computational 
difficulties  because  of  the  requirement  for  spectral  convolution  of  the 
angular  rate  spectra  and  because  of  the  requirement  for  knowledge  of 

the  cross  spectra  in  addition  to  the  co-spectra  of  the  RGG  angular 
rates . 
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3. 


Tran„ver„^  Differentt.l 


The  component  ol  mans  nnbalance  of  each  arm  of  the  RGG 
.tonal  accelcrattone  producer  time  varying 

prtnc.pai  tore.onaI  axis  of  the  RGG.  The  difference  of  theee  arm- 
mass-unbalance  "body  torques"  nf  f, 

of  fV,  >^orques  of  the  two  arms  is  a  driving  function 

_  the  differential  mode  of  the  RGG  and,  as  such,  is  capable  of  produc 
rrors  in  t  e  sensible  output  of  signals  of  the  RGG.  The  intent  of 
IS  memorandum  is  to  derive  analytical  expressions  which  relate 
these  output  errors  to  the  input  translational  accelerations  and  ^  the 
ansverse  differential  arm-mass-unbalance  of  the  RGG. 

■^.ief  Descriptions  of  Arm-Mass-Unbalano. 


general,  the  acceleration  induced  "body/torque" 

actmg  .  ,  ^  ^  centSf^mas,  is 

r  s  t  a  ;  r 

body  -PPor.  and  m  is  dte  mass  of  the 


L  =  m  a  X  r 

d 


ne  "body  torque"  of  each  arm  of  the  RGG  achng  about  its  principal 
Stonal  axts  (z)  may  be  stated  from  eq.  (1)  as  (2)  and  (3). 


T,  =  T 


•  h  n  'ij 


T,  =  z  .  [m^  a^  X  r  J 
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In  the  case  of  pure  translational  acceleration  of  the  RGG,  the  "center 
of  support"  accelerations  are  equal,  and  the  effective  torque  acting 
on  the  differential  mode  of  the  RGG  may  be  expressed  as  the  difference 
of  eqs.  (2)  and  (3). 

To  the  first  order,  the  RGG  arm  masses  are  equal,  and  the 
differential  mode  unbalance  driving  torque  may  be  expressed 

as  eq.  (4). 


T 


d 


(4) 


where 


h  = 


Expansion  of  eq.  (4)  in  RGG  cease-fixed  coordinates  (xyz)  yields 
eq.  (5). 


T 


d 


m  (a  h 
X  y 


a  h  ) 

y  X 


(5) 


The  differential  arm  mass  unbalance  components,  h^  and  h^, 
are  time  variable  due  to  the  spatial  modulation  produced  by  the  RGG 
rotation  frequency  (w^),  i*  e. 


h  (t)  =  h  cos  (w  t  +  4> ) 

X  s 

h  (t)  =  h  sin  (u)  t  +  <t> ) 

y  ® 

where  <j)  is  the  fixed  mechanical  phase  of  the  differential  arm-mass- 
unbalance  vector  (h)  relative  to  the  principal  axes  of  the  RGG  arms. 
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For  analytical  convenience  here,  this  mechanical  phase  will  be 
assumed  zero  such  that  eq.  (5)  may  be  expressed  as  eq.  (8). 


mh 


cos  CO  t 

s 


a 

y 


sin  CO  t 
s 


(8) 


At  this  point  it  is  convenient  to  normalize  the  error  torque  of 
eq.  (8)  with  respect  to  the  input  signal  (2  F.y  of  the  RGG  signal  process. 
This  normalized  error  function  is  shown  as  eq.  (9),  where  q  is  the  arm 
inertial  efficiency  and  C  is  the  principal  polar  inertia  of  one  arm. 


^d 


qC 


a  a 

X  .  V  . 

—  cos  CO  t - ^  sin  (o  t 

g  s  g  s 


(9) 


b*  Signal  Processing 

For  the  purposes  of  this  analysis  it  is  convenient  to 
determine  the  power  spectrum  of  the  normalized  error  function  of 

eq.  (9).  This  may  be  accomplished  by  first  deriving  its  autocorrela¬ 
tion  function  as  defined  by  eq.  (10). 

* 

Vr‘2'  "  =  [{'J  f'l' Vz*]  ('O' 


Assuming  the  statistics  of  the  input  accelerations  are  wide  sense 
stationary,  the  error  function  autocorrelation  may  be  expressed  as 
eq.  (11). 


Rd(T)  = 


mngl 

qCj 


1/2  (r  (t)  +  R  (t)1  COS  CJ  T 

'  X  y  *  s 

+  1/2  (t)  -  I 


(11) 
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When  the  input  accelerations  are  uncorrelated,  their  cross-correlations 
are  zero,  and  eq.  (11)  simplifies  to  eq.  (12). 


Rd(T) 


=  1/2 


R,(T) 


V^'] 


cos 


U)  T 

s 


(12) 


The  Fourier  transform  of  eq.  (12)  yields  the  power  spectrum  of  the 
error  function  for  the  assumed  conditions. 


If  we  assume  that  the  spectra  of  the  accelerations  on  each  axis  are 
identical,  eq.  (13)  simplifies  to  eq.  (14). 


S,(f) 


2 

= 

\(f  +  f^)  +  SJf  .  ij 

L  c  J 

2 

(14) 


Equation  (14)  shows  that  the  error  function  spectrum  is  proportional 
to  the  frequency  shifted  acceleration  spectra.  This  spectrum  is 
filtered  by  the  narrow  band  process,  Hj(s),  of  the  RGG  defined  by 
eq.  (15). 


where 


Hj(s)  = 


u 


w  =  2w 
o  s 


(15) 
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The  resultant  spectrum  after  the  narrow  band  filter  process  is  stated 
as  eq.  (16). 


S^(f)  =  |Hj(j2Tf)|  2  S^(f)  (16) 

After  the  narrow  band  filter  process,  the  filtered  error  signal  is 

demodulated  into  two  channels  by  the  linear  functions,  2  cos  w  t  and 

o 

2  sin  w^t,  to  yield  the  spectra  of  eqs.  (!7)  and  (18). 


S  (f) 
c 


S^(f) 


Se(f  +  y  +  S^(f  .f^)J 

Se(f  +  V  +  S^(f  -  f^)  j 


(17) 


(18) 


After  demodulation,  the  signal  in  each  channel  is  filtered  by  the  low 
pass  filter,  H^(s),  to  remove  the  second  harmonics  of  (f  ).  The 
resultant  spectrum  of  each  channel  is  stated  by  eq,  (19). 


S 


'I' 


Ho(j2 


s,(f- 


(19) 


The  variances  of  the  output  errors  may  be  obtained  by  integration  of 
eq.  (19)  over  all  frequencies  as  in  eq.  (20). 

Ho(j2Trf)|  j^S^(f  +  f^)  +  s^(f  -  f^)j  df  (20) 
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This  integral  may  be  evaluated  to  sufficient  ac 

the  previously  described  processes  graphically  T’'  ^ 

spectrum  of  eq.  (16).  This  is  sketched  as  fJ' 6 

ts  determined  from  eq.  ,14,  as  stated  in  eq  (gl|  '  ^ 


) 

e  s 


■(■S') 


Sa(3f  1  +  S  (f  ) 


The  demodulation  processes  of  eos.  (I  7,  .  a  „  c, 

fig.  6  by  ±  f  .  xiji  ,  "  )  shift  the  spectrum  of 

O  result  IS  shown  as  Fig.  7. 

The  low  pass  filter  R 

fo  ftom  the  spectrum  of  Fig.  °  end  °f 

frequency  centered  part  as  shown  in  Fig.°7  '  ^ero 

The  equivalent  bandwidth  of  the  spectrum  „f  F-  o  • 

"lately  the  same  as  the  equivalerf  •  u  8  is  approxi- 

-.-process,  xhus. 

terms  of  an  -effective  acceleration-  variance  (a  ^ 

Fig.  8  as  eq.  (22).  '  a  ^  eq.  (21)  and 


where 


O'  ^  =  0-  2  ~  ( mhg  1  ^ 

s  ~  y 
a"  “  K-fs'  +  S^(fe)J  "eq 


Care  should  oe  observed  in  evaluating  the  "effect- 
variance  to  ensure  that  the  "two-sided"  m  ■  ^"“l«etion" 

spectra  are  used  and  that  thea  egmhldes  of  the  acceleration 

of  "g^/Ha.  "  magnitudes  are  expressed  in  the  units 
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The  standard  deviations  of  the  gradient  tensor  elements  may¬ 
be  determined  from  eq.  (22)  with  one  exception,  i.  e. ,  the  trace 
elements  of  the  gradient  tensor  are  obtained  from  pairs  of  RGG  out¬ 
puts*  Common  accelerations  acting  on  the  mass  unbalance  of  each 
RGG  produce  correlated  errors  in  the  gradient  tensor  trace  elements 
depending  on  the  relative  mechanical  phases  of  pairs  of  mass  unbalance 
vectors.  The  resultant  trace  element  variance  may  be  stated  as 
eq.  (23). 


2 
(T  . 

11 


1/9 


h. 


hih2 


cos  ((j>2  -  «j>j 


(23) 


The  subscripts  are  used  to  denote  the  parameters  of  each  RGG  in  a 
particular  pair.  The  cross-gradient  variances  are  a  function  of  one 
instrument  only,  and  these  variances  are  defined  by  eq.  (24). 


(24) 


To  avoid  confusion,  the  variances  of  all  of  the  gradient  tensor  elements 
are  expressed  as  eqs.  (25  through  (30)  using  the  spin  axis  coordinate 
as  the  parametric  identification. 


I  h  h  ^  +  h  h 

L  y  z  y  Z 

h^  +  h  ^  +  h  h 

X  z  X  z 

jh  ^  +  h  ^  +  h  h 

L  X  y  X  y 


cos  ((j> 


y 


cos  (4> 


X 


cos  (6 

X 


(25) 


(26) 


(27) 
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^  j 

't]  ' 


Note  that  the  mass  unbalance  parameters  (h.)  are  the  total  magnitu  Jes 
of  the  mas^s  unbalance  vectors  of  each  RGG  whose  spin  axis  coincides 
with  the  (i  )  reference  coordinate. 

Error  Variance  Evaluation 

The  first  step  in  evaluating  the  error  variances  is  the 
evaluation  of  the  "effective  acceleration"  variance  defined  in  eq.  (22). 
Previous  estimates  of  the  acceleration  spectra  on  a  "one-sided'  basis 
yield  the  following  values  at  one  and  three  times  the  spin  frequency. 

2S  (f  I  5  4.5x10-’8^/Hz  ,3,, 


2S^(3tJ  3  8.9x10-‘%^/Hz 

The  "effective  acceleration"  variance,  based  on  an  equivalent  band- 
width  of  =  0.05  Hz,  is  stated  as  eq.  (33). 


j  2  «  4.  5  X  10'^  t 
a  “  2 


8.9  X  10'^^  ['^^J  “ 


11. 6  X  10"  g  (33) 


The  standard  deviations  of  the  gradient  tensor  elements  will  be  evaluated 
assuming-  that  all  mass  unbalance  magnitudes  are  equal  and  that  the 
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correlating  mechanical  phase  aneles  are  zem  tv.  o 

^  angles  are  zero,  Th,s  assumption  yields 

e  worst  case,  and  the  results  are  stated  ae  eqs.  (34)  and  (35). 


nC 


(34) 


lie 


(35) 


Summary  and  Conclusions 

The  RGG  output  error  due  to  differential  arm-niass- 
unbalance  driven  by  translational  acceleration  was  investigated  in 
tern-.s  of  the  power  spectra  of  the  input  accelerations  and  the  filter 
process  of  the  gradiometer.  The  results  show  that  this  RGG  error 
variance  is  proportional  to  input  acceleration  spectral  magnitudes  at 
one  and  three  times  rotation  frequency  (as  anticipated). 

An  "effective  acceleration"  variance  was  defined  for  this  error 
source  in  terms  of  the  acceleration  spectral  magnitudes  at  one  and 

three  times  rotation  frequency  and  of  the  equivalent  bandwidth  of  the 
RGG  signal  process. 

Analytical  expressions  were  derived  for  the  gradient  tensor 
error  variances,  and  it  was  shown  that  the  trace  element  variances 
are  functions  of  the  phase  differences  of  the  mass  unbalance  vectors 

of  the  pairs  of  instruments  whose  outputs  are  used  to  compute  the 
particular  trace  elements. 


RQC  Torques  Due  to  Paramagnetic  Met.ri.l. 

Most  m,iterials  suitable  for  use  in  the  construction  of  the  RGG 
are  paramagnetic  and  magnetic  fields  can  introduce  torques  in  these 
paramagnetic  materials  at  the  even  harmonics  of  the  rotation 
frequency.  (The  harmonics  higher  than  the  second  will  be  rejected  by 
the  signal  processor  and  will  be  disregarded  in  this  analysis. ) 
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If  a  paramagnetic  bar  is  introduced  into  a  uniform  magnetic 
field  magnetic  poles  will  be  induced  in  the  bar  and  the  bar  will  tend  to 
align  Itself  with  that  field.  If  it  is  assumed  that  all  of  the  flux  from 
the  induced  poles  emanates  from  the  end  of  the  bar  then  estimates  can 
be  made  of  the  induced  torque  on  the  bar.  This  is  an  oversimplification 
but  it  provides  simple  equations  with  practically  acceptable  accuracies 
for  relatively  long  bars  and  the  assumption  will  be  used  in  this  analysis. 


pole  is 


Magmtic  Fields  and  PoIpb 


In  the  ogs  electromagnetic  system  the  force  on  a  magnetic 


F  =  Hm 


where: 


F  =  force  on  pole  in  dynes 
H  =  field  strength  in  oersteds 
m  =  pole  strength  in  emu. 


The  flux  emanating  from  a  magnetic  pole  is 

<}>  =  4Trm  (maxwells)  . 

If  a  paramagnetic  bar  is  introduced  into  a  magnetic  field  the 
magnetic  flux  balance  in  the  bar  will  be 


*b  ■  *o  = 
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where: 


<|)^  =  flux  in  the  bar 

<|)^  =  the  flux  in  the  volume  occupied  by  the  bar 
with  the  bar  removed 

4'n'm  =  flux  due  to  the  induced  poles  in  the  bar. 
Dividing  (3)  by  the  cross  section  area  of  the  bar,  A  we  obtain 


*^b  *^o  _  4TTm 

A  ■  A  ^  A 


(4) 


where 


and 


and 


<t>^/A  —  =  flvix  density  in  the  bar 


(t>^/A  —  =  flux  density  in  free  space 


B^  —  H  in  the  emu  system. 


Rewriting  (4)  using  these  definitions,  we  obtain 


H  + 


4Trm 

A 


(5) 


Thus,  the  magnetic  flux  density  in  the  bar  "lan  be  calculated  if 
the  field  strength  and  the  induced  pole  strength  are  known. 
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b.  Magnetic  Field  Due  to  a  Current  Carrying  Conductor 

The  magnetic  field  strength  due  to  current  flow  in  an 
infinitely  long  conductor  is,  at  a  distance  x; 

2^ab 

H  =  —  oersteds  (6a) 


and  the  magnetic  field  gradient  at  thif  point  is: 


dH 

dx 


2  lab 

- 2“  oersteds/cm 


(6b) 


where: 


lab  “  abamperes  =  10  amperes  practical 
X  =  distance  in  centimeters 

For  short  lengths  of  wire,  the  above  calculations  are  accurate,  to  a 
few  percent  if  the  length  is  10  times  the  distance  x. 

c.  Permeability  and  Susceptibility 

These  tern-is  are  used  to  describe  the  magnetic 
characteristics  of  materials.  Permeability  (p)  is  defined  in  the  emu 
system  by  (7) 


b  _  ,  4'n'm 

H  "  ^  Ha 


(7) 


It  is  seen  that  the  right  hand  side  of  (7)  can  be  obtained  by  dividing  (5) 
by  H. 
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Susceptibility  (X)  is  defined  by  (8)  for  susceptibility  per  unit 

mass. 

X  =  K/p  (8) 


Susceptibility  may  also  be  defined  on  the  basis  of  molecular  weight 
(Xm)  atomic  weight  (x^).  This  section  will  be  based  on  the  unit 
mass  values  commonly  tabulated  for  metals. 

In  equation  (8)  p  is  the  density  and  K  is  a  constant  defined 


Thus,  equation  (7)  can  be  written 


|i  =  1  +  4ttK 


(9a) 


or 


p  =  1  +  4TrXp  (9b) 

The  permeability  and  susceptibility  of  some  cqi^mon  materials 
are  shown  in  Table  II  from  Refs.  1,  3,  and  5. 

d.  Torque  on  a  Bar  Magnet 

If  a  bar  magnet  length  i,  with  pole  strengths  +m  and  -m, 
is  introduced  into  a  uniform  field  of  strength  H  at  an  angle  0  to  the  field, 
we  have  from  eq.  (1)  and  the  geometry 

Torque  =  mfH  sin  0  .  (10) 
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TABLE  II 


Magnetic  Properties  of  Pure  Materials 


X 

X  10^ 

P 

P 

Aluminum 

+  0.  63 

2.  68 

1. 0000212 

Beryllium 

-1.0 

1.  85 

0. 9999767 

Copper 

-0. 086 

8.  96 

0.  9999900 

Tungsten 

+0.  28 

19.  3 

1. 0000675 

Uranium 

+  2.  6 

19.  1 

1. 000624 

Elgiloy 

1. 00004 

Air 

1. 00000036 

Mercury 

-0.  18 

13.  6 

0. 999969 

Crown  Glass 

-0.  90 

2.  47 

0.  999972 

Permalloy 

100,  000 

Iron,  technically  pure 

6,  500 

Cast  iron 

600 

T724 


Linear  Magnetic  Field  Torque  on  a  Paramagnetic  Bar 

jfgff 

It  is  seen  from  eq.  (3) 

■t'b  -  *0  =  ( 

that  poles  of  strength  are  induced  in  a  bar  when  it  is  introduced 
into  a  field.  F rom  (7)  for  a  field  of  H  oersteds 


or 


4Trm£ 

HA£ 


1  + 

H  Vol 


or 


H(ji  -  1) 


4'iTm^ 

Vol 


!? 


and 


ml  =  <l^-‘>HVol 
4it 


(11) 


This  magnet"  of  pole  strength  m  and  length  H  can  be  substituted 
into  (10)  and 


Torque  =  (ij-.-D  Vol  sin  6 
^  4'it 


(12) 


Two  things  have  not  been  taken  into  account  in  the  derivation  o'  (12). 
These  are:  the  demagnetizing  effect  and  the  reduced  projected  area  of 
the  ends  of  the  bar  when  it  is  not  parallel  to  the  field.  The  reduced 
area  coefficient  is  simply  cos  6. 

When  the  bar  is  introduced  into  the  field  the  induced  poles  tend 
to  counteract  the  initial  field.  Thus,  the  effective  field  H  is 

H  =  H'  -  NKH  =  H'  -  NxpH  (13 


where 


H  =  effective  field 
H'  =  initial  field  in  a  vacuum 
N  =  demagnetizing  factor. 
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N  varies  from  nearly  zero  for  long  thin  rods  to  4^  for  flat  plates.  For 
paramagnetic  materials  with  mass  susceptibilities  in  the  ranges  shown 
in  Table  II,  we  can  say 


H  «  H' 

even  if  N  is  as  large  as  4tt,  its  maximum  value. 

Thus,  the  final  torque  equation  for  a  paramagnetic  bar  intro¬ 
duced  into  a  magnetic  field  of  strength  H  is 


(u-1)  H  Vol  sin  0  cos  0 
Torque  =  - 


(14) 


or 


Torque  =  XP^  ®  ® 


(15) 


where 


Torque 

X 

P 

H 

Vol 


dyne  cm 
emu  units 
grams/cc 

field  strength  in  oersteds 
cubic  cm  of  the  bar 


f.  Forces  on  a  Paramagnetic  Mass  in  Nonuniform  Field 

The  force  on  a  homogenous  mass  M  in  a  nonuniform 
magnetic  field  is  derived  by  Williams  (ref.  1): 


dH 

f  =  )(MH  dx 


(16) 
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where 


f  =  dynes 
X  =  emn 
M  =  mass 
H  =  oersteds 
dH/dx  =  oersteds/cm  . 

The  torque  on  an  arm  with  a  concentrated  mass  M/2  at  each  end 
and  an  arm  length  of  2r  in  a  gradient  field  is 


but 


Torque  =  =  (fj-f2)rcose 


H 


2 


H, 


r  sin  0 


dH 

dx 


for 


constant 


dH 

dx 


and 


Torque  =  y  ^  r  cos  0  j  +  r  sin  0  ^  j 


Y  M  2  dH  ^  „ 

X  y  r  ^  sin  0  cos  0  (18a) 
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This  can  be  written  from.  (9b)  as: 


Torque  = 

Stt  / 


(18b) 


where 


fx  =  permeability  of  the  end  masses 

c 

^m  ~  volume  of  both  end  masses 

r  =  radius  to  the  center  of  the  end  masses 
dH 

dx  "  oersteds/cm  gradient 


Summary  of  Magnetic  Field  Induced  RCG  Torques 

Torques  «re  introduced  in  the  RGG  due  to  formation  of  magnetic 
poles  in  the  arms  and  end  masses  and  by  the  magnetic  field  gradient 
acting  on  the  end  masses.  The  gradient  torque  on  the  arm  itself  is 

negligible  compared  to  the  torque  due  to  the  end  masses.  These  torques 
are  from  eqs.  (14)  and  (18) 

(K  -  1)  V 

Torque  - - - L  ^  ^ 


where: 


+ 


-  1)  V 


m 


Sir 


2 

r 


sin  6  cos 


6 


(19) 


l^a  -  a  composite  value  of  permeability  for  the  arm  and 
end  mass  as  a  whole  =1.  00005 

Va  =  volume  of  arm  and  end  masses  160 
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=  permeability  of  the  end  masses  si.  00007 
r  =  radius  from  pivot  to  center  of  the  end  masses  a;5  cm 
~  volume  of  end  masses  sSOcc 

Inserting  the  above  approximate  values  the  peak  torque  equation  becomes 


Torque 


3.  1  X  lO"'^  +  2.  8  X  10' 


cm  peak 


(20) 


h.  Magnetic  Fields  and  Gradients  Expected 


The  magnetic  field  of  the  earth  has  a  field  strength  of  about 
0.  7  oersted  and  a  maximum  gradient  of  approximately 
3.  3  X  10  ^  oersted/cm.  The  magnetic  field  of  a  long  straight  wire 
carrying  one  ampere  is  0.  02  oersted  at  10  cm  and  the  gradient  at  this 
point  is  0.  002  oersted,  .it  is  believed  that  all  magnetic  fields,  except 
that  due  to  the  earth,  can  be  kept  below  about  0.  1  oersted  and  the 
gradients  below  0.  01  oersted/cm.  Thus  it  is  evident  that  only  the 
earth's  field  is  significant. 

The  RGG  torque  introduced  by  the  earth's  field  could  be  as 
large  as: 


Peak  Torque  =  3.  1  x  lO"^  x  (0.  7)^  +  2.  8  x  lO’^  x  (0.  01)^ 


_4 

=  1.  5  X  10  dyne  cm  (due  to  earth's  magnetic 

field)  (21) 


The  peak  torque  in  the  RGG  due  to  a  gravity  gradient  is 


Peak  Torque  =  qC 


=  1.  5  X  10  d-cm/EU 
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Thus  the  peak  torque  at  the  sensor  resonant  frequency  due  to  the  earths 
magnetic  field  will  be  approximately  10  EU.  However,  this  is  constant 
in  direction  and  magnitude  over  distances  of  hundreds  of  miles  except 
in  areas  with  iron  or  nickel  ore  concentrations.  Thus,  the  induced 
magnetic  torques  can  be  considered  as  a  bias  and  compensated  with 
trivial  error. 


i.  Conclusions 


The  earth's  magnetic  field  will  introduce  significant  bias  terms 
in  the  RGG  gradient  tensor  outputs.  However,  these  will  be  constant 
over  long  periods  of  time  and  over  large  distances  and  can  be  com¬ 
pensated  with  negligible  error. 

Torques  due  to  nearby  magnetic  devices  and  current  carrying 
conductors  can  be  held  to  trivial  values.  Shielding  of  tne  RGG  will  not 
be  necessary  in  general. 

C.  DESIGN  CONCEPTS  AND  PARAMETERS 

1.  Selection  of  resonant  frequency  and  Q  for  the  RGG 

2.  RGG  Sensor  Electronics 

3.  Arm  Mass  Unbalance  and  Arm  Inertia 
Unbalance  Mechanical  Adjustments 

4.  Vibrational  Motion  Environment  Estimates 

5.  Preliminary  Gradiometer  Designs 

6.  Isoelasticity  of  Spherical  Gimbal  Support  with 

Hydros..atic  Bearing  Pads 

7.  Spin  Bearing  Requirements 

a.  Gradient  Errors  Due  to  Sum-Mode  Mismatch 

b.  Torque  Variations  with  Eccentricity  in 
Fluid  Journal  Bearing 
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c.  Phase  Error  Due  to  Spin  Bearing 
Disturbance  Torque 

d.  Spin  Bearing  Procurement  Specification 

8.  RGG  Baseline  Parameters 

These  sections  treat  many  aspects  of  the  RGG  design  and 
establish  values  for  some  of  the  parameters.  Future  design  and 
analysis  will  be  based  on  these  parameters. 

1.  Selection  of  R(esonant  Frequency  and  Q  for  a  RGG 


a.  Introduction 

Inspection  of  the  factors  involved  in  the  selection  of 
resonant  frequency  and  Q  demonstrates  that  a  single  equation 

cannot  be  written  that  will  provide  a  precise  determination  of  the 
optimum  values  for  these  parameters.  Instead  each  factor  must  be 
evaluated  in  terms  of  its  effect  on  a  complete  gravity  gradient 
measuring  system  and  compromises  made  as  required.  The  factors 
to  be  considered  are; 

(1)  The  thermal  noise  in  the  gravity  gradient  tensor 
elements  has  a  specific  limit. 

(2)  The  transducer  output  signal  level  must  be  at 
least  a  few  lO's  of  nanovolts  per  EU  signal. 

(3)  For  ease  in  fabricating  low  noise  electronics 

the  signal  frequency  should  be  above  appioximately 
10  Hz. 

(4)  The  signal  phase  shift  uncertainty  must  not 
exceed  0.  735  m^J  (0.  0042  deg)  total  through 
the  system. 

(5)  The  sensor  time  constant  (  t  =  2Q/w^)  cannot  be 
allowed  to  be  equal  to  the  total  allowed  system 
time  constant  since  some  filtering  must  be 
provided  for  preamplifier  noise. 
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(6)  The  frequency  of  the  sensor  must  be  high  enough 
that  the  mechanical  parts,  particularly  pivots, 
are  not  too  flimsy  to  work  with, 

(7)  The  piezo-electric  transducer  output  impedance 
must  be  reasonable,  i.  e.  not  too  high. 

(8)  The  uq  selected  should  avoid  any  peaks  in  the 
environmental  spectrum. 

(9)  The  natural  frequency  must  not  be  so  high  that 
the  losses  in  the  spin  bearing  become  intolerable. 

(10)  The  centrifugal  acceleration  on  the  arm  mass 
balance  adjustment  tubes  must  not  exceed 
approximately  50  g's. 

b.  Baseline  Parameter  Selection 

After  reviewing  these  factors  it  was  decided  to  make  a 
preliminary  selection  of  the  limiting  sensor  and  system  parameters 
and  determine  if  this  allowed  a  viable  detailed  baseline  RGG  sensor 
and  system  design.  Accordingly  the  following  parameter  limits 
were  established. 

c.  Preliminary  Baseline  Parameter  Limits 

(r—  =  1/3  EU  maximum  RMS  noise  signal 

for  each  of  the  elements  of  the  compmted 
gravity  gradient  tensor 

Q  =  300  max.  The  Q  of  the  sensor  difference 

mode. 

f  =35  Hz  nominal  undamped  resonant  frequency 
of  the  sensor  in  the  difference  mode. 

T.  =  10  sec,  system  (sensor  plus  filters) 

equivalent  integration  time  defined  by  the 
AFCRL  specification. 

The  rationale  for  the  selection  of  each  parameter  is  discussed 
in  the  various  sections. 
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(1)  Thermal  Noise  Cons ide rations  —  It  has  recently  been 

shown  that  the  one  sigma  value  of  the  noise  output  of  the  sensor,  a 
demodulate r-RC  filter  system  and  a  RSS  sine  and  cosine  summing 
circuit  is 


4 

n  v/^ 


(1) 


where: 


T  =  the  sensor  time  constant 
s 

=  equivalent  integration  time  of  the  sensor  and  the 
filter;  (t^  “  ^  ^f^  which  is  specified  as  10  sec¬ 

onds  for  this  system. 

q  =  inertia  efficiency  of  the  sensor  arms. 


C^_  =  polar  moment  of  inertia  of  one  sensor  arm. 
z  z 

k  =  Boltzman  constant 

T  =  temperature  °K  (290°K  =  IT^’c  assumed) 


The  above  equation  can  also  be  vrritten 


‘"ij 

II 

/  2kT 
'  ^zz  ^ 

1^  „  T. 

y  2q 

and 

O'.  . 

I'll 

IZkT 

/CJ  T. 

1  /  °  ^ 

11 

n3T. 

V  Czz 

V  2Q 
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where 


<r. . 
ij 


<r.. 

11 


=  one  sigma  value  of  thermal  noise  in  the  off  diagonal 
elements  of  the  gradient  tensor.  Specified  in  the 
preliminary  parameters  as  1/3  EU  maximum. 

=  one  sigma  noise  in  the  trace  elements. 


Since  <r^j  is  slightly  larger  than  and  if  (t.j  is  held  to  1/3  EU,  the 
system  will  be  within  the  required  limits  on  the  other  elements  of 
the  tensor. 

The  thermal  noise  analysis  also  shows  that: 


N 


where; 


0"^^  =  one  sigma  of  the  thermal  noise  in  the  root-sum- 
square  sensor  output. 

cr^  =  one  sigma  of  the  thermal  noise  in  the  sine  channel. 

=  one  sigma  of  the  thermal  noise  in  the  cosine  channel. 


Thus,  with  a  1/3  EU  system  limit  on  the  limiting  value  of  is 
0.  946  EU  due  to  thermal  noise. 

The  parameters  in  eq.  (2)  can  be  considered  as  follows.  It  is 

always  a  design  goal  to  make  the  sensor  inertia  efficiency  q  as  high 

as  possible  with  practical  values  from  0.  7  to  0.  97.  The  required 

value  for  o',,  has  just  been  shown  to  be  0.  946.  The  sensor  must  be 

operated  near  300°K  and  the  integration  time  i  .  is  specified  by  the 

AFCRL  contract.  Since  and  Q  have  been  sel-jcted  on  a  preliminary 

basis  considering  all  of  the  pertinent  factors  the  only  remaining 

factor  in  eq.  (2)  is  C 

^  zz 
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Evaluating  cr.j  =1/3  EU,  in  eq.  (2)  we  find 


fo  r 


T. 

1 


T 

s 


r^C 


■f 


zz 


10  sec  specified  =  Tg  + 

2Q/u  =  2.  73  seconds 
o 

7.  27  sec 

26,  400  gram  cm^  (2.  64  x  lO”^  kg  m^) 


f  =  35  Hz 
o 

=  220  rad/sec  (actual  219.91) 
kT  =  4x10^^  ergs  (4  x  lO”^^  joules) 

T  =  290°K  (17°C) 

Q  =  300 

All  of  these  values  seem  to  be  within  practical  design  limits 
and  the  thernal  noise  consideration  is  satisfied. 

(2)  Signal  Level  Considerations  —  The  signal  level  must  be 
at  least  a  few  lO's  of  nanovolts /EU.  Signal  levels  below  10  nanovolts 
at  frequencies  below  10  kHz  are  very  difficult  to  work  with.  They 
are  difficult  to  generate  for  test,  difficult  to  shield  and  difficult  to 
amplify  linearly.  In  addition,  the  one  EU  signal  level  .should  be  much 
higher  than  the  preamplifier  noise  so  that  the  preamplifiv'ir  noise 
does  not  contribute  significantly  to  the  output  signal  noise.  However, 
if  the  signal  level  exceeds  a  few  hundred  nv/EU  littlo  is  gained  by 
forcing  the  pie  zoelectric  transducer  design  to  attain  the  highest  possible 
signal  level  since  this  level  is  limited  to  about  1000  nv/E  J  by  other 
considerations.  A  signal  level  of  about  1000  nv/EU  requires  about 
the  same  circuitry  and  the  same  care  as  a  signal  level  of  a  few 
hundred  nv/EU.  A  final  reason  for  not  attempting  to  obtain  a  very 
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large  signal  out  of  the  piezoelectric  material  is  that  this  would  be 
a  deviation  from  Hughes  experience.  Hughes  has  not  studied  or 
experimented  with  transducer  characteristics  at  high  signal  levels. 

The  peak  differential  angle  of  the  sensor,  from  page  380  of  the 
AFCRL  proposal, is 


e 


peak 


(3) 


Since  the  output  signal  is  a  linear  function  of  differential  angle  (all 
other  parameters  remaining  constant)  it  is  apparent  that  the  largest 
output  voltage  would  be  obtained  by  making  Q  large  and  q  small. 
However  this  statement  does  not  allow  establishing  an  optimum  ratio  for 
Q  and  or  an  optimum  value  for  either  one. 

The  peak  signal  strain  energy  in  the  sensor  springs  is  from 
page  381  of  the  AFCRL  proposal. 

€  -  Z2  /ti  Q 

peak  mech  4  {  w 


where 


peak  mechanical  energy  stored  in  the  sensor  springs 

inertia  efficiency  ratio  =  0.  861 

the  undamped  natural  frequency  of  the  sensor 

equivalent  gravity  gradient  input  =  10"^  sec"^  =  1  EU 

polar  moment  of  inertia  of  one  sensor  arm 
3.  561  X  10  ^  kg  m^ 

mechanical  Q  of  the  sensor  =  300 

sensor  resonant  frequency  =  220  rad/sec 


P 

T  = 


OJ  = 

o 

3r  = 


c 

zz 


Q  = 


CJ 

o 
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£  mech  =  1,  23  X  10”^^  joules  (v/att  sec)/EU  (5) 

P 

If  1.  0  %  of  this  mechanical  energy  is  absorbed  by  the  piezo¬ 
electric  transducer,  part  of  it  is  converted  to  electrical  energy  and 
this  is  determined  by  the  bending  coupling  factor  k^.  The  relation 
for  the  transducer  is 

peak  mech 

*peak  elect  100 


where: 


and  k^p  one  of  the  piezoelectric  constants  for  the  material  being 
considered  is 


k31  » 


Therefore 


k^  »  0.  065 


and 


peak  elect 


=  8  X  10 


•25 


joules 


The  piezoelectric  transducer  acts  as  a  capacitor  and  a 
reasonable  value  for  the  capacity  is 


C  =  2  nF 


(6) 
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This  provides  a  transducer  output  impedance  of 


~  C  ~  ohms 

o 

which  is  satisfactory,  see  Section  7. 

The  peak  transducer  output  voltage  will  be 


*  peak  ~  2.83  X  10  ^  =  28.3  nanovolts/EU 

This  signal  level  is  acceptable  in  a  practical  sensor  design. 

In  general,  a  higher  percentage  of  the  sensor  energy  would  be  stored 
in  the  transducer  and  a  lower  impedance  could  result.  Thus  the  trans¬ 
ducer  output  can  be  expected  to  be  in  the  low  hundreds  of  nanovolts 
per  EU  and  have  few  megohms  output  impedance. 

From  the  point  of  view  of  the  signal  level  the  preliminary  design 
parameters  are  satisfactory. 

Signal  Frequency  Considerations  —  The  sensor  signal 
frequency  must  be  above  about  10  Hz  to  allow  reasonable  component 
size  and  values  in  the  preamplifier.  Up  to  perhaps  100  Hz,  component 
sizes  decrease  inversely  with  increasing  frequency.  This  consideration 
indicates  that  the  frequency  should  be  as  high  as  possible.  The 
preliminary  specification  for  f^  =  35  Hz  is  a  convenient  working 
frequency  and  thus  satisfies  this  requirement. 


AFCRL  Proposal. 


(4)  Signal  Phase  Shift  Considerations  —  It  is  shown  on  page 
462  of  the  AFCRL  proposal  that  the  error  introduced  into  the  off 
diagonals  off  the  gradient  tensor  by  a  phase  error  (|)  is  given  by 

Zr. .  error  =  -SF  gin  2  4* 
ij 


where 


r=  GM/R^ 

4)  =  the  total  phase  error  of  the 
sensor  output  with  respect 
to  the  reference  coordinate 
system 


For  the  small  phase  angles  involved  this  can  be  assumed  to  be 
linear  and  thus  the  phase  accuracy  required  is  the  average  accuracy 
over  the  specified  10-second  integration  time. 

For  small  angles  and  for  a  sensor  operating  with  its  spin 
axis  horizontal  at  the  earth's  surface 


3  GM 

r. .  error  =  rr — 5-^9 
ij  2  1,3 


=  2250A(|)  EU 


If  it  is  desired  to  limit  F. .  error  to  1/3  EU  due  to  this  cause  then 

ij 


_5 

A4)=  7.41  X  10  radians  =  0.0042  degrees 

The  sensor  A<^  error  is  due  to  (1)  gradiometer  servo  including  the 
frequency  reference;  (2)  sensor  phase  error  due  to  a  shift  in  its 
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natural  frequency;  (3)  phase  pick-off  and  electronics.  If  these  are 
assumed  to  be  statistically  independent,then  0.  00242  degrees,  1  sigma, 
can  be  allowed  for  each  source.  This  is  the  error  that  will  be 
tentatively  allowed. 

The  phase  error  due  to  the  servo  and  frequency  reference  or 
due  to  a  change  in  sensor  resonant  frequency  is  shown  on  page  471  of 
the  AFCRL  proposal  to  be 

A4>  =  Q degrees  (3) 

o 

where  theAoj  may  be  either  the  shift  in  the  servo  and  reference 
frequency  or  the  shift  in  the  sensor  resonant  frequency.  Limiting 
the  phase  error  due  to  each  of  these  causes  to  0.  00242  degrees, 
the  frequency  control  accuracy  requirement  becomes  for  a  Q  of  300 

CjJ  1  l  A  ^ 

-  =  7.  1  X  10 

0) 

o 

The  phase  shift  error  is  a  linear  function  of  the  sensor  Q  and 

thus  a  low  Q  is  desirable  from  this  point  of  view.  However,  if  Q  is 

made  small  the  thermal  noise  equation  (1)  shows  that  must  be 

increased  to  compensate  for  any  reduction  in  Q  to  maintain  the  same 

thermal  noise  limit.  Increasing  C  requires  an  increase  in  the 

z  z 

physical  size  of  the  sensor  and  this  is  undesirable.  The  value  of  Q 
selected  (300)  provides  a  reasonable  physical  size  and  a  manageable 
phase  shift  problem.  Since  frequency  always  enters  the  phase  shift 
problem  as  the  ratioAoj/oj^  the  actual  frequency  selected  is  not 
important. 

The  basic  frequency  reference  stability  is  not  a  problem. 

General  radio  frequency  synthesizers  are  on  hand  that  are  good  to 
9 

one  part  in  10  .  A  preliminary  design  for  a  digital  servo  that  can 
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maintain  the  required  accuracy  has  been  made.  This  servo  will  be 
analyzed  in  detail  and  presented  in  the  final  report. 

A  preliminary  analysis  of  the  temperature  sensitivity  of 
the  sensor  resonant  frequency  has  also  been  made  and  the  required 
accuracy  appears  to  be  state  of  the  art.  This  analysis  will  be  com¬ 
pleted  and  presented  in  the  final  report. 

(5)  Sensor  Time  Constant  Considerations  —  The  thermal 
noise  equation  for  the  sensor  can  be  written -as 


where: 

Tg  =  sensor  time  constant 
Tj.  =  filter  time  constant 

(t^  +  T^)  =  T.  =  system  integration  time  specified  as 
10  seconds 

Consideration  of  this  equation  shows  that  the  smallest  value  of  (r 
results  if  IS  made  equal  to  10  seconds  and  the  filter  time  constant  is 
reduced  to  zero.  If  this  were  to  be  done  it  would  have  three  adverse 
effects.  These  are; 

•  The  sensor  Q  would  have  to  be  raised  which 
would  increase  the  signal  phase  shift  problem. 

•  The  would  have  to  be  lowered  which  would 
increase  the  electronic  circuit  problems  and 
increase  the  size  of  the  sensor  mounted  electronics. 

•  No  filtering  would  be  available  for  the  noise 
introduced  by  the  signal  preamplifier. 
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The  first  tv/o  iuc^rns  have  been  discussed  in  detail  in  previous 
paragraphs  and  the  assumed  values  for  Q  and  u  ^  have  been  deemed 
acceptable.  The  adequacy  of  the  filtering  after  the  preamplifier  will 
now  be  examined. 

Only  the  first  stages  of  a  low  level  —  low  frequency  amplifier 
need  be  considered  as  noise  sources  since  if  these  stages  have  a 
gain  of  10  or  more,  the  noise  contribution  of  the  following  stages  s 
negligible. 

As  shown  in  the  section  on  RGG  Sensor  Electronics,  the 
first  stages  cf  the  amplifier  to  be  used  consist  of  two  low-noise  field 
effect  transistors  and  a  low-noise  operational  amplifier.  The  gain  of 
this  combination  is  20  and  the  estimated  noise  will  be  17  nanovolts 
per  square  root  Hertz.  The  "square  root  Hertz"  refers  to  the 
equivalent  bandwidth  of  the  data  handling  system.  Since  the  preampli¬ 
fier,  amplifier,  frequency  converter  and  transmitter  are  all  relatively 
wide  bandwidth,  the  preamplifier  noise  bandwidth  will  be  determined 
by  the  equivalent  noise  bandwidth  of  the  output  filter. 

It  is  shown  on  pages  390  and  392  of  the  Hughes  AFCRL 
proposal  that  equivalent  noise  bandwidth  for  a  linear  filter  is 


noise  2 

The  previously  selected  parameters  provide  a  =  7.  27  seconds  and  the 
root-mean-square  amplifier  noise  will  be 

amplifier  =  4.  5  nanovolts  rms 

If  this  4.  5  nanovolt  amplifier  noise  signal  is  root-sum-squared  with  a 
signal  of  approximately  30  nv  for  one  EU  signal  the  4.  5  nv  is  seen  to 
be  negligible. 
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,  *  ’  Component  Size  and  Riridity  Consid.raH^.,  -  ,> 

frequency 

(P)  and  the  difference  mode  frequency  (o,^)  are  approximately: 

„2  _  Kq  K, ) 

O  C 


54  K, 


where: 

-  difference  mode  frequency  =  220  rad/sec 
P  =  sum  mode  frequency 

"  svs^im''®  t”'*'  °f  'he  difference  mode  pivot 

Itae  dTaig^”"  P'““  'ransducers  in  th?  base 

Kj  =  torsional  stiffness  of  each  pair  of  end  pivots  Ther, 
are  lour  "end  pivots"  in  the'^base  line  d'is.>'  " 

^  hg‘m2  of  •>'0  -nsor 

^  ?t53rx"Jo-3  “hgtT"  -nsor  stator 

The  sum  mode  frequency  should  be  lower  than  the  dif- 

rUiIitT^  f-ators 

imit  the  arm  inertia  mismatch  problem  as 'shown  in  another  section 
this  report.  In  order  to  demonstrate  size  and  rigidity  of  parts  it  is 

TanTd""  '  ',7  "’*'’0  and  that  the 

transducers  will  store  20  percent  of  the  peak  system  energy.  These 

th?e°o  ”7  ^  assumption 

the  equation  for  can  be  written 

,,2  3.08  K 

^o  =~-^  andp2,-- - E. 

‘-zz 
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where  is  the  torsional  stiffness  of  each  individual  pivot. 

The  torsional  stiffness  of  a  cylindrical  pivot  is 


rG  d 
32 


where: 

G  =  the  shear  modulus,  4.  5  x  10^^  for  beryllium  copper 
d  =  pivot  diameter 
t  =  pivot  length. 


Combining  the  last  two  equations 


A  32  w  ^  C 
a  _ _ o  zz 

f  ~  7,5  TT  G 


If  it  is  assumed  that  f  =  2  d,  this  becomes 


d 


3 


2.72  u  ^  C 

_ O  2Z 

G 


meters 


3 


and 


d 


=  2.15x10  ^ 


meters 


or 


d  =  0.0845  inches 
f  =  0.169  inches 
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It  is  reasonable  to  make  pivots  of  this  size  with  good 
control.  ® 

The  cross-section  area  of  each  pivot  will  be 


A  -  ^  (0.0845)^  -  5.6  x  10”^  in.^ 


Since  only  the  end  pivots  support  the  load, 
arms  of  1.563  kg  x  2.2  =  3.44  lb  each.  The 
pivot  will  be  at  Ig. 


then  4  pivots  support  two 
shearing  stress  in  each 


c _ 2  X  3.44 

4  X  5.6  X  10 


This  IS  very  safe  loading  in  any  expected  environment. 

Based  on  this  same  assumption  that  all  of  the  pivots  are 
Identical  and  20%  of  the  difference  mode  sensor  energy  is  stored  in  the 
piezoelectric  transducer  the  sum  mode  frequency  is 

/C  X  i.QS\ 


t 


Bum 


22.4  Hz. 


Thic  IS  a  satisfactory  sum  mode  frequency. 

These  computations  demonstrate  that  the  selected  w 

and  Q  allow  aafe  and  reaaonable  aenaor  pivot  designs  and  a  salisfao°ony 
sum  mode  frequency. 

ll^^sducer  Output  Impedance  Consideration.^  -  Tr, 
previous  paragraph  on  signal  level  considerations  it  was  shown  that 
satisfactory  signal  levels  were  obtained  if  only  1  percent  of  the  peak 
sensor  energy  was  stored  in  the  piezoelectric  transducer  and  the 
transducer  had  an  output  capacity  of  2  nF.  In  the  final  design  it  is 
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expected  that  approximately  20%  of  the  peak  system  energy  will  be 
stored  in  the  two  transducers,  10%  in  each.  With  the  higher  energy 
available,  the  signal  level  can  be  raised  above  the  previously 
calculated  or  the  output  capacity  can  be  increased  and  the  output 
impedance  lowered.  The  final  parameter  balance  will  be  determined 
by  the  design  of  the  piezoelectric  transducer. 

It  was  previously  shown  that  the  2  nanofarad  output  capacity 
gave  an  output  impedance  of  2.  27  x  10  ohms.  This  is  connected  to 
the  preamplifier  input  impedance  of  about  10^  ohms  shunted  by  a  few 
picofarads  input  capacity.  The  phase  shift  and  signal  loss  is  negligible 
for  this  situation  and  will  have  even  less  effect  in  the  final  design. 

The  transducer  output  impedance  is  satisfactory  using  the 
selected  sensor  parameters. 

(8)  Environmental  Considerations  of  coq  Selection  —  The 
linear  and  angular  vibration  spectrum  in  any  laboratory  or  vehicle 
will  contain  peaks  at  certain  frequencies.  These  peaks  are  usually 
due  to  nearby  rotary  equipment  operating  at  nearly  constant  speed.  In 
the  laboratory  the  peaks  are  la.gely  due  to  1725  RPM  (28.  8  Hz) 
induction  motors  driving  compressors  and  pumps.  In  an  airplane  they 
occur  at  engine  rotation  speed. 

The  selected  (35  Hz  or  2Z0  rad/sec)  is  well  removed  from 
the  expected  environmental  peaks  and  is  high  enough  that  effective 
vibration  isolation  is  relatively  easy. 

(9)  Spin  Bearing  Loss  Lionsiderations  —  The  sensor  spin 

bearing  drag  is  essentially  a  linear  function  of  The  spin  bearing 

specification  allows  the  mean  running  torque  to  be  as  high  as  5  x  10^ 
dyne  cm  at  17.  5  Hz  (110  rad/sec).  The  power  loss  in  the  bearing 
under  these  conditions  will  be 

W  =  Torque  X  Wg  =  0.  55  watts 
This  is  considered  to  be  an  acceptable  spin  bearing  loss. 
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(10)  Upper  Acceleration  Limit  —  The  upper  acceleration 
limit  that  the  arm  mass  balance  adjustment  devices  can  withstand 
without  danger  of  malfunctioning  is  approximately  50  g's.  Assuming 
that  the  balance  devices  can  be  located  within  1.  5  inches  of  the 
center  of  rotation  then  the  acceleration  is 


g's 


980 


50 


Thvs  the  selected  sensor  resonant  frequency  is  satisfactory  from  the 
balance  device  acceleration  limit  point  of  view. 
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2. 


RGG  Sensor  Electronics 


Electronic  circuits  have  been  designed  to  perform  all  of  the 
signal  and  logic  functions  required  at  the  sensor.  These  functions  are 
divided  into  the  following  groups: 

•  Piezoelectric  Transducer  Amplifier  and  Telemetry- 
Circuit 

•  Powe  r  Supply 

•  ®^lance  and  Gain  Control  Logic  Circuit 

Each  function  is  described  in  detail  in  the  following  paragraphs  and  the 
schematic  diagram  is  shown  in  Fig.  9. 

Piezoelectric  Transducer  Amplifier  and  Telemetry 

Circuit  - - - — ^ 

P^scription  —  The  voltage  out  of  the  piezoelectric  trans¬ 
ducer  undergoes  two  stages  of  amplification  by  bandpass  amplifiers 
Aj  and  A^,  and  drives  the  voltage-controlled  oscillators  (VCO).  The 
VCO  has  a  center  frequency  of  200  kHz  and  deviates  ±100  kHz  for  the 
largest  drive  signals.  The  VCO  output  is  buffered  by  and  drives 
the  output  telemetry  capacitor. 

Preamplifier  -  Amplifier  A^  and  field-effect  transistors  Qj 
and  form  a  high-input-impedance  operational  amplifier  with  an 
open-loop  gain  of  several  million.  Feedback  elements  R^,  C  .  R  . 
and  C^  establish  a  closed-loop  gain  of  20  at  35  Hz  with  6  dB  pU  octave 
rolloffs  above  and  below  this  frequency. 

Changes  in  DC  level  of  the  input  signal  at  P^  by  up  to  1.  2  volts 
are  tolerated  and  matched  by  a  corresponding  change  in  the  voltage  of 
P^.  Constant-current  diode  ensures  that  changes  in  DC  1>  vel  do 
not  affect  the  Ql  and  operating  current.  Transistor  provides 
46  dB  of  power  supply  isolation  for  the  preamplifier. 

P^ost  Amplifier  -  Field-effect  transistor  is  normally 
turned  off.  In  this  case,  amplifier  A^,  with  feedback  network  C  , 

Ry  R^.  has  a  gain  of  R^/30  K  at  35  Hz  and  a  6  dB  per  octave  roLff 
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EUectxDaic  Circ'+lt  Schetn4Hc 


below  35  Hz.  is  chosen  to  optimize  the  gain  of  the  system  so  that 
the  maximum  expected  input  signal  drives  the  VCO  at  its  maximum 
frequency  deviation,  ±100  kHz.  The  amplifier  also  may  be  operated 

in  a  low-gain  mode  by  turning  on  Q  . 

5 

^5'  ^5’  ^6'  ^6’  ^7  ^e-establish  an  accurate  DC  level 

(9.  2  volts)  for  the  VCO  input,  and  provide  additional  6  dB  per  octave 
rolloffs  above  and  below  35  Hz. 

Voltage- Controlled  Oscillators  -  VCO  with  timing  resistor  R, 
and  timing  capacitor  C.^,  produces  a  200  kHz  square  wave  when  the 
drive  voltage  at  pin  5  is  at  the  center  of  its  range,  9.  2  volts.  Deviations 
from  this  value  of  ±0.  15  volt  produce  frequency  deviations  of  ±20  kHz. 

Emitter-follower  buffers  the  VCO  output  providing  a  5  volt 
ieak-to-7)eak  200  ±100  kHz  square  wave  for  the  rotating  plate  of  the 
output  telemetry  capacitor. 

(2)  Specificaticns  —  Noise  is  generated  in  the  preamplifier 
primarily  by  transistors  and  Q^,  resistors  R^,  R^,  and  R^^,  and 

amplifier  A^.  The  equivalent  input  noise  voltage  per  root  Hz,  en,  is 
given  by 


2 

c  n 


4kT  Rj  +  4kT  R^^  + 


4kT 

+ 


+  2 


9 
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“  Reactance  of  piezoelectric  sensor 

'Q  =  Equivalent  gate  noise  Voltage  per  root 
Hz  of  Qj  and 

A1  ■  equivalent  input  noise  voltage  per  root  Hz 
of  amplifier  A^. 

^FQ~  t^ansconductance  of  Q^,  Q^. 

For  our  devices,  X_  =  2.  5xl0^« 

Lr 

6x10*^  volts/yHz 
8x10*^  volts/v/Hz 

1400n 

Rjj  =  7500n 

Rg  =  lo^n 

Rj  =  5000n 

Rjq=  looon 


and  therefore  -  17  nv/N/Hz.  This  corresponds  to  0.  1  EU  for  a 

sensor  with  45  nv  per  EU  using  a  7.  27  second  output  filter  integration 
time. 


Gam  and  Phase  -  Shift  Drift  -  Since  the  open  loop  voltage  gains 
of  A^  and  A^  are  greater  than  10^,  gain  drift  and  phase  drift  in  the 
preamp  and  post  amp  will  only  be  due  to  changes  in  value  of  the  feedback- 


loop  and  network  capacitors  and  resistors.  R  .  C  R  C  C  R 

1  1’  2’  2'  ^3'  '^3' 

^4>  ^5»  Cg,  C^,  R^,  and  R,^.  Since  these  feedback  networks  have 
Q's  near  unity,  the  percent  phase  drift  and  percent  gain  drift  are 


roughly  equal  to  the  sum  of  the  percent  drifts  in  values  of  these 
capacitors  and  resistors.  Metal  film  resistors  are  available  with 
drifts  in  value  less  than  25  parts  per  million  per  degree  Centigrade 
(ppm/°C).  Capacitors  (such  as  metalized  film)  are  available  with 
drifts  less  than  25  pprn/°C  and  with  leakage  resistance  greater  than 
5  X  10^°  n. 

In  an  environment  temperature-controlled  to  1°C,  the  phase 
and  gain  of  the  preamp  and  post  amp.  will  drift  less  than  0.  5  milliradian 
and  0.  05  percent  respectively. 

The  frequency  drift  of  the  VCO  is  about  100  ppm/°C.  The  FM 
gain  drift  will  be  less  than  one-fourth  this  value  since  our  peak  excur¬ 
sions  are  ons-fourth  the  total  oscillator  range.  This  will  add  another 
0.  0025  percent  gain  drift  to  the  system. 


Linearity  —  Because  of  the  high  open-loop  gains  of  amplifiers 
and  A_,  and  the  high  quality  of  available  feedback  elements^  non- 
linearity  will  arise  primarily  in  the  voltage -controlled  oscillator,  VCO. 
For  excursions  of  ±20  kHz  about  200  kHz,  the  voltage-to-frequency 
non-linearity  is  about  0.  2%.  The  FM  detector  outside  the  gradiometer 
will  be  a  phase-locked  loop  using  an  identical  VCO  as  reference 
oscillator.  This  should  reduce  the  overall  voltage-to-frequency-to- 
voltage  non-linearity  to  less  than  0.  05%. 


Power  SuT 


Power  for  the  sensor  electronics  is  brought  in  continuously 

by  a  rotating  power-input  capacitor.  The  outside  fixed  plate  is  driven 

with  80  volts  rms  at  500  kHz.  The  inner  rotating  plate  charges  storage 

capacitor  Cg  through  rectifiers  and  D^. 

The  total  power  supply  current  is  16  milliampere.  To  supply 

this  current,  the  power-input  capacitors  must  be  100  to  300  pf.  This 

will  charge  capacitor  Cg  to  24  volts  DC.  Current  diode  Dg  and  the 

24 -volt  zener  clamp  on  Cg  create  a  16  ma  supply  to  Q,^  that  is  nearly 

independent  of  variations  in  the  capacitors  of  the  power-input  capacitor. 

_3 

A  and  Q_  form  a  shunt  regulator  with  an  impedance  of  less  than  10  D. 

^  '  5 

Since  the  parallel  resistance  of  is  greater  than  10  D,  the  supply 
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voltage  will  vary  less  than  10  %  for  a  one  %  variation  in  capacitance 

of  the  power-input  capacitor.  Use  of  the  shunt  regulator  rather  than  a 
series  regulator  keeps  changing  circuit  current  requirements  from 
affecting  the  reaction  force  applied  to  the  power-input  capacitor. 

c.  Balance  and  Gain- Control  Logic  Circuit 

(1)  Description  —  Sixteen  binary  bits  of  information  are  hold 
by  the  sensor  electronics  package.  These  bits  control  the  operation 

of  the  balance  tubes  and  adjust  the  sensor  amplifier  gain.  The  user 
may  supply  a  new  set  of  sixteen  bits  at  any  time  by  telemetering  them 
in  serially  on  the  500  kHz  power  supply  input. 

Each  successive  bit  is  entered  by  briefly  interrupting  the 
500  kHz  power.  A  100  |jLsec  interruption  is  a  logical  one  and  300  |jisec 
interruption  is  a  logical  zero.  Successive  binary  bits  are  shifted  down 
a  sixteen-bit  shift  register.  After  sixteen  new  bits  are  supplied,  these 
are  available  on  the  sixteen  shift-register  parallel  outputs. 

The  first  twelve  bits  select  one  of  twelve  balance  tubes  and  bits 
thirteen  and  fourteen  select  the  polarity  of  the  current  (or  no  current 
at  all)  supplied  to  the  selected  tube.  Bit  fifteen,  if  set,  lowers  the 
gain  of  the  sensor  amplifier  by  turning  on  Qg. 

(2)  Logic  Elements  —  The  elements  G  and  S  are  of  a  new 
Motorola  family  of  low-power  MOS  logic.  They  consume  only  ten 
nanowatts  per  package  and  operate  well  on  10  volts.  The  logic  levels 
are  roughly  10  volts  for  logic  one  and  roughly  ground  for  logic  zero. 

The  outputs  will  sink  or  source  several  milliamperes. 

(3)  Control  Logic  —  When  500  kHz  power  is  supplied, 
point  Pg  (at  the  input  to  oscillates  between  ground  and  10.  5  volts. 
The  output  of  word  gate  Gj  also  oscillates,  and  through  D6  charges 
P4  to  10  volts. 

When  the  500  kHz  is  interrupted,  point  P^  stays  at  10.  5  volts, 
the  output  of  Gj  stays  at  ground,  and  P^  goes  to  ground.  Thus,  a 
negative  pulse  occurs  at  P^  when  the  500  kHz  is  interrupted. 
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Inverter  and  ccnverts  this  to  a  positive  pulse  at  Pg  and 
again  a  negative  pulse  at  P^. 

Nand  gate  G^  outputs  a  logical  one  (±10  volts)  except  fo)  the 
brief  interval  (30  usee)  following  a  positive  transition  of  P^.  The 
30  psec-pulse  initiates  the  200  psec  one-shot  multivibrator  formed 
by  Gg  and  G^.  The  trailing  edge  of  the  one-shot  output  at  Pg  clocks 
the  edge-clocked  shift  registers  Sj,  S^,  S^,  and  S^.  These  shift 
registers  are  connected  in  series  as  one  16  bit  serial-in,  parallel  out 
shift  register.  The  input  is  at  line  PD. 

The  shift  register  is  clocked  200  fisec  after  each  interruption 
of  the  500  kHz  supply.  If  the  interruption  is  greater  than  200  fisec, 

Pg  is  at  logical  zero  (ground)  when  it  is  clocked  into  the  shift  register 
(through  line  PD'.  If  the  interruption  was  less  than  200  ^isec,  Pg  has 
returned  to  one  when  it  is  clocked  into  the  shift  register. 

(4)  Balance  Driver  —  A  balance  tube  is  driven  by  setting  its 
bit  (1  through  12)  to  one.  The  direction  of  current  flow  is  selected  by 
setting  either  bit  13  or  bit  14  to  one.  This  sets  point  P^  for  that  tube 
at  10  volts,  and  sets  point  P^q  or  P^  ^  to  ground.  One  milliampere 
then  flows  through  the  tube  in  the  direction  selected  by  bits  13  and  14. 

If  bite  13  and  14  are  both  zero,  then  P^^  and  P^  ^  are  both  at  10,  5  volts 
and  no  current  flows. 

Since  the  outputs  of  G.^  and  Gg  will  not  saturate  when  sinking 
two  milliampere s,  if  more  than  one  balance  tube  is  to  be  driven  at  one 
time,  the  current  amplifiers  in  bo>.  B  must  be  included. 

Transistor  Q.^  and  word  gates  G.^  and  Gg  interrupt  balance-tube 
current  for  12  millisec  after  each  bit  is  entered.  This  suppresses  all 
the  current  while  the  sixteen  new  bits  are  entered.  These  are  entered 
at  a  2  kHz  rate  (500  psec  between  bits). 

3,  Sensor  Arm  Mass- Unbalance  and  Inertia-Unbalance  Mechanical 

Adjustment 

It  is  the  purpose  of  this  section  to  estimate  the  limits  to  which 
sensor  arm  mass  and  inertia  unbalance  can  be  adjusted  and  the  stability 
limits  of  that  adjustment. 
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a.  Mechanical  Mass  Balance  Adiustment 

Experience  has  shown  that  a  0.  25  gram,  80  pitch,  screw 
can  be  manually  adjusted  to  within  *1/10  turu.  This  provides  an  incre¬ 
mental  mass  balance  adjustment  of 

_4 

(Amh)  =  8x10  gram  cm.  . 


If  the  nonrotating  sensor  is  vibrated  a.c  Zuy  (=uj  )  the  output  signal  in 

S  G 

terms  of  equivalent  gravity  gradient  is 


^Amh) 


=  cos 

nc 


(w^t  +  O') 


where: 


A 

g 


Amh 

qC 


peak  acceleration  in  "g" 
gravity 

magnitude  of  the  sensor  mass  unbalance  vector 

sensor  arm  inertia  efficiency  x  arm  polar 

4  2 

inertia  =  3.  066  x  10  gm  cm 
220  rad/sec 
phase  angle 


If  the  nonrota*-ing  sensor  has  a  residual  mass  unbalance  of 

-4  -3 

8x10  gm  cm  and  it  is  vibrated  at  5  x  10  g  the  sensor  output  will 

be 


^(Amh)  =  129  cos  (w^t  +  a  )  (EU) 

This  is  an  easily  readable  signal  and  thus  the  manual  adjustment  of 
mass  unbalance  can  be  made  as  proposed. 
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b. 


Balance  Tube  Mass  Unbalance  Adjustment 


Each  of  the  mercury  filled  balance  tubes  has  a  total 
range  of  1.4  x  10“^  gm  cm.  If  the  gap  is  initially  centered,  the  adjust¬ 
ment  range  will  be  ±7  x  10“'^  gm  cm.  It  is  proposed  to  use  8  balance 
tubes  in  each  direction  on  each  arm  giving  a  balance  range  with  the 
mercury  tubes  of 

Amh  =  ±56  X  10'^  gm  cm/arm/axis,  available 

The  mechanical  mass  unbalance  adjustment  is  well  within  the 
range  of  the  balance  tube  system  and  assembly  of  the  sensor  can  be 
completed  and  final  balance  accomplished  on  the  rotating  sensor.  The 
balance  tube  gap  position  can  be  adjusted  to  a  small  fraction  of  a  milli¬ 
meter  on  a  rotating  sensor  by  the  electronic  circuits  discussed  in 
another  part  of  this  report. 

It  will  be  shown  in  a  subsequent  paragraph  of  this  section  that 
the  arm  mass  unbalance  allowed  is 

(Amh)  =  4  X  lO""^  gm  cm  maximum  allowed. 


This  resolution  is  easily  obtainable  by  the  balance  tubes. 

The  mass  unbalance  adjustment  signal  for  a  rotating  sensor 
can  be  calculated  from  the  equations  on  page  337  of  the  AFCRL  proposal 
after  the  equations  are  corrected  as  shown  below. 


r 


accel 


Ag(Amh) 

Zr]C 


cos  (2w  t  +  p) 
s 


when  vibrated  at 


Since  it  is  desired  to  place  the  mass  center  of  the  arm  at  the  center 
of  rotation  of  the  sensor  the  concept  of  simply  matching  the  two  arm 
unbalances  and  thus  reducing  (Amh)  to  zero  can  not  be  used. 
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If  we  assume  that  we  would  like  to  balance  to  10%  of  the  maximum 
desired  unbalance,  and  further  assume  that  we  can  vibrate  the  table  at 

_3 

5  X  10  g,  then 

A  =  5  X  10"^g 
(Amh)  =  4x10^ 
r]C  -  3.  066  X  lO"^ 

r  ,  =  3.2  cos  (2a)  t  +  P  )  EU 
accel  s 

This  is  a  good  working  level  signal. 

c.  Inertia- Unbalance  Meehan. cal  Adjustment 

It  is  shown  in  paragraph  7a  of  this  section  entitled  "Gradient 
Errors  Due  to  Sum- Mode  Mismatch"  that  the  sum-mode  mismatch 
coefficient  is  defined  by  the  equation 


where: 


K,  A, 
1  A  1 


DjS  +  Kj 


^2®  "  "^2 


P  ^1 


kp 

^1 

°1 

K, 


sum-mode  mismatch  coefficient 
end  pivot  stiffness  for  arm  1 
end  pivot  damping  for  arm  1 
end  pivot  stiffness  for  arm  2 
end  pivot  damping  for  arm  2 
polar  inertia  of  arm  1 
polar  inertia  of  arm  2 
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If  lihe  damping  is  assumed  to  be  zero  and  Kj  equal  to  then 


where; 

AI  =  the  difference  in  the  polar  moment  of  inertia  of  the 
two  arms 

I  =  the  nominal  polar  inertia  of  either  arm 


Al  =  4  r  (m  Ar  ) 
s  s  s  s 

If  the  same  0.  25  gm  -  80  pitch  balance  screws  are  assumed  to  be 
placed  at  a  radius  of  2  cm  and  rotated  1/10  turn,  then  for  the  arms 
with  a  polar  moment  of  inertia  of  35,  610  gm  cm^ 

k  =  -  4  X  2  X  0.  25 

P  I  "  80  X  10  X  35,  610 

=  0.  7  X  10'^ 


a 
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It  is  shown,  in  the  previously  referenced  section  of  this  report  that 
a  =  10“^  is  desired  In  this  sensor  to  limit  the  2u  torque  from  the 
bearing  from  coupling  into  the  sensor  and  to  limit  signal  overloading 
due  to  torques  at  the  spin  frequency,  We  can  achieve  an  inertia 

unbalance  adjustment  to  1/10  that  required  without  the  use  of  the 


balance  tubes. 

It  is  also  shown  in  the  previously  referenced  section  that  if  a 
torque  is  applied  to  the  sensor  rotor  at  the  frequency  a  gradient 

error  signal  appears  at  the  sensor  output  due  to  the  arm  inertia 
mismatch.  The  peak  magnitude  of  this  error  is 


or 


since  T./I 
a  c 


is  the  angular  acceleration  of  the  rotor  case  a. 
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where: 


=  sensor  output  error  in  EU 
kp  =  s\un-mode  mismatch  factor  AI/I,  lO"^ 

=  sensor  difference  mode  frequency 
=  sensor  sum-mode  frequency  ' 

T)  =  sensor  arm  inertia  efficiency,  0.  861 

I,  =  polar  moment  of  inertia  of  one  sensor  arm, 

^  2 
35610  gm  cm 

I  =  rotor  case  polar  moment  of  inertia  not  including 

^  2 
the  arms,  132700  gm  cm 

a  =  angular  acceleration  of  the  rotor  case 

If  it  is  assumed  that  the  sensor,  nonrotating,  is  oscillated  on  a  test 
fixture  so  that  accelerometers  1.  5  feet  from  the  center  have  an  output 
of  5  X  10  g  at  a  frequency  of  then  the  angular  acceleration  at  the 
center  of  the  test  fixture  will  be 

o  =  =  0.  107  rad/sec^ 

Notice  that  the  same  drivers  and  accelerometers  are  assumed  that 

were  proposed  for  the  mass  unbalance  tests. 

If  it  is  further  assumed  that  oj  /6  =2,  which  is  near  the 

o  *^0 

design  value,  the  output  signal  will  be 

r  =  27.  6  EU 
c 


to  k 


P 


This  is  a  substantial  signal  and  the  sensor  can  easily  be  adjusted 
10  ^  as  previously  suggested. 
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d.  Masa  —  Unbalance  Limits  Based  on  Linear  Vibration 
Estimates 

Theoretically,  the  only  vibration  frequencies  that  introduce  2u 

s 

torques  due  to  arm  mass  unbalance  of  the  RGG  are  those  at  and 

3  u  . 
s 

The  RGG  errors  due  to  differential  arm-mass-unbalance  are 
derived  in  detail  in  another  section  of  this  report.  Using  this  analysis 
and  the  best  current  estimates  of  the  vibration  power  spectral  density 
(PSD)  stable  platform  we  find 

PSD  at  17.5  Hz  =  4.  5  x  ] g  ^/Hz 
PSD  at  52.  5  Hz  =  8.  9  x  10"  g  ^/Hz 

Based  on  our  contractually  specified  integration  time  of  10  seconds 
the  equivalent  noise  bandwidth  is  0.  05  Hz.  Substituting  these  into  the 
reference  equation  33,  the  standard  deviation  of  the  acceleration  is 

0-  ,  =  1.  16  X  lO’^g 

accel  * 

This  vibration  level  estimate  is  quite  low  and  we  believe  the 
RGG  sensor  should  operate  in  a  more  severe  environment.  We  estimate 
that  the  vibrational  power  spectral  density  could  be  as  much  as  10  times 
as  high  as  the  values  shown  and  thus  the  acceleration  could  be 


“^accel  possible 


3.  68  X  10”^g 


The  reference  shows  in  equations  34  and  35  that  the  gradient 
tensor  error  terms  due  to  arm  mass  are,  with  specific  but  reasonable 
assumptions: 

g(Amh) 

=  "3^0 

g(Amh) 

_4 

Assuming  a  mass  unbalance  term  of  (Amh)  =  4x10  gm  cm;  qC  =  30, 
2 

660  gm  cm  ;  and  the  hypothesized  worst  case  vibration  environment 
the  tensor  error  verms  are: 

0-..  =  0.  235  EU 
ij 


0-..  =  0.  272  EU 
11 

_4 

Thus  a  (Amh)  mass  unbalance  of  4  x  10  gm  cm  in  the  worst 
case  hypothesized  environment  with  the  RGG  mounted  on  the  specified 
vibration  isolation  system  produces  an  acceptable  uncertainty  level. 

e.  Arm  Mass  Balance  Stability  Expected 

The  arm  mass  balance  can  be  expected  to  change  slightly 
due  to  temperature  shifts,  temperature  cycling,  temperature  gradients, 
and  creep.  It  has  been  shown  that  the  mass  unbalance  can  be  adjusted 
to  values  well  below  the  design  limit  of  4  x  lO”"^  gm  cm.  The  following 
paragraphs  discuss  the  problems  of  maintaining  the  unbalance  to  this 
limit. 

If  the  halves  of  the  arms  grow  unsymmetrically  and  differ  by 
lO"^  inches  the  unbalance  would  be 

-4 

(Amh)  =  20  X  10  gm  cm 
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Thus,  the  balance  lubes  can  accomnnndate  only  1.4  microlnches 
dlHerential  arm  growth.  The  importance  of  material  slabihnat.on 
before  assembly  is  obvious. 

Little  is  known  about  material  creep  in  the  fraction  of  micro¬ 
inch  region,  but  we  would  not  expect  such  stability  over  periods  of 
more  than  a  few  days.  Also  with  each  cool-down  and  warm-up  we  can 
expect  shifts  of  this  magnitude.  Thus,  we  can  expect  to  do  a  mass 

unbalance  trim  on  each  warm-up. 

It  should  be  recognized  that  it  is  quite  unlikely  that  we  can  mass 

balance  the  sensor  by  vibrating  when  it  is  nonrotating  and  have  this 
same  mass  unbalance  when  rotating  to  better  than  about  ±4  x  10  gm  cm 
One  guide  to  the  type  of  mass  stability  that  we  can  expect  can  be 
derived  from  the  characteristics  of  the  Autonetics  0-9  gyro.  This  gyro 
has  a  rather  large  free  rotor  on  a  spherical  air  bearing.  The  rotor  .s 
beryllium  with  a  high  density  ring.  In  addition,  one  copper  cup  and 
one  manganin  cup  are  attached  to  the  rotor.  This  rotor  is  at  least 

comparable  to  the  RGO  rotor.  Autonetics  quotes  the  axial  mass 

unbalance  repeatability  of  this  gyro  as  0.  01°/hr/g.  Other  data  are 


I  =  1650  gm  cm 

6  2 . 

H  =  2.  5  X  10  gm  cm  /sec 
k  s  .2  cm  (estimated) 
mass  =  I/j^?  =  1650/4  «  400  gms 

Unbalance  torque  =  H  x  drift  rate 

=  0.  125  dcm 
-4 

=  1. 25  X  10  gm  cm 
in  a  one  g  field 
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The  mass  unbalance  ration  will  be  approximately  a  linear  function  of 
mass  for  a  device  of  this  type  so  we  can  expect  about  4  times  the 
unbalance  they  expect.  That  is,  the  mass  unbalance  stability  we  can 
expect  in  our  sensor  comparing  it  to  the  Autonetics  G-9  gyro  is 

_4 

5x10  gm  cm 

Thus  we  seem  to  have  reasonable  expectations  for  the  RGG. 

Another  comparison  was  made  with  a  Bell  Ae  rosy  stems  gyro. 
This  is  also  a  relatively  large  gyro.  In  addition  to  the  inertia  wheel 
the  Bell  gyro  has  a  spin  motor  stator  on  the  sensitive  element.  Bell 
quotes  0.  05  /hr/g  shifts  from  warm-up  to  warm-up  but  states  that 
the  gyro  will  not  exceed  0.  003  /hr/g  from  its  initial  value  in  any  one 
warm-up.  Thus  Bell  has  high  initial  shifts  but  expects  excellent 
stability  for  each  warm-up. 

Another  calculation  that  is  important  in  the  mass  unbalance 
stability  estimation  is  that  due  to  the  effects  of  temperature  gradients. 
The  effective  length  of  the  arms  is  approximately  5  inches.  The  total 
cross  section  is  approximately  0.  5  x  2.  125  in.  sq.  The  thermal 
expansion  coefficient  for  aluminum  is  approximately  20  x  10”^cm/cm/°C 
and  its  thermal  conductivity  is  0.  37  cal/cm°C  sec.  The  thermal 
resistance  of  the  arm  is 


th 


=  ai  = 

A 


1.  2  C/watt 


If  there  is  a  uniform  temperature  gradient  along  the  arm  the  mass 
unbalance  due  to  a  1°C  difference  will  be  that  due  to  the  average 
temperature  of  one  half  the  arm  being  1/2°C  hotter  than  the  other 
half.  The  hot  half  will  expand  and  shift  1/2  the  mass  of  the  arm  to 
produce  an  unbalance  of 

=  I  f 


-4 


800  X  10  gm  cm 
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If  we  wish  to  limit  this  mass  shift  to  2  x  10“*  gram  cm  then  the 
temperature  difference  between  the  two  ends  must  not  exceed 

AT  =  0.  0025°C 

Since  the  thermal  resistance  is  1.2°C/watt  this  limits  the  power  Tow 
through  the  arm  to 


P  =  0.  002  W 

We  can  expect  the  power  flow  into  the  balance  tubes  {~0.  008  W)  to 
cause  a  small  thermal  mass  unbalance  transient. 

f.  Conclusion 

It  has  been  shown  that  the  arm  mass  and  inertia  unbalance 
can  be  adjusted  to  acceptable  values  for  the  worst  case  expected  environ¬ 
ment.  It  is  shown  that  similar  devices  in  similar  environments  main¬ 
tain  mass  unbalance  shifts  to  the  limits  required  by  the  RGG.  Thus 

the  RGG  is  expected  to  meet  its  accuracy  specification  in  the  expected 
environment. 

Vibrational  Motion  Environment  Estimates 


Estimates  of  both  the  translational  and  angular  vibration  environ 
ment  are  necessary  for  determining  design  parameters  of  both  the  RGG 

sensor  and  the  VIALS  (Vibration  Isolation,  Alignment  and  Leveling 
System). 

As  stated  in  our  proposal,  we  have  selected  the  airborne  appli¬ 
cation  as  being  the  most  critical.  The  choice  has  been  somewhat 
arbitrary,  however  it  is  felt  that  the  vehicle  motion  environment  for 
the  airborne  application  is  generally  likely  to  be  as  severe  as  for  a 
surface  ship  or  a  submarine.  Differences  between  each  vehicle's 
environment  may  exist  as  regards  the  specific  details  such  as  the  fre¬ 
quency  distribution  of  vibrational  power. 
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The  RGG  sensor  is  sensitive  to  specific  narrow-band  frequencies 
of  translational  vibration  which  are  related  to  its  design  rotational  speed. 
The  sf^sor  speed  can  be  varied  sornewhat  to  avoid  operation  at  or  near 
large  resonant  peaks  in  a  particular  vehicle's  vibration  spectrum.  The 
sensor  is  also  sensitive  to  any  and  all  frequencies  of  angular  vibration 
normal  to  its  spin  axis  and  is  independent  of  its  design  rotational  speed. 
However,  the  angular  vibration  transmitted  to  the  sensor  via  its  angular 
isolation  system  is  somewhat  unrelated  to  the  vehicle's  environment. 

The  purpose  of  the  VIALS  is  to  provide  attitude  stabilization  and  i'> 
isolate  angular  rate  inputs  to  the  sensor.  Its  own  design  characteristics 
are  as  dominant  in  determining  the  angular  vibration  inputs  as  is  the 
vehicle's  exact  angular  rate  power  spectrum. 

We  have  attempted  to  obtain  motion  environment  data  for  sub¬ 
marines  and  ships  and  as  yet  have  not  been  able  to  obtain  such  data. 

We  have,  therefore,  proceeded  to  generate  design  requirements  for 
the  sensor  and  VIALS  based  on  vibration  data  for  the  KC- 135  aircraft. 

It  is  felt  that  this  approach,  while  partially  one  of  expediency,  is 
probably  conservative  and  will  result  in  a  sensor  and  VIALS  design 
which  will  certainly  typify  designs  for  other  vehicle  applications.  Our 
goal  has  been  and  will  continue  to  be  one  of  providing  a  design  suitable 
to  any  application.  Realizing  that  any  design  must  be  tailored  to  some 
extent  to  its  particular  application,  we  will  have  the  flexibility  to  vary 
the  sensor  spin  speed  to  accommodate  different  vibration  spectrums. 

In  addition,  the  VIALS  system's  ability  to  provide  appropriately  low 
levels  of  angular  vibration  is  to  8o:ne  extent  independent  of  the  environ¬ 
ment.  As  additional  vibrational  motion  data  for  other  vehicles  becomes 
available,  it  will  be  compared  to  our  KC-135  data  and  the  system's 
ability  to  accommodate  such  other  environments  will  be  determined. 

a.  KC-135  Vibration  Data 

Several  references  have  been  used  to  glean  the  necessary 
translational  and  angular  vibration  spectral  data.  They  are: 

(1)  "Summary  of  Power  Spectral  Results  —  KC-135 

AF55-3118  Flight  Fatigue  Program,  "Boeing  Co. 

Report  D6-7102,  June  1962. 
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(2) 


(3) 

(4) 


Environmental  Vibration  Measurements  on 
Components  of  the  AN/ASQ-92  Subsystem  on  the 
Kt.-135  Airplane"  Boeing  Co.  Report  T6-3604 
August  1966. 


''KC-135  Power  Spectral  Vertical  Gust  Load 
Analysis,  "  Volume  I,  Boeing  Co.  Report 
D6- 18252  July,  1966  (AFFDL  TR-66-57). 


'^Measurements  of  the  Motions  of  a  Large 
Sw'jpt-Wing  Airplane  in  Rough  Air"  by  R.  H.  Rhvne 
NACA  TN4310  September  1958.  ’ 


The  data  of  reference  (2)  above  was  used  to  obtain  both  trans¬ 
lational  and  angular  vibration  data  in  the  range  of  frequencies  above 
20  Hz.  This  report  documents  results  of  a  flight  test  evaluation  of 
a  camera  mount  subsystem  for  the  reconnaissance  airplane.  The  flight 
condition  was  a  level  flight  at  M  =  0.  80  and  35,  000  ft  altitude.  The 
test  was  conducted  "during  a  photographic  flight  condition  considered 
typical  by  flight  crews"  according  to  the  report. 

The  data  of  reference  (1)  above  was  used  to  obtain  the  transla¬ 
tional  vibration  powe.r  spectrum  in  the  0  to  10  Hz  range.  Since  the 


purpose  of  these  tests  was  to  measure  flight  fatigue  properties  of  the 
aircraft,  severe  flight  conditions  were  sought  for  the  tests.  The  r.  m.  s. 
gust  velocities  for  the  tests  from  wh'ch  our  data  is  taken  was  5.  38  ft/sec. 
Hence  the  power  spectral  densities  were  scaled  down,  per  the  tech¬ 
niques  described  in  the  foregoing  reference  (4),  to  simulate  a  2.  0  ft/sec 
r.  m.  s.  gust  velocity.  The  2  ft/sec  r.  m.  s.  gust  velocity  is  classified 
as  a  "mild  turbulence"  condition  which  should  typify  flight  conditions 
for  airborne  gravimetry.  This  scaled  down  translational  vibration 
data  was  also  used  to  generate  angular  vibration  data  in  the  0  to  10  Hz 
frequency  range.  This  was  done  by  making  use  of  the  aircraft's  mass 
and  aerodynamic  properties  of  reference  (3)  above  and  techniques  of 
reference  (4)  above  to  transform  aircraft  c.  g.  normal  acceleration 
into  aircraft  pitch  rate  data. 


Generally,  the  high  frequency  data  (<  20  Hz)  showed  similar 
angular  and  translational  data  for  all  three  (x,  y,  z)  axes.  Only  air 
craft  normal  (vertical)  translational  acceleration  data  was  available 


from  reference  (1)  above.  It  is  normal  to  expect  that  aircraft  lateral 
and  longitudinal  translational  accelerations  in  the  lower  frequency  range 
will  be  no  more  severe  than  normal  acceleration.  Similarly,  low 
frequency  yaw  and  roll  motion  i  \mld  be  expected  to  be  no  more  severe 
than  pitch  motion  for  an  autopilot-controlled  large  aircraft  such  as 

'.heKC-135.  Hence,  for  purposes  of  gradiometer  and  VIALS  design, 

we  have  assumed  the  same  vibrational  power  spectrums  for  all  three 
angular  and  translational  axes.  This  assumption  is  also  justified 
because  the  gradiometer  must  perform  in  any  spin-axis  orientation. 

The  translational  and  angular  vibration  power  spectrums  from 
references  (1)  and  (2)  above  were  plotted  and  the  spectrums  for 
frequencies  between  10  and  20  Hz,  for  which  no  data  was  available  from 
the  references,  were  approximated  by  fitting  smooth  curves  to  the 
overall  data.  The  resultant  data  are  shown  as  the  solid  lines  of 
Figs.  10  and  11. 

Inertial  Platform  Base  Motion  Spectra 

Having  established  estimates  for  the  vehicle  motion 
environment,  estimates  were  also  needed  for  the  vibration  levels  at  the 
base  of  the  inertial  platform.  As  per  our  September  1971  proposal,  a 
passive  vibration  isolation  mount  is  envisioned  to  serve  primarily  the 
purpose  of  attenuating  the  translational  vehicle  vibration  and  secondarily 
to  attempt  to  isolate  angular  vibrations  into  the  inertial  platform.  Pre¬ 
liminary  estimates  of  vibration  attenuation  performance  of  the  base 
motion  isolation  system  were  made.  The  resulting  translational 
and  angular  vibration  power  spectra  occurring  at  the  mounting  base  of 
the  inertial  platform  are  shown  as  dashed  lines  on  Fig.  10  and  Fig.  11. 
Note  that  at  the  sensor  critical  frequencies,  its  spin  frequency 
(17.  5  Hz)  (Wg  and  are  the  sensor  critical  frequencies  as  regards 
Its  error  caused  by  arm  mass  unbalance)  the  vibration  power  level  is 
attenuated  by  approximately  2.  5  orders  of  magnitude.  Due  to  mass 
unbalance  coupling  (i.  e.  center  of  m='^s  of  base  motion  isolation 
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Fig.  11.  Angular  Rate  Power  Spectrum 


system's  payload  not  coincident  with  its  elastic  spring  center)  the 
angular  rate  power  a.t  the  inertial  platform  base  is  slightly  above  the 
aircraft  power  levels  up  to  8  Hz,  however  at  the  higher  frequencies, 
angular  rate  power  is  attenuated.  This  result  is  not  considered  to 
create  a  serious  problem  since  the  region  of  amplification  is  within 
the  servo  bandwidth  of  the  inertial  platform. 

c.  Inertial  Platform  Stable  Element  Motion  Spectra 

Estimates  of  the  vibration  environment  of  the  platform 
stable  element  are  also  required.  For  a  first  preliminary  estimate, 
the  translational  v'bration  spectra  for  the  stable  element  has  been 
assumed  to  be  the  same  ar  for  the  base.  Obviously,  this  is  not  pre¬ 
cisely  true  since  the  platform  gimbal  structure  will  have  vibration 
resonances  which  will  alter  the  base  spectra  somewhat.  Any  such 
resonances  will  have  to  be  controlled  by  specifying  appropriate  vibra¬ 
tion  transmission  characteristics  of  the  inertial  platform. 

A  preliminary  estimate  of  the  angular  vibration  spectra  of 
the  stable  element  has  been  made.  Several  sources  of  angular 
rate  excitation  were  investigated.  These  included  angular  disturbance 
torques  caused  by: 

(1)  Gimbal  mass  unbalance 

(21  RGG  Sensor  case  mass  unbalance 

(3)  Friction  torque  disturbance. 

Our  estimates  for  both  of  items  1  and  2  above  resulted  in 
negligible  angular  rate  contribution.  However,  the  friction  torque 
disturbance  did  prove  to  be  a  problem. 

First,  there  is  no  known  simple  analytical  method  for  determin¬ 
ing  the  angular  transfer  function  between  base  and  stable  element  angu¬ 
lar  rate  power  spectra.  This  is  due  to  the  stick- slip  nonlinearity 
characteristic  associated  with  the  conventional  ball  bearing  supported 
platform  gimbals.  We  have,  for  a  preliminary  estimate,  used  a 
simplified  method  which  isolates  this  non-linearity  outside  the  plat¬ 
form's  servo  stabilization  loop,  which  is  not  exact.  Other  methods  of 
analyzing  this  problem  exist  but  are  very  time  consuming  and  expensive 


108 


since  they  all  would  involve  setting-up  special  corr^uter  simulations. 

The  results  of  our  simplified  analysis  are  shown  as  the  dotted  line  in 
Fig.  11.  The  accuracy  of  this  result  is  not  known  but  it  is  prob¬ 
ably  somewhat  unconservative.  It  is  interesting  to  note  that  if  the 
RGG  sensor  were  hard  mounted  to  the  platform  stable  element,  the 
resulting  rotational  field  error,  based  on  our  simplified  estimate, 
would  be  only  2.  5  EU. 

d.  Additional  Angular  Isolation  Needed 

Based  on  the  above  simplified  analysis,  it  is  obvious  that 
additional  angular  isolation  is  needed  to  meet  the  overall  specified 
accuracy  of  1  EU.  There  are  two  basic  approaches  to  providing  such 
additional  angular  isolation.  One  would  involve  providing  an  additional 
stage  of  isolation  between  the  platform  stabilized  element  and  the  sensor. 
Several  configurations  of  systems  which  could  perform  such  a  function 
are  described  elsewhere  in  this  report.  The  second  approach  v/ould  be 
to  improve  upon  the  isolation  performance  of  the  gimbaled  inertial 
platform.  The  above  simplified  estimates  were  based  on  ball  type 
gimbal  bearings,  whose  payload  capability  was  assumed  to  be  200  lbs 
and  whose  friction  torque  level  was  10  inch-ounces.  If  the  friction 
torque  were  reduced  by  utilizing  hydrostatic  gas  gimbal  bearings,  it 
appears  that  a  suitable  reduction  in  friction  could  be  achieved  v/hich 
would  provide  adequate  performance. 

5.  Preliminary  Gradiometer  Designs 

Prior  to  and  during  the  preparation  of  our  September  1971 
proposal,  much  consideration  was  given  to  configuring  a  gradiometer 
which  would  provide  for  isolation  of  angular  vibrations  from  the  rotating 
sensor  arm  pair,  thus  providing  a  reduction  of  the  inherent  rotational 
field  errors  of  the  gradiometer.  The  need  for  such  isolation  had  been 
anticipated  at  that  time.  Our  analyses  of  the  vehicle  motion 
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environment  in  Section  III-C-4  demonstrates  the  desirability  of  providing 
such  isolation.  After  considering  many  alternative  schemes  for  supporting 
the  rotating  arm  pair  in  such  a  way  as  to  provide  adequate  angular 
vibration  isolation,  we  recommended  the  neutrally  buoyant  rotating 
sphere  in  our  September  1971  proposal.  However,  practical  details  of 
the  design  had  not  been  studied  in  depth. 

Realizing  that  many  of  the  other  alternate  schemes  also  appeared 
practical  and  potentially  feasible,  we  first  decided  to  study  some  of  the 
important  design  details  of  the  neutrally  buoyant  rotating  sphere  to 
ascertain  if  it  was  truly  a  practical  and  feasible  mechanization 
approach.  The  study  of  such  details,  e.  g.  ,  ease  of  manufacturing  and 
assembly,  centering  of  the  sphere,  signal  and  power  routing,  etc. 
resulted  in  a  preliminary  design  of  the  configuration.  These  design 
details,  including  sketches  of  the  system  are  summarized  at  the  end 
of  this  section  as  configuration  "A”  (Fig.  12).  Briefly,  it  became 
apparent  that  this  configuration,  although  feasible  to  build,  would  have 
some  very  difficult  assembly  and  fluid  filling  problems.  Also,  although 
our  brief  computer  simulation  results  indicated  it  would  provide  what 
appeared  to  be  adequate  angular  isolation,  we  were  reticent  to  fully 
accept  these  results  since  there  is  general  lack  of  understanding  of  the 
behavior  of  rotating  fluids.  For  the  above  reasons,  we  then  considered 
other  alternatives  and  prepared  similar  preliminary  designs  for  two 
additional  sensor/angular  isolation  configurations.  These  two  config¬ 
urations,  "B"  and  "C"  (Figs.  13  and  14),  have  a  common  concept  which 
differs  from  configuration  "A".  In  configuration  "A",  the  sensor 
spherical  "float"  rotates  and  thereby  the  angular  isolation  occui-s  in 
the  spinning  frame.  As  such,  the  long  time  average  angular  orienta¬ 
tion  of  the  sensor  arms  is  maintained  in  precise  alignment  with  the 
sensor  spin  bearing  via  the  self  aligning  characteristics  of  the  .system. 

That  is,  the  sensor  float  is  designed  to  have  its  axis  of  maximum 
moment  of  inertia  aligned  along  the  torsional  axis  of  the  sensor  arms. 
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As  such,  the  gyroscopic  and  viscous  torques  of  this  system  cause  the 
sensor  arms  to  rotate  about  this  preferred  axis. 

In  configurations  "B"  and  "C",  the  sensor  arms  and  case  rotate 
in  the  spin  bearing  and  the  angular  isolation  is  provided  by  supporting 
the  spin  bearing  stator  on  a  two  degree- of- freedom  mounting  system. 
Thus,  the  angular  isolation  occurs  in  the  non- spinning  frame.  Gyro¬ 
scopic  torques  still  tend  to  provide  angular  stabilization  but  the  system 
has  nutation  instabilities  and  is  subject  to  mass  unbalance  and  other 
disturbance  error  torques.  Thus,  some  form  of  angular  restraint  is 
necessary  to  provide  the  required  long  term  angular  stabilization. 

In  configuration  "B”,  the  two-axis  suspension  is  provided  by  a 
ring  gimbal  very  similar  to  that  of  a  two  degree- of- freedom  gyro. 

Both  flexural  pivots  and  journal-type  hydrostatic  gas  bearings  were 
considered  for  the  gimbal  bearings..  It  was  found,  however,  that 
flexural  pivots  were  not  appropriate.  This  was  because  the  torsional 
spring  rate  had  to  be  so  low,  to  accomplish  the  required  angular  isola¬ 
tion,  that  pivots  could  not  be  designed  to  support  the  load.  Sketches 
depicting  configuration  "B"  are  shown  at  the  end  of  this  section. 

In  configuration  "C",  the  two-axis  suspension  is  provided  by  a 
set  of  four  spherical- segment  hydrostatic  gas  bearing  thrust  pads. 

The  pads  are  located  at  the  surface  of  a  spherical  housing  encasing  the 
sensor  and  spin  bearing  and  positioned  at  the  corners  of  a  circum¬ 
scribed  equilateral  tetrahedron.  This  suspension  provides  an  iso¬ 
elastic  support  for  the  sensor.  The  suspension  is  restrained  to  have 
only  two  rotational  degrees  of  freedom  about  axes  normal  to  the  sensor 
spin-axis  by  a  system  of  restraint  wires  as  depicted  in  the  sketches  of 
configuration  "C"  at  the  end  of  this  section. 

A  comparision  chart,  shown  in  Table  III,  summarizes  the  various 
tradeoffs  considered  in  comparing  these  preliminary  designs.  For 
purposes  of  program  continuity  and  organization,  we  have  continually 
carried  along  a  particular  sensor  configuration  as  the  baseline  to 
provide  a  common  ground  for  analysis  and  design  of  the  various  sensor 
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subsystems  and  serve  as  a  basis  for  comparison  of  alternative 
configurations.  At  the  inception  of  the  program  the  baseline  was  con- 
figuration  "A”.  As  we  studied  other  configurations,  it  became  apparent 
that  configuration  "C"  was  preferred  and  it  was  tentatively  adopted  as 
the  baseline.  At  the  time  oi  changing  from  "A"  to  "C,  it  was  thought 
that  "C"  could  provide  the  necessary  angular  isolation  while  maintaining 
the  required  apin-axis  alignment  by  using  only  passive  spring  cvnd  viscous 
restraints.  We  have  since  performed  further  analyses  on  configuration 
"C"  and  found  that  an  active  servo  control  system  will  be  required. 

These  analyses  are  not  yet  completed  or  documented.  The  specific 
mechanization  which  the  servo  control  system  must  take  is  still  under 
study.  Thus,  the  selection  of  a  baseline  configuration  must  await  com¬ 
plete  knowledge  as  to  the  various  tradeoff  considerations. 

Some  thought  has  been  given  to  the  design  of  the  sensor  arms, 
support  pivots,  and  transducer  mounting.  Detail  studies  and  analyses 
are  now  in  progress,  but  are  not  yet  sufficiently  complete  to  state  a 
firm  recommended  design.  A  design  similar  to  the  one  illustrated  in 
the  configuration  "C”  sketch  is  now  under  consideration. 

Neutrally  Bouyant  Sphere  Rotating  Gravity  Gradiometer 

Notes  for  Preliminary  Sketch  (General):  Almost  entirely 
6061  aluminum.  Stress  relieved  before  final  machining.  Most  screws 
and  balance  weights  brass  which  has  almost  the  same  coefficient  of 
expansion  as  aluminum. 

(1)  Alternate  screws  and  taper  pins.  Approximately  24 
around  outer  case.  Approximately  6  around  air 
bearing. 

(2)  O-rings.  Used  as  temporary  seal  only. 

(3)  Brass  balance  screws.  Their  coefficient  of 
thermal  expansion  almost  identical  to  aluminum 
and  can  get  adequately  low  magnetic  susceptibility. 

Will  not  gall  with  aluminum. 


114 


(4)  Spanner  nuts  each  end.  Note  that  lower  one 
provides  a  gimbal  tilt  stop  that  just  matches 
the  pivot  to  float  clearance  at  top.  Select 
top  nut  for  balance.  Proper  seals  not  yet 
provided  at  these  points. 

(5)  Outer  case  of  float.  6061  aluminurrj. 

(6)  Flex  leads  —4  required. 

(Pivot  =  Gnd);  1  +  Battery;  1  -  Battery;  1  Signal;  ]  Logic. 

(7)  Bellows.  Can  only  install  at  one  end.  Bellows 
shown  is  not  adequate. 

(8)  Spanner  nut. 

(9)  Insulator  if  needed  or  desired. 

(10)  Beryllium  copper  pivot  with  aluminum  or  brass 
extender.  Extender  keeps  float  control  cen¬ 
tered  as  temperature  varies.  This  may  not  be 
necessary.  Pivot  itself  can  be  shortened.  This 
will  help  reduce  suspension  point  shift  as  temper¬ 
ature  varies. 

(11)  Air  bearing  rotors. 

(12)  Flotation  fluid.  Carbon  tetrachloride  in  present 
design  but  there  are  a  number  of  other  possibilities. 

(13)  May  have  to  change  slightly  so  that  we  havf'  a  half 
sphere  as  a  cup  for  half  flotation  duriz-g  a  s^embly. 

Have  not  provided  a  good  seal  yet  at  this  point. 

(14)  Central  plate  and  central  post.  6061  aluminum 
of  course.  Make  these  first  and  assemble  and 
balance  arms  and  add  internal  guts  to  these. 

Will  want  to  make  a  handling  fixture.  The  central 
plate  can  have  large  holes,  bosses,  posts,  etc,  as 
required.  Internal  structure  not  yet  designed. 

(15)  Fill  port  with  non-locking  taper  plug.  Held  in 
place  by  spanner  nut  with  hole  to  accomodate 
extractor  post  of  plug. 

(16)  Non-locuing  taper,  centers  and  levels  pivot  post. 
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(17)  Batteries.  12  —  0.  225  ampere  hour,  1.  3  volts /cell. 
Nickel -cadmium.  Must  be  provided  with  case 
fixed  magnetic  shield. 

(18)  Battery  carrier.  Visualize  batteries  epoxy  set 
in  an  aluminum  ring. 

(19)  Power  input  capacitor  plates. 

(20)  Stand-off  insulators  for  19. 

(21)  Exhaust  air  deflector  attached  to  main  frame. 

(22)  Air  bearing  stators  attached  to  main  frame. 

(23)  Backup  motor  stator  iron.  Overlaps  statur  stack 
to  act  as  magnetic  shield. 

(24)  Motor  stator.  TranCore  T,  0.  007  in.  thick. 

Should  be  able  to  buy  standard  punching. 

(25)  Mafi[petic  leakage  shield. 

(26)  Motor  stator  mount  attached  to  main  case. 

(27)  Epoxy  seals 

(28)  Float  angle  pickoff. 

(29)  Pivot  nut  with  Allen  Wrench  socket. 

Neutrally  Bouyant  Sphere  Rotating  Gravity  Gradiometer 
Preliminary  design  data: 

(1)  Float 


Material 

6061  aluminum 

Outside  Diameter 

(6.  50  in.  ) 

16.  51  X  10  ^  m 

Inside  Diameter 

(6.  00  in.  ) 

15.  24  X  10"^  m 

Displaced  Vol 

(2356  cc) 

2.  356  X  10  ^  m^ 

Average  density 

(1.  6  gm/ cc) 

1.6x  10  ^kg/m 

Mass  (3770  gms) 

3.7  kg 
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Polar  Inertia 

(1.  02  X  10® 
gm  cm^) 

1.  02  X  10"^ 
kg/m2 

Diameter  Inertia 

(8.  0  X  10^ 

gm  cm^) 

8.  00  X  10"® 
kg/ m2 

Arms 

Mass  (Mallory 
1000  &  Al) 

(700  gms) 

0.7  kg 

Inertia 

(20,  000  gm  cm^) 

2.  0  X  10"® 
kg/ m2 

Inertia  Efficiency 

0.  7r 

Peak  Torque, 

1  EU  Field 

7.  0  X  10"^®  Nm 

Balance  Screv's 

Brass 

Central  Pivot 

Material 

Beryllium  Copper 

Length,  Active 

(0.  100  in.  ) 

Diameter 

(0.010  in.  ) 

K  ,  Torsional  St 
Stiffness 

(7.  2  X  lO”  dcm/ 
rad) 

7.  2  X  10  ®  Nm/rad 

K^,  Bending 
Stiffness 

(9.  4x10^  dcm/ 
rad 

9.  4  X  10”®  Nm/rad 

Stress  when 

Lilting  Float 

(106,  000  psi) 

Torsional 

Freedom 

(2.  39  deg) 

±0. 0417  rad 

Bending  Freedom 

(2.39  deg) 

±0.  0417  rad 

Max  Tension 

Stress  Bend 

(34,  400  psi) 

Max  Shear 

Stress  Torsion 

(  ~12,  000  psi ) 

Hang-Off  at  20°/hr  (0.  6l  deg) 

Input  Rate 

0.  010  rad 
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Fluid 


Carbon  Tetrachloride  (tentative  selection) 


Density,  20°C 

(1.6  gms/cc) 

1.6x10  ^  kg/m^ 

Viscosity,  |i 

(0.  009  poise) 

9. Ox  10  ^  Nsec/m 

Gernaanium  Tetrachloride 

Density 

(1. 84  gms/cc) 

1. 84  X  lO"^  kg/m^ 

Viscosity, 

Damping  Gap,  h 

Physical 

(0.  025  in.  ) 

6.  35  X  10"^  m 

Effective, 

10  Hz 

2.  7  X  10  ^  m 

Effective, 

20  Hz 

1.  9  X  10  ^  m 

Effective, 

40  Hz 

1.35xl0^m 

Damping  Coefficient 

Low  Frequency 

(5,  950  dcm  sec) 

5.  95  X  lO”^ 

Nm  sec 

10  Hz 

(8.  340  dcm  sec) 

«.  34  X  10“'* 

Nm  sec 

20  Hz 

(11,  800 
dcm  sec) 

11.  80  X  10""* 

Nm  sec 

40  Hz 

(16,  680 
dcm  sec) 

1 6.  68  X  1 0“^ 

Nm  sec 

Volume  of 

Fluid 

(55  cc) 

5.  5  X  10"^  m^ 

Coefficient  of 
Expansion 

(1.  2  X  10"^  °C) 

» 

Net  Volume 
Change 

(1.8  cc) 

!.  8  X  10"^  m^ 

AT  =  (124  -  70) 
=  54°F 
=  30OC 


Natural  Frequenciea  and  Damping  Ratios 


Torsional 

Undamped 

(0.  135  Hz) 

0. 35  rad/sec 

Bending 

Undamped 

(0.  172  Hz) 

1 . 08  rad/ sec 

Damping  Ratio 
Torsion 

0.  024  (low 
frequency) 

Damping  Ratio 
Bending 

4 

0.  034  (low 
frequency) 

Spin  Motor  Estimates 

Excitation 

Frequency 

180  Hz 

Number  of  Poles 

12 

Synchronous  Speed 

30  RPS 

Stator  and  Rotor 
Material 

Tran-Core  T 
(.  007  ) 

Mass  of  Iron 

(1.  83  lbs) 

0.  828  kg 

Mass  of  Copper 

(0.  60  lbs) 

0.  414  kg 

Power  During 

Kun  Up 

20  watts 

Power  During 
Operation 

10  wattj 

Locked  Rotor 
Torque 

(3.  0  in.  -oz) 

2.  12  X  10"^  Nm 

Running  Torque 

(0.  5  in.  -oz) 

3.  54  X  10"“^  Nm 

c. 

(1) 


Basic  Data  For  Baseline  Configuration 
Mass  and  Inertia  Characteristics  of  Arms.  Rotor, 
and  Stator 


Wt 

gms 

I 

XX 

2 

gm-cm 

I  ■ 

'7  2 

giri-cm 

I 

7,Z 

2 

gm*cm 

Arm  No.  1 

Arm  No.  2 

Rotor  (all  except  arms) 

1563 

1563 

2750 

4,  990 

35,  610 

103,  400 

35,610 

4,990 

103,400 

35,  600 
3  5,  600 
132,700 

Total  Rotor 

Stator 

5876 

3766 

144, 000 

289,  200 

144,000 
289; 200 

203, 900 
252,600 

Total  Rotor  plus  Stator 

9642 

433,  200 

433,200 

456, 500 

Rotor  plus  Stator  Mass  Unbalan^c^  »  1 

gram-cm 

Spin  Frequency  =17.5  Hz 

=110  rad/ sec 

On- Axis  Sprint!  Restraint  of  Air  Pad  ' 

'Gimbal"  -  Such 

that  20°/hr  transport 
iecrce. 

ana  eartn  a  - 

(K  =  2  X  10^  dync-cm/rad  =  Hn/0) 

Daiviping  Coefficient  of  Air  Pads  =  500  dyne-cm-sec 


Tnnrtial  Platform  Friction  =  10  in-ozs 


Tnertial  Platform  Mass  and  Inertia  Properties 


Payload  (incl.  stable  element) 

Middle  Gimbal 
Outer  Gimbal 

Core  (50  lbs  plus  top  half  of 
isolation  mount  50  lbs) 

T 

f  _  10  in-ozs 

Mass  unbalance  aim  4 q q  qzs 


Radius  of 
Gyration 


ZOOJ 

8" 

40# 

10" 

60# 

12" 

100# 

14" 

001 56  in 

Base  Motion  Isolation  System 

Spring  arm  of  movint  =  17.  5 
Damping  ratio  for  translation  =  0.4 
Damping  ratio  for  rotation  =  1.0 
Mass  unbalance  arm  =  0.025  ft  =  0.  3  in 


Sensor  Parameters 

Difference  mode  resonant  frequency  =  35  Hz 
Difference  mode  Q  =  300 
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Isoelasticity  of  Spherical  Gimbal  Support  With  Hydrostatic 
Bearing  Pads 


One  possible  configuration  of  the  prototype  moving-base  gravity 
gradiometer  employs  a  set  of  four  pressurized  hydrostatic  bearing  pads, 
arranged  on  the  surface  of  a  sphere,  to  support  the  gradiometer  stator 
with  two  degrees  of  angular  freedom  for  angular  rate  isolation  from  the 
base  to  which  the  gradiometer  case  is  mounted.  It  is  desired  that  the 
hydrostatic  support  system  be  isoelastic  under  linear  acceleration. 

It  is  the  purpose  of  this  analysis  to  show  that  a  particular 
spherical  configuration  of  four  hydrostatic  bearing  pads  will  provide 
isoelastic  support.  The  four  pads  will  be  found  to  be  located  at  the 
comers  of  an  equilateral  tetrahedron  inscribed  within  the  sphere. 

Consider  a  pair  of  initially  concentric  spheres  of  radius  R  and 

R  -  h  ,  with  centers  at  the  origin  of  coordinates,  where  h  «  R. 

o 

Now  let  the  inner  sphere  undergo  a  displacement: 

6  =  16  ■^16  +k6 

X  y  z 


Then  the  change  in  clearance,  Ahp,  between  the  spheres,  at 
any  point  P(R,  6,  <())  on  the  outer  sphere,  will  be  given  by  the  negative 
projecuon  of  6  on  the  radius  vector  to  the  point  P. 


(V  6  + 

'px  X 


Y  6 

py  y 


+  Y 


pz 


Where  the  are  the  direction  cosines  of  the  point,  P,  with  respect 
to  the  coordinate  axes. 

If  we  now  assume  that  a  hydrostatic  bearing  pad  centered  at  any 
point  P  produces  a  radial  (i.  e. ,  normal  to  the  sphere  surface)  restoring 
force,  F  ,  directly  proportional  to  the  change  in  clearance  at  P,  then: 

F  =  -  T  K  Ah 
P  r  p 


IZQ 


Resolving  this  radi&l  force  along  the  coordinate  axes: 


F  =  -  KAh  {i  Y  +  i  Y  +  k  V  ) 

P  p  'px  ^py  ’pz' 


Substituting  for  from  above  and  writing  the  result  in  matrix  form: 


/  2 

/  Y  Y  Y 

/  px  *px  ’py 

Y„^  Y  A 
px  pz  \ 

2 

f 

Y  Y  Y 

py  px  Vy 

Y  Y 

py  ^pz 

i 

5 

y 

\  Y  Y  Y  Y 

\  pz  ‘px  ipz  ’py 

We  now  have  an  expression  for  the  force  due  to  a  single  pad  in  terms 
of  its  direction  cosines.  This  may  be  written  as: 


Y  •  Y  •  6- 

pi  pj  j 


Then,  for  a.  system  of  N  pads,  at  arbitrary  locations: 


Where  the  p  and  j  summatipns  have  been  interchanged.  We  can  now 
define  a  set  of  matrix  elements: 


N 

~  k  ''^pi  ''^pj 

P 


Such  that: 


N 

F.  =  F  .  =  Kk  Q..  6. 
1  Lu  pi  ij  j 

p  j 

or: 


for  the  resul'iant  foice  of  N  pads.  The  are  the  direction  cosines  of 
the  force  vector  with  respect  to  the  displacement.  The  isoelasticity 
condition  requires  that  the  force  vector  be  in  the  direction  of  the  dis¬ 
placement.  That  is,  for  isoelasticity, 

F  =  Kk6 

which  requires  that: 
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or 


Y  •  Y  • 
pi  PJ 


6..  = 

ij 


1  for  i  =  j 
0  for  i  j 


where  6^^.  is  the  kronecker  delta. 

fac.  ,  co„fig„aa.io„  of  four  pads  on  the  aur 

point  Vt"' 

points,  p,  are  as  shown  below: 


Pad 

X 

a 

^py 

^pz 

1 

sin  ij; 

cos  ijj 

0 

2 

sin  ijj 

-cos  ijj 

0 

3 

-sin  ijj 

0 

cos  (|f 

-sin  i|; 

0 

-cos  i|j 

As  shown  in  the  sketch  below,  the  fou 
inscribed  tetrahedron  but  the  value  of 


r  points  are  at  the  corners  of  an 
ijj  is  as  yet  unspecified. 


Forming  the  indicated  products  and 
expression; 

N 


a. . 
ij 


summing  over  p  in  the 

k  ^pi  "'^pj 
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The  a-matrix  becomes: 


^IZ 


"Z2 


32 


4  ,  2 
j^sin 


0 


0 


The  i,  diagonal  aa  a  noauit  oi  a  iortuiloua  choice  o,  coordinatea 

The  re,u.rementa  of  the  iaoelaatlcity  condition  thnn  ■ 

identically  if,  and  only  if:  ’ 


1 


or: 


F rom  which: 


and: 


4  sin  ijj  =  2  cos^  =  k 

tan^  ijj  =  i 

tan  ijj  =  gjjj  45OJ 

sin  ijj  =  ^  (=  ^an  30°) 
cos  ijj 

k 

ijj  ^  35°  16' 
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The  resultant  restoring  force  becomes; 


?  =  |kt 

the  vairlfVf  '"r  k'  (see  preceding  sketch,  that 

he  value  of  ,  found  to  be  reguired  for  isoelasticity  is  identically  egual 

to  that  required  for  the  tetrahedron  to  be  equilateral. 

under  th"  "  -‘l-ogonal  set  of  six  pads  is  also  isoelastic 

onal  t  '  “  (s  i-'eresting  to  note  that,  if  the  orthog- 

on  e  h  ‘'d'“'  “  suoh  that  three  of  the  six  pads  are 

on  each  side  of  the  plane,  and  such  that  each  of  the  six  radius  vectors 

orms  an  angle,  e.  with  the  symmetry  plane,  then  the  value  of  e  wi 
e^identical  to  t  e  value  of  ,  found  above  for  the  tetrahedral  set  ,1.  e. , 

.  16  ).  For  this  sui-pad  set,  the  force  is  given  by; 


2K  6 


The  above  results  can  be  extended  to  show  that  any  set  of  2q 

pa  s  (w  ere  q  ,s  any  integer  greater  tha;.  or  equal  to  2).  with  proper 

symmetry  about^a  symmetry  plane,  and  with  each  radius  vector  foxing 

n  ang  e  4,  (-35  16')  with  the  symmetry  pla,-ie,  will  be  isoelaslic.  It 

can  further  be  shown  that  the  resultant  force  for  any  such  system  wi 
be  giv>n  by;  ' 


F  =  fK6 


In  the  limit,  as  q  the  system  would  become  a  pair  of  rings 

of  radius  R  cos  4,  and  separated  by  a  distance  2R  sin  4,. 
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7. 


Spin  Bearing  Requirements 


Three  of  the  following  subsections  analyze  the  gravity  gradient 
errors  that  could  be  introduced  by  the  sensor  spin  bearing.  The  fourth 
subsection  is  the  spin  bearing  specification  that  is  based  on  the  analysis, 
practical  adjustment  limits  and  the  expected  environment. 
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Gradient  Errors  Due  to  Sum-Mnde  Mismatch 

In  the  rotating  gravity  g radiometer  the  arms  may  rotate 
in  either  of  two  fundamental  modes:  (1)  the  sum  mode,  in  which  the 
arms  move  together  relative  to  the  sensor  case;  and  (2)  the  difference 
mode,  in  which  the  arms  rotate  relative  to  each  other.  It  is  the  differ¬ 
ence  mode  which  is  excited  by  the  external  gravity  gradient  and  which 
genercites  an  output  signal. 

In  an  ideal  sensor,  spin  axis  disturbance  torques  excite  only  the 
sum  mode  and  produce  no  spurious  output  signal.  However,  to  the 
extent  that  imperfection  exists  in  the  matching  of  the  sum- mode  inertia- 
spring  constants  of  the  arms  (sum-mode  mismatch),  spin  axis  distur¬ 
bance  torques  will  excite  the  difference  mode,  resulting  in  a  gradient 
measurement  error. 

The  equations  of  motion  are  written  as: 

'l®  +  I  =  (B-A)r.. 

A  1  C  IJ 

I2S  ^2+  +  ®i)+ (D2S+ K  )(82-  e  )  =  (A-B)r. . 

c  c  IJ 

IcS^0^  +  (DiS+Ki)(0^-  0j)  +  (D2S  +  K2)(e^-02)  = 

for  arm  1,  arm  2  and  the  sensor  case,  c,  respectively.  D  ,  K  refer 
to  the  difference  mode  pivot  and  D^.  and  D^.  refer  to  the^two 
sum  mode  pivots.  A.  B  are  the  arm  transverse  inertias,  T..  the 

gradient  input  and  the  spin-axis  disturbance  torque  acting^on  the 
sensor  case. 

Defining: 


^o  ^  V  ^  ^o'  ^1  -  +  Kj;  A^  ^  D^S  +  K 
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and  writing  \';he  equations  in  matrix  form: 


These  equations  are  solved  for 


=  sum  mode  frequency 
=  difference  mode  frequency 

Qp  =  sum  mode  Q 
Q  =  difference  mode  Q 
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The  normalized  equivalent  gradient  output  in  the  presence  of 
disturbance  torque,  T^(S),  is  ;hen 


and  .he  equivalent  gradient  output  due  to  apin  axis  dieturbanee  torque 
and  sum  mode  mismatch  alone  is 


Defining 


,  A 

Kp  -  sum  mode  mismatch  coefficient 


This  expression  gives  the  frequency  resuanse  nf  fho 

^  ^  response  of  the  sensor  output  sienal 

to  aptn  axte  diehtrbanoe  torquea  in  terms  of  an  equivalent  gradient 
amphtude.  Only  that  portion  of  the  output  within  the  narrow  pass  band 
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at  the  difference  mode  frequency,  is  interpreted  as  gradient  error. 
To  find  the  gradient  error  amplitude  due  to  disturbance  torque  at  the 
difference  mode  frequency  we  evaluate  jr^l  at  S  = 


where  Qp  »  1  has  been  assumed.  This  value  of  differs  from  previ¬ 
ously  published  values  by  the  factor  (I^/B-A).  Representative  parame¬ 
ter  values  for  the  baseline  sensor  are: 


2 

Ij  =  35,  600  gm-cm 

2 

A  =  4,  990  gm-cm 

2 

B  =  35, 610  gm-cm 

2 

=  132,  700  gm-cm 

=  220  rad/sec 

P  =  110  rad/sec 

o 


for  whicL 


r  = 

f 


E.U. 


where  T^  is  in  dyne-cm. 

Then,  in  order  to  limit  sum-mode  mismatch  errors  to  less  than 
0.1  EU,  spin  axis  disturbance  torque  uncertainties  should  be  held  to 
less  than  50  dyne-cm,  at  the  difference  mode  frequency. 
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Since  the  sensor  response  characteristic  also  has  a  high-Q  peak 
at  the  sum  mode  frequency,  it  is  necessary  to  consider  limits  oh  spin- 
axis  disturbance  torque  at  the  sum  mode  frequency.  The  response  at 
the  sum  mode  frequency  is: 


Then  the  two  response  amplitudes  are  in  the  ratio 


For  representative  values: 


Q  =  300 
=  1000 

^  =  13.3 


for  a  given  disturbance  torque  amplitude. 
Or 


r 

P 


E.U. 
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If  we  wish  to  limit  the  output  at  the  sum  mode  frequency  to  the 
equivalent  of  25  EU  due  to  sum-mode  mismatch,  then  we  can  allow 
up  to  1000  dy;ie-cm  disturbance  torques  at  the  sum  mode  frequency. 


b.  Torque  Variation  With  Eccentricity  in  a  Fluid  Journal 
Bearings 

In  any  spin  bearing  now  being  considered  for  the  rotating 
gravity  gradiometer,  a  thin  fluid  tilm  (liquid  or  gas)  will  exist  between  the 
bearing  surfaces  and  this-  fluid  film  will  produce  a  viscous  drag  torque 
on  the  rotor.  Furthermore,  this  drag  torque  will,  in  general,  vary 
with  the  load  applied  to  the  bearing.  Clearly,  spin  axis  torque  varia¬ 
tions  can  result  in  gradiometer  measurement  errors.  It  is  the  purpose 
of  this  analysis  to  define  the  viscous  torque  variation  with  bearing 
eccentricity  ratio.  For  convenience  a  cylindrical  journal  bearing  is 
assumed  of  sufficient  length,  L,  that  end  effects  may  be  neglected. 

For  any  Newtonian  fluid  the  shear  stress  on  unit  surface  is 
proportional  to  the  velocity  gradient  normal  to  the  surface 

dF  ^  ^ 

dA  ^  dz 


where  the  constant  of  proportionality  is  the  viucosity,  |x.  The  viscous 
torque,  dT,  on  an  element  of  surface,  dA,  of  a  journal  of  radius,  R., 
is 


dT  =  hiR  ^  I  dA 
R. 


HIT 

where  |  is  the  velocity  gradient  at  r  =  R.  When  tho  lubricant  film 
R 

thickness,  h,  is  very  much  less  than  R,  the  velocity  gradient  is 
approximately 


dU  ^  U  _  Rn 
dr  ~  h  h 


(h«R) 


where  is  the  journal  spin  velocity. 
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Then 


dT  =  de 

h 


When  the  journal  is  displaced  from  its  concentric  position, 


h  =  h  (1  +  «  cos  6) 
o 


where  «  is  the  eccentricity  ratio, 


dT  = 

o 


de 


l  +  t  cos  6 


and  the  total  torque  on  the  journal  is 


T  = 


3 

nR  LS2  I  _ 

h  I  1  + 

O  “'O 


de 


e  COS  e 


“  h 


and,  if  we  write  the  torque  at  zero  eccentricity  as 


o  h 

o 


then 


T  = 
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Fo^  small  eccentricities 


T  =  T 


approx,  for  e  -►.0 


which  shows  that,  for  a  periodic  variation  in  eccentricity,  ®  ^  ®  q  sin  wt, 
the  change  is  torque 


AT  =  T  -  T  =  (1  -  cos  2  ut) 

o  4  ' 


consists  of  a  constant  term  plus  a  double-frequency  term. 


Returning  to 


T  = 


we  have 


d€ 


T  6 
o 


{1-€2) 


3/2 


and,  letting  a  =  acceleration. 


da  d€  da 


T  6 
o 


{1-€2) 


_ 

3/2  da" 


where  de  /da  is  in  the  form  of  a  compliance.  For  the  present  discussion 
we  may  assume  a  constant  compliance 


da 


-r-  =  k  =  constant 


«  =  ka 


144 


Then 


§lE 

da 


T  k^a 
o  _ 


(1-kV) 


371 


and  we  see  that  the  sensitivity  of  bearing  torque  to  changes  in 
acceleration  is  a  function  of  the  total  acceleration.  Current  repre 
sentative  values  of  parameters  in  the  above  equation  are: 

4 

T  =10  dyne -cm 
k  =  0.2  per  g 

a  =  I  2 

which  gives: 

J  rn 

—  =  425  dyne-cm/g 


as  the  sen?itivity  to  changes  in  acceleration  when  operating  in  a 
nominal  I  g  field.  Gradiometer  measurement  errors,  resulting  from 
spin  axis  torque  variations,  are  discussed  elsewhere  in  this  report. 
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c»  Phase  Errors  Due  to  Spin  Bearing  Disturbance  Torques 

The  speed  control  servo  of  the  rotating  gravity  gradiome- 
ter  really  serves  two  basic  functions  -  (1)  It  maintains  the  sensor  spin 
velocity,  at  precisely  one-half  the  tuned  resonant  frequency,  co 

of  the  sensor  difference  mode  (velocity  servo)  and  (2)  it  maintains  a 
precisely  constant  mechanical  phase  angle,  4),  between  the  rotor  and 
the  rotating  measurement  coordinate  reference  frame  (position  servo). 
These  two  functions  are,  of  course,  related  by  the  fact  that,  by 
definition: 


A  d<t> 

‘^s  -  dt 

Spin  bearing  disturbance  torques  can  produce  speed  control 
servo  errors  which  result  in  gradiometer  measurement  errors  by 
intioducing  errors  in  the  orientation  of  the  rotating  sensor  coordinate 
reference  frame  relative  to  the  measurement  coordinate  reference 
frame.  Because  of  the  high  "Q"  of  the  difference  mode  resonance  the 
error  due  to  spin  velocity  error  is  by  far  the  dominant  one  of  the  two 
above. 

The  slope  of  the  sensor  signal  phase  at  resonance  is  given  by: 

^  2Q 
dwiu)  "  ‘*’0 

o 


where 


e  =  signal  phase  angle 
Wq  -  difference  mode  frequency 
Q  =  difference  mode  Q 


146 


or,  in  terms  of  spin  frequency,  u  : 

s 

d6  ^ 

du  b) 

8  S 

For  this  analysis  the  speed  control  servo  transfer  function  is 
assumed  to  be  given  by 


<|»(S) 


Si3a 
(S+a)^  ^ 


where 


<|>  =  servo  phase  error 

o  -  servo  corner  frequency 
=  disturbance  torque 
J  =  rotor  polar  inertia 

then,  from  above 

d9  _  d9  ^‘*’1 

dt  dh)  dt 

s 

and 

0  _  2Q  ^ 

w  dt 
s 


d9  d  <>  ^ 

dw  ,.2 
a  dt 


^  d  «{> 

w  j  .2 

s  dt 
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or 


0  (S) 

0  (S) 


s<j>(s) 


2Q  S(S+3q) 

Joj  .3 

s  (S+q) 


If  we  assume  a  si  ep  input  disturbance  torque. 


T 


d 


0  (S) 


Jco 

s 


S+3q_ 

(S+Q)^ 


and  in  the  time  domain 


0  (t) 


ZQ^ 

s 


te  °'^(1+Qt) 


Because  of  the  integration  function  inherent  in  the  gradiometer 
filter  process,  we  are  interested  not  in  the  instantaneous  value  of  9  but 
in  the  mean  value,  0,  over  the  integration  interval,  x.  For  the  worst 
case,  where  the  step  input  occurs  at  t  =  0  and  remains  through  x. 


,  rT  2QT 

=  7  "'‘W'  =  -JZTT 

Jo  s  Jo 


-atj, 
e  dt 


G  =  - 


6QT 


Jw  xa 
s 


/  2  2\ 
1  -  |l  +  Qx  + 


For  representative  values. 


a  =  10  rad/sec,  x  =  10  sec 
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this  redt'ices  to 


e 


Jw  TO 

s 


and  with 


Q  =  300 

5  2 

J  =2x10  gm-cm 

oj  =  110  rad/sec 
s 

0  =  8.  2  X  10"®  T,  rad 
d 

=  4.  7  X  10  ^  Tj  deg 

If  we  wish  to  restrict  9  to  a  maximum  value  of  0.  0013  degrees,  which 

corresponds  to  0.  1  E.  U.  maximum  gradient  error  in  the  earth's  field, 

then  T ,  must  be  limited  to  less  than  280  dyne-cm. 
d 

d.  Spin  Bearing  Specification 

The  spin  bearing  specification  appears  on  the  following 

pages . 
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PROCUREMENT  SPECIFICATION 
FOR 

SPIN  BEARINGS 
FOR 

PROTOTYPE  ROTATING  GRAVITY  GRADIOMETER 

1.  0  SCOPE 

This  specification  covers  the  requirements  for  a  spin-axis 
support  bearing  for  a  prototype  moving-base  rotating  gravity 
gradiometer,  which  may  be  herein  referred  to  as  the  "sensor.  " 
The  sensor  constitutes  the  basic  sensing  element  in  a  system 
designed  to  precisely  measure  gradients  of  the  gravitational 
field  from  a  moving  vehicle.  The  requirements  and  environ¬ 
mental  conditions  associated  with  the  sensor  necessitate  that 
the  spin  bearing  perform  its  function  with  great  precision  as 
well  as  being  rugged,  reliable  and  reproducible, 

2.  0  APPLICABLE  DOCUMENTS 

The  following  documents  in  their  latest  issue  at  contract  date 
form  a  part  of  this  specification  to  the  extent  specified  herein. 

Specifications  — 
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3.  0  REQUIREMENTS 


3.  1  Conflicting  Requirements  —  Any  conflicting  requirements 
arising  between  this  specification  and  any  specifications  or 
drawings  listed  herein  shall  be  referred  in  writing,  to  the 
Hughes  Research  Laboratories  (HRL)  for  interpretation  and 
clarification. 

3.  1.  1  Request  for  Deviation  —  Any  deviation  from  the  require¬ 
ments  specified  herein  shall  be  considered  a  deviation  and  shall 
not  be  allowed  except  by  written  authorization  from  HRL. 

3.2  Materials.  Parts  and  Processes  —  Materials,  parts  and 
processes  used  in  the  design,  fabrication  and  assembly  of  the 
products  covered  by  the  specification  shall  be  in  accordance 
with  sound  and  proven  engineering  and  manufacturing  practices. 
The  manufacturer's  selection  shall  assure  the  highest  uniform 
quality  and  conditions  of  the  product,  suitable  for  the  intended 
use.  The  bearing  parts  and  their  fasteners  shall  have  a 
magnetic  permeability  not  greater  than  1.  010  cgs  units. 

3.  3  Gravity  Gradiometer  Description  —  A  conceptual  design 
sketch  of  the  baseline  configuration  Rotating  Gravity  Gradio¬ 
meter  is  attached  as  Fig,  3.  3.  The  basic  gravity  gradient 
sensor  consists  of  a  crossed  pair  of  mass  quadrupoles  coupled 
by  a  torsional  spring  and  enclosed  in  a  sealed,  evacuated  case. 
This  case  is  then  rotated  at  a  spin  frequency  which  is  adjusted 
to  precisely  one  half  the  inertia- spring  resonant  frequency  of 
the  coupled  mass  quadrupoles.  The  spinning  system  is  enclosed 
within  a  nominally  spherical  shell  which  is,  in  turn,  suspended 
within  its  mounting  frame  with  two  degrees  of  angular  freedom 
for  base  motion  isolation.  The  baseline  sketch  shows  a  hydro¬ 
static  gas  spin  bearing  with  precsurized  gas  supplied  through 
hydrostatic  gas  gimbal  bearings.  The  preferred  configurations 
of  both  spin  and  gimbal  bearings,  however,  are  yet  to  be 
finalized. 
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3.  3.  1  Sensor  Physical  Parameters  —  Preliminary  design 
estimates  of  various  sensor  physical  parameters,  which  influ¬ 
ence  spin  bearing  selection  and  detailed  design,  are  as  follows: 


Rotor  Mass 

7000  gms 

Rotor  Polar  Inertia 

2.  5  X  10^ 

2 

gm-cm 

Rotor  Transverse  Inertia 

2.  0  X  10^ 

2 

gm-cm 

Rotor  Shell  Diameter 

15.  5  cm 

Stator  shell  Diameter 

18  cm 

Stator  Mass 

5000  gms 

Stator  Polar  Inertia 

3.  0  X  10^ 

2 

gm-cm 

Stator  Transverse  Inertia 

3.  5  X  10^ 

2 

gm-cm 

Rotor  Spin  Speed 

1050  rpm 

Motor  Specifications 

Two- Phase  Servo 
Drag- Cup 

Locked  Rotor  Torque  5  x  10“*  dyne-cm  (max) 

4 

Running  Torque  5x10  dyne-cm  (max) 

3.  4  Spin  Bearing  Performance  Requirements  —  Gravity  gradient 
sensor  performance  requirements  impose  specific  performance 
requirements  on  the  sensor  spin  bearing. 

The  sensor  spin  bearing  shall  be  capable  of  meeting  the  perfor¬ 
mance  requirements  set  forth  herein  when  supporting  the  rotor 
mass  and  moment  of  inertia  load  specified  in  paragraph  3.  3.  1 
and  rotating  at  the  spin  speed  specified  in  paragraph  3.  3.  1  while 
the  sensor  is  operating  under  the  environmental  conditions 
defined  in  paragraph  3.  5.  1  after  the  bearing  is  thermally  stabi¬ 
lized.  The  required  thermal  stabilization  time  shall  not  exceed 
the  following  limits. 

Stabilization  Time 

3  Hrs  [1  Hr] 

2  Hrs  [1/2  Hr] 

5  Min  1/2  Min] 
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The  spin  bearing  shall  be  capable  of  meeting  these  performance 
requirements  following  exposure  to  the  conditions  set  forth  in 
paragraph  3.  5.  2  or  3.  5.  3. 

In  addition,  the  spin  bearing  shall  meet  these  performance^ 
requirements  throughout  a  10,  000  hour  operating  life  with  a 
minimum  of  500  rotor  start-stop  cycles.  Furthermore,  the 
bearing  shall  meet  the  performance  requirements  at  any  time 
during  a  minimum  one  year  period  following  assembly  into  the 
sensor. 

In  addition  to  the  values  assigned  to  the  performance  require¬ 
ments,  desire  goals  are  indicated  by  values  in  brackets  [  ]  . 

3.  4.  1  Second  Harmonic  Torque  Ripple  —  Torque  oscillations 

about  the  spin  axis  in  a  narrow  frequency  band  centered  at  twice 

the  spin  frequency  (2w^)  may  cause  significant  errors  in  the 

sensor  output.  To  the  extent  *-hat  these  torque  oscillations  are 

deterministic,  they  can  be  compensated,  however,  the  random 

portion  of  these  torque  oscillations  cannot.  The  deterministic 

portion  is  made  up  of  oscillation  occurring  at  exactly  Zu  and 

s 

whose  phase  is  precisely  fixed  relative  to  the  mechanical  phase 
of  the  spin  bearing.  It  is  required  that  the  magnitude  of  the 
deterministic  torque  oscillation  not  exceed  the  values  specified 
in  paragraph  3.4.  1.  1.  Random  torque  variations  are  character¬ 
ized  by  variations  in  both  amplitude  and  phase  relative  to  the 
above  defined  deterministic  torque  oscillation.  As  a  consequence, 
the  allowable  random  torque  variations  must  be  specified  in 
terms  of  the  magnitude  of  two  mutually  orthogonal  components. 

It  is  requirtd  that  the  standard  deviation  of  the  magnitude  of 
either  of  these  orthogonal  components  within  a  narrow  frequency 
band  centered  at  not  exceed  the  value  specified  in 
paragraph  3.  4.  1.  2. 

3.4.  1 .  1  Deterministic  Torque  Oscillation  —  The  deterministic 
value  of  the  torque  oscillation  shall  be  defined  as  the  average 
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value  over  a  ten  (10)  hour  operation  following  thermal 
stabilization.  This  value  shall  not  exceed  1000  dyne-cm.  In 
addition  the  mean  value  of  the  deterministic  torque  oscillation 
when  averaged  over  the  first  hour  of  operation  following  thermal 
stabilization  shall  not  differ  from  the  ten  (lO)-hour  mean  value 
by  more  than  50  dyne-cm. 

3.4.  1.2  Random  Torque  Variation  —  The  standard  deviation  of 

either  orthogonal  component  of  the  random  torque  variation 

within  a  0.  1  Hz  wide  frequency  band  centered  at  2w  shall  not 

s 

exceed  50  dyne-cm  over  a  ten  (10)  hour  operation  following 
thermal  stabilization. 


4.  2  Torque  Oscillations  at  Other  Frequencies  —  The  root- 
mean- square  value  of  spin-axis  torque  oscillations  within  any 
0.  1  Hz  wide  frequency  band  outside  the  band  specified  in  para¬ 
graph  3.4.  1 .  2  shall  not  exceed  500  dyne-cm  [50  dyne-cm] . 


3,  4.  3  Bearing  Torque 


a. 


Mean  running  torque  shall  not  exceed  5  x  lo'*  dyne-cm 

4 

[2x10  [when  the  mean  applied  load  is  15  lbs. 


b.  Running  torque  sensitivity  to  variation  in  applied 
load  shall  not  exceed  300  dyne-cm/^ound  [lO] 
when  the  mean  applied  load  is  15  lbs. 

c.  Bearing  torque  shall  not  exceed  1  x  10^  dyne-cm 
while  starting  or  during  acceleration  to  the  speci¬ 
fied  spin  speed  with  a  mean  load  of  1  5  lbs. 


3.4.4  Spin  Bearing  Load  Capacity 


a.  Axial  and  radial  load  capacities  shall  not  be  less 
than  50  pounds  [200]. 


Torsional  loa^capacity  shall  not  be  less  than 
50  in-pounds  [lOO]. 
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3.4.5  Spin  Bearing  Compliance 

a.  Axial  and  radial  compliances  shall  not  exceed 
5  X  10  cm/dyne  [lO"^^]. 

b.  Axial  and  radial  compliances  shall  not  differ  by 
more  than  20  percent  flO]. 

c.  Torsional  compliance  shall  not  exceed 

-12 

10  rad/dyne-cm. 

3.4.6  Rotor  Vibration 

a.  The  integrated  power  spectrum  of  bearing- 

induced  rotor  translational  acceleration  shall 

not  exceed  0.  01  cm/sec  rms  fo.  OOl]  in  the 

frequency  range  1/2  (A)  to  4  w  .  Outside  this 

s  s 

frequency  range,  the  power  spectral  density  shall 

'y  A  * 

not  exceed  1. 0  cm  /sec^/Hz  fO.  Ol] 

b.  The  integrated  power  spectrum  of  bearing-induced 
rotor  angular  rate,  normal  to  the  spin  axis,  shall 
net  exceed  5  x  10~^*  sec  [l0"^M, 

3.  4.  7  External  Energy  Requirements  —  If  an  external  energy 
supply  is  required  to  provide  bearing  support,  the  following 
limits  shall  apply: 

a.  For  a  hydrostatic  gas  bearing; 

(1)  Inlet  pressure  shall  not  exceed  100  psig. 

(2)  Total  volume  flow  rate  shall  not  exceed 
1.  0  standard  cubic  feet  per  minute. 

b.  For  an  electrically  actuated  bearing; 

(1)  Voltage  shall  not  exceed  100  volts. 

(2)  Current  shall  not  exceed  1.  0  amp. 

(3)  Power  shall  not  exceed  0.  5  watt  per  pound 
of  load. 
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3.  4.  8  Spin- Axis  Alignment  Reference  -  A  spin-axis  alignment 
reference  shall  be  provided  which  permits  location  of  the  spin 
reference  axis  to  an  accuracy  of  one  mrad. 

3.  4.  9  Electrical  Insulation  —  Though  not  an  absolute  requirement, 
it  is  highly  desirable  that  — [it  shall  be  possible  to  electrically 
insulate  each  of  the  bearing  components  from  its  mechanical 
mounts.  The  insulation  resistance  shall  not  be  less  than  5  meg¬ 
ohms  at  500  volts  and  60  Hz]. 

3.  5  Environmental  Conditions  —  During  the  specified  life 
requirements  of  paragraph  3.4,  and  while  the  spin  bearing  is 
supporting  the  mass  and  moment  of  inertia  rotor  load  specified 
in  paragraph  3.  3.  1,  the  spin  bearing  may  be  subjected  to  the 
following  environmental  conditions.  These  conditions  are  sum¬ 
marized  in  Table  3.  5. 

3.  5.  1  Operating  Performance  Condition—  This  condition 
represents  the  most  extreme  environments  under  which  the 
bearing  is  required  to  operate  and  meet  the  specified  perfor¬ 
mance  of  paragraph  3.4.  When  subjected  to  the  environments 
of  this  condition,  the  bearing  will  be  energized  (if  applicable) 
and  rotating  at  the  specified  spin  speed. 

3.  5,  Z  Operating  Standby  Condition  —  This  condition  represents 
the  most  extreme  environments  under  which  the  bearing  is 
required  to  operate  and  survive  without  damage.  When  subjected 
to  the  environments  of  this  condition,  the  bearing  will  be  ener¬ 
gized  (if  applicable)  and  maybe  either  not  rotating  or  rotating 
at  speeds  up  to  5,  000  rpm. 

3.  5.  3  Non-Operating  Condition  —  This  condition  represents  the 
most  extreme  environments  to  which  the  bearing  may  be  sub¬ 
jected  while  in  a  non-operating  state.  It  must  survive  these 
environments  without  damage.  TJig  non-operating  state  is 
defined  as  not  rotating  and  not  energized  (if  applicable).  A 
rotation  locl<ing  device  may  be  employed  if  necessary. 
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3.  6  Reliability  Objectives  —  The  sensor  spin  bearing  shall  have 
a  reliability  objective  of  0.  995  when  operated  under  the  environ¬ 
mental  conditions  set  forth  in  paragraph  3.  5.  1  at  any  time 
during  10,  000  hours  of  operation. 

4.  0  QUALITY  ASSURANCE  PROVISIONS 

Assurance  that  the  sensor  spin  bearing  meets  the  performance 
requirements  set  forth  in  paragraph  3.4  will  be  provided  by 
means  of  (a)  the  Vendors  Quality  Control  Program;  (b)  an 
adequate  Testing  Program;  and  (c)  a  Reliability  Verification 
Program. 
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Fig.  3.3.  Conceptual  Design  Sketch  For  Rotating  Gravity 
Gradi  ometer. 
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RGG  Baseline  Parameters 


The  follo’ving  parameters  define  the  characteristics  of  the 
Moving  Base  Rotating  Gravity  Gradiometer  as  of  31  July  1972.  All  of 
these  are  subject  to  minor  adjustments  in  the  final  design  and  a  few 
may  have  to  be  changed  significantly  as  the  analysis  progresses. 

a.  System  Parameters 


(1) 

Spin  frequency,  f 

B 

17.  5  Hz 

3 

110  rad/sec 

(2) 

Sensor  difference 
mode  resonant 
frequency,  f^ 

35  Hz 

U) 

o 

220  rad/sec 

(3) 

Sensor  sum  mode 
resonant  frequency 

142  rad/sec 

(4) 

System  integration 

time,  T. 

1 

10  sec 

a.  Sensor  integration 
time,  T 

s 

2. 73  sec 

b.  Filter  integration 
time, 

7 . 27  sec 

(5) 

Se’isor  Q 

300 

(6) 

Inertia  efficiency  of 
each  arm,  q 

0.861 

(7) 

Polar  inertia  of  each 
arm,  C  or  Ij 

35,  610  gm-cm^ 

(8) 

Peak  signal  energy 
stored  in  sensor,  € 

p 

1.23  X  10"^^  joules/EU 

(9) 

Peak  arm  torque, 

1.  53  X  10"^  dcm/EU 

t 
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b. 


Rotor  Parameters 


(1)  Materials 

•  Sensor  arms 

•  Sensor  masses 

•  Sensor  pivots 

•  Rotor  frame 

•  Rotor  end  bells 
Spin  bearing 

•  Balance  screws 


6061  A1 

Mallory  1000 

Be-Cu 

6061  A1 

6061  A1 

not  specified 

6061  A1  and  Be-Cu 


(2)  Rotor  Data: 


Mass 

gms 

XX 

yy 

2 

gm  cm 

zz 

Arm  No.  1 

1,  563 

4,990 

35,  600 

35, 

610 

Arm  No.  2 

1,  563 

35,  600 

4,990 

35, 

610 

Rotor 

(except  arms) 

2,  750 

103,400 

103,400 

132, 

700 

Total  Rotor 

5,  876 

143,  990 

143, 990 

203, 

920 

Stator 

3.  766 

289, 200 

289, 200 

252, 

600 

Rotor  and  Stator 

9,  642 

433,  190 

433,  190 

452, 

500 

c.  Design  Limits 

(1)  Temperature 

•  Operating  control  point 

•  Operating  control  range 

(2)  Arm  mass  unbalance,  Amh 

•  Balance  tube  adjustment 
range  available 


130°F 

±0.005°F 

_4 

±4x10  gm  cm 
±28  X  10"^  gm  cm 
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(3) 

Arm  inertia  unbalance  ratio. 

-6 

kp 

10 

(4) 

Maximum  signal  provision 

dbl  0,000  EU 

(5) 

Signal  level  (not  firm) 

20  to  200  nv/EU 

(6) 

Peak  signal  energy  stored 
in  transducer  -  percent  of 

10  to  20% 

total 

D.  MATERIALS  SELECTION 

1.  Material  Characteristics 

2.  Paramagnetic  and  Diamagnetic  Material  Tests 

These  sections  provide  the  background  upon  which  the  materials 
are  selected  for  construction  of  the  RGG. 
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1. 


Material  Characteristics  (Preliminary) 


a.  General 

Materials  used  in  the  construction  of  rotating  gravity 
gradiometers  (RGG)  should  have  the  following  characteristics; 

Dimensional  stability 

Linear  and  stable  thermal  coefficients 

•Psotropic 

High  anelastic  limit  (high  strength) 

High  density  and  low  density  selection 
High  mechanical  Q  (low  damping) 

Low  magnetic  susceptibility  (nonferromagnetic) 

Workability 

High  and  low  thermal  and  electrical  conductivity  selection. 

A  review  of  Table  IV  Ref.  (4),  (5)  and  (6)  demonstrates  that  no 
materials  have  all  of  the  desired  characteristics  when  examined  min¬ 
utely  and,  furthermore,  data  or  analysis  on  some  of  the  characteristics 
are  not  readily  available.  Many  of  the  parameters  of  many  of  the 
materials  have  not  been  studied  in  depth.  Parameters  of  high  purity 
materials  can  not  be  extrapolated  to  commercial  pure  grades  and 
various  alloys,  fable  IV  Ref.  (5)  in  particular  demonstrates  that  almost 
all  metals  creep  with  age,  stress  and  temperature  cycles.  That  they 
are  all  anisotropic;  the  temperature  coefficients  vary  with  stress,  alloy 
and  operating  temperature;  and  the  elastic  limit  is  largely  dependent 
on  the  accuracy  of  the  measuring  equipment. 

In  view  of  this  situation  the  best  RGG  design  will  result  if  all 
materials  are  viewed  (on  the  microscopic  scale)  as  ferromagnetic 
silly  putty"  until  they  are  proven  otherwise  in  a  commercially  avail¬ 
able  or  high  purity  grade.  The  situation  is  not  actually  as  bad  as  it 
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might  at  first  appear  since  instrument  designers  have  long  been 
concerned  with  basic  material  parameters.  Some  materials  and  pro¬ 
cesses  have  been  tested  extensively  and  most  of  these  are  included  in 
Table  IV,  Ref.  4. 

Table  IVprovides  a  summary  of  many  of  the  characteristics  of 
several  materials.  Since  some  of  the  data  was  obta'ned  from  several 
sources  the  data  are  not  always  completely  consistent.  However,  the 
data  are  believed  to  be  representative  and  usable  for  RGG  design.  The 
tabulations  "Dimensional  Stability  Tests"  and  "Anelastic  Limit  o'^" 
require  comment. 

Dimensional  stability  test  data  was  taken  from  Table  IV  Ref.  (4). 

The  notation;  RT  .  -1  mo-3  mo-12  mo  and  160°F  ^  -1  mo-3  mo-12 

store  store 

mo  refer  to  the  changes  in  the  strain  ina4inchlong  sample  when  stored  at 
room  temperature  and  at  160°F.  The  measured  strains  are  reported  to 
the  nearest  5  pm/m.  Unfortunately,  it  is  not  completely  clear  whether 
the  12  mo  figure  represents  the  change  from  the  initial  value  or  the 
change  from  the  3  mo  value.  However,  a  study  of  the  data  on  many 
materials  indicates  that  the  quoted  value  is  the  change  from  the  initial 
value.  Note  that  a  dash  in  the  data  indicates  an  absence  of  data  and  not 
a  zero  change. 

The  temperature  cycling  tests  associated  with  the  Dimensional 
Stability  heading  represent  10  cycles  between  70°F  and  -95°F  with  a 
30  minute  holding  period  at  -95°F.  The  changes  marked  with  are 
for  10  times  from  +200  to  -100°F. 

The  term  "anelastic  limit"  is  taken  from  Table  IV  Ref.  (13)  and 
the  meaning  is  illustrated  in  Fig.  15.  A  sequence  of  stress  strain 
curves  are  run  to  increasingly  higher  stresses.  The  first  curve  that 
shows  an  open  loop,  but  no  residual  strain  is  called  the  elastic  limit. 

The  first  loop  that  shows  a  residual  plastic  strain  of  10"^  m/m  is  called 
the  anelastic  limit.  The  anelastic  elastic  limit  is  synonymous  with 
"Precision  Elastic  Limit  (PEL),"  "Micr^yield  Strength  (MYS)"  and 
"Micro  Offset  Yield  Stress  (MOYS)"  H  ths-  same  residual  plastic  defor¬ 
mation  base  is  used  in  all  cases.  Co;'i9ideraV*le  care  is  required  in  the 
interpretation  of  data  since  various  authors  are  not  consistent. 
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TABLE  IV 

Properties  of  Materials  (Sheet  1  of  3) 


Properties  of  Materials  (Sheet  2  of  3) 


GOOD  MAXERIAU 

AL*PEL=.  MVS  »MOYS 
^IP  SAME  PLASTIC  STRAIN)  LIMIT 
IS  USED  IW  ALL) 

AMELASTIC 


-^STRAIN)  RESIDUAL 

PLASTIC 

STRAIN) 


POOR  naaterial 


AMELASTIC  LINAIT 


Fig.  15.  Anelastic-Elastic  Limit. 
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b.  Comments  and  Notes  on  Various  Materials  and 

Processes 

The  following  remarks  and  notes  on  various  materials 
and  processes  have  been  obtained  from  a  wide  variety  of  sources  and 
experiences.  They  should  be  considered  in  the  selection  of  instrument 
materiaiE . 

Aluminum;  Aluminum  is  a  generally  acceptable  instrument 
material.  Its  greatest  drawback  is  that  it  is  usually  rated  only  fair 
for  dimensional  stability. 

Beryllium:  Beryllium  is  a  generally  acceptable  instrument 
material.  Its  drawbacks  are  a  health  hazard  during  machining  and 
low  mechanical  Q. 

Stainless  Steel:  Stainless  steel  is  not  acceptable  if  magnetic 
field  sensitivity  is  a  problem.  310  is  the  best  (p  =  1.  0018)  and  the 
permeability  is  not  significantly  increased  by  cold  working. 

Brass;  Brass  has  been  used  in  various  types  of  instruments 
for  years  and  is  generally  considered  to  be  "nonmagnetic".  Most 
brasses  if  made  from  pure  element  materials  do  have  a  permeability 
of  1.  00001  or  less.  However,  nearly  50%  of  all  brass  products  are 
made  from  reclaimed  scrap  and  thus  the  impurities  are  poorly  con¬ 
trolled.  Brass  has  little  to  recommend  it  as  an  instrument  material 
except  its  easy  availability.  Its  creep  strength  is  poor,  its  stability  is 
poor  and  government  specifications  do  not  provide  good  imputity  con¬ 
trol.  Table  IV  Refs.  (16)  and  (6)  provide  good  insight  to  the  magnetic 
properties.  Beryllium  copper  should  generally  be  used. 

Uranixim;  Uranium  is  notorious  for  its  anisotropic  moduli  and 
temperature  coefficients.  Investigate  carefully  if  it  is  considered. 

Tungsten  Alloys;  Txingsten  has  a  great  affinity  for  iron  and  its 
alloys  are  frequently  slightly  ferromagnetic.  Tungsten  alloys  should 
have  the  highest  possible  purity  and  should  be  batch  tested. 

Ni-Span  C;  Ni-Span  C  is  a  material  that  has  been  especially 
developed  to  have  a  constant  modulus  of  elasticity  over  a  relatively 
wide  temperature  range.  It  has  a  significant  magnetic  permeability 
and  is  usable  only  in  special  cases  in  precision  instruments. 
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Impurities  in  Alloys;  Commercial  and  government  alloy 
specifications  are  checked  for  the  desired  alloy  element  limits  only. 
Impurities  are  not  checked  unless  contamination  is  suspected.  Thus 
commercial  or  government  specification  alloys  may  contain  unknown 
impurities.  In  critical  applications  the  most  economical  approach  may 
be  to  purchase  to  the  government  specification  and  then  have  the  batch 
checked  for  the  impurity  level. 

Impurities  on  Alloys:  Machine  tool  wear  leaves  significant 
traces  of  ferromagnetic  impurities  on  the  surface  of  otherwise  para¬ 
magnetic  or  diamagnetic  materials.  In  critical  applications  all  finished 
parts  should  be  given  a  preferential  etch  to  remove  this  contamination. 

Residual  Stresses  in  Stock  Material:  Many  metals  are  hardened 
and  strengthened  by  severe  cold  working  and/or  by  a  solution  heat  treat 
and  quench.  If  a  large  billet  is  cold  reduced  75%  by  rolling  into  a  bar, 
the  residual  stress  near  the  surface  may  be  twice  that  near  the  center. 
If  the  top  half  of  the  bar  stock  is  milled  away  the  finished  piece  wil 
curl  significantly  immediately  and  continue  to  curl  with  time.  The 
same  type  of  creep  can  be  produced  by  unequal  stresses  introduced  by 
quenching  during  a  solution  heat  treatment.  The  most  stable  parts  =»re 
produced  by  symmetrical  machining  -^f  the  stock  and  by  minor  cold 
working  (usually  stretching)  to  equalize  the  thermal  induced  stresses. 

Stability  with  Age  and  Temperature  Cycles:  Many  metals  are 
hardened  and  strengthened  by  precipitation  of  alloy  elements  into  the 
base  material  crystalline  structure.  Precipitation  hardening  may  be 
used  alone  after  a  solution  heat  treat  (beryllium  copper)  or  in  addition 
to  cold  work  and  solution  heat  treat  (2024  aluminum).  The  precipitation 
may  occur  at  room  temperature  and  unless  artifically  aged  it  can 
continue  for  years  at  room  temperature.  Precipitation  for  maximum 
hardness,  ductility  and  strength  usually  occurs  at  a  specific  temper¬ 
ature  and  the  material  can  be  held  at  this  optimum  temperature  for 
extended  periods  of  time  with  little  further  change  in  characteristics. 
However,  if  the  material  is  hardened  at  a  temperature  higher  than 
optimum,  the  strength,  hardness  and  ductility  will  at  first  increase  as 
a  function  of  treatment  time,  peak  out  and  then  decrease  as  a  function 


of  time.  This  is  called  overaging.  Experience  has  shown  that  overaged 
precipitation  hardened  materials  are  more  stable  with  time  and  temper¬ 
ature  cycling  than  those  given  an  optimum  strength  heat  treat. 

c .  Details  on  Selected  Materials 

The  following  paragraphs  provide  detailed  information 
for  the  materials  listed  in  Table  IV. 

(1)  Aluminum; 

Aluminum  can  be  obtained  in  the  following  mill  forms: 

•  Unalloyed  ingot 

•  Casting  alloy  ingot 

•  Alloyed  ingot 

•  W  rought  alloys 

•  Heat-treatable 

•  Non-heat-treatable 

Unalloyed  aluminum  is  very  difficult  to  machine  and  has  low 
strength.  It  is  not  appropriate  for  precision  instrument  fabrication 
except  possibly  as  an  electrical  conductor. 

Casting  alloys  do  not  offer  any  superior  characteristics  for 
instrument  manufacture  and  will  not  be  considered  further. 

Alloy  ingot  is  the  base  for  the  wrought  alloy  products.  Except 
in  unusual  circumstances  it  would  be  more  cost  effective  to  purchase 
the  wrought  alloy  form  (sheet,  rod,  bars,  etc.  )  directly. 

The  numerous  wrought  alloys  of  aluminum  have  been  developed 
to  obtain  specific  characteristics  such  as  high  strength,  corrosion 
resistance,  weldability  and  workability.  Since  most  instruments  do 
not  require  especially  high  strength  or  other  unusual  characteristics 
the  parameters  of  chief  concern  are: 

•  Dimensional  stability 

•  Elastic  moduli  stability 


•  Thermal  expansion  stability 

a  Low  magnetic  susceptibility 

•  High  anelastic  limit 

•  Workability. 

The  non-heat-treatable  alloys  attain  their  strength  and  work¬ 
ability  by  severe  cold  working.  The  resulting  internal  stresses  can 
not  be  relieved  by  heat  treatment  and  because  of  this  they  are  not 
particularly  dimensionally  stable.  These  non -htat-treatable  alloys 
are  considered  to  be  generally  unsuitable  for  i  nstrument  manufacture. 

The  heat-treatable  wrought  alloys  usually  attain  their  best 
strength,  workability  and  other  characteristics  by  a  combination  of 
cold  working  find  heat-treatment.  Alloy  AA-6061  requires  little  cold 
work  and  responds  well  to  heat  treatment.  It  has  beer,  selected  as  one 
of  the  best  for  precision  instrument  work. 

The  chemical  composition  limits  for  AA-6061  from  Table  IV 
Ref.  (1)  and  Federal  Specification  QQ-A-250/11E  (27  Aug  71)  are: 


Silicon 

0.40-0.  80 

Iron 

0.  7 

Copper 

0.  15-0.40 

Manganese 

0.  15 

Magnesium 

0.  80-120 

Chromium 

0.  15-0.  35 

Zinc 

0.  25 

Titanium 

0.  15 

Other  elements: 

Each 

0.  05 

Total 

0.  15 
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Aluminum  Remaiiider 


Nctes:  Composition  is  in  percent  maximum  unless  a 

range  is  shown.  Analysis  is  regularly  made 
only  for  the  elements  for  which  specific  limits 
are  shown. 

The  best  temper  and  stabilizing  heat-treatm':'nt  for  6061  has 
not  yet  been  definitely  established.  However,  the  following  should  be 
near  optimum.  The  material  should  be  purchased  in  the  T-451  temper. 

In  this  temper  it  has  been  solution  heat  treated  at  985±10°F  for  4  hours, 
quenched,  stretched  1  to  3%  and  naturally  aged  (stored  at  room  temper¬ 
ature).  After  finish  machining  it  should  be  overaged  by  precipitation 
hardening  at  400±10°F  for  24  hours.  The  overaging  will  reduce  the 
ultimate  tensil  strength  and  yield  strength  slightly  below  the  maximum 
possible  values  but  it  will  also  improve  the  creep  stability.  The  X51 
in  the  temper  designation  designates  the  stretching  to  help  equalize  the 
cold  work  and  quenching  strains.  Finally  the  part  shall  be  temperature 
cycles  10  times  from  +200°F  to  -95°F  with  a  holding  period  of  not  less 
than  10  minutes  at  each  temperature. 

The  magnetic  susceptibility  of  the  aluminum  alloys  is  nearly  the 
same  as  that  of  pure  aluminum.  As  shown  on  page  176,  Vol  I,  of 
Table  IV  Ref.  (2),  some  of  the  alloy  materials  reduce  the  susceptibility 
slightly  while  others  increase  it  slightly.  Iron,  if  not  in  excess  of  a 
few  percent,  combines  with  the  aluminum  to  produce  F^Af  ^  which  is 
paramagnetic  to  about  the  same  extent  as  pure  aluminum.  Mr.  W.  C. 
Slepy  of  Alcoa  Research  Laboratories  states  in  a  private  communication, 
"Atypical  U.  S,  Government  magnetic  property  specification  reads: 

"The  magnetic  permeability  of  the  material  is  not  to  exceed  1.  001  cgs 
units  in  a  30  orsted  field.  "  All  of  Alcoa's  commercial  alloys  meet 
this  specification  quite  easily;  6061,  for  example,  has  a  magnetic 
permeability  of  1.  00002  cgs  units.  None  of  Alcoa's  commercial  alloys 
are  ferromagnetic.  They  are  weak  paramagnetic  materials  whose 
permeabilities  are  in  the  range  1.00002  -t  0.  000005  cgs  units.  Thus 
the  only  problem  might  be  the  accidental  contamination  of  a  batch  or 
billet  not  tested.  For  extremely  critical  and  high  cost  experiments 
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the  most  reasonable  procedure  would  be  to  purchase  Government 
specification  material  (which  could  be  as  high  as  1.  001)  and  have 
samples  tested  for  iron'^nd  susceptibility. 

The  maximum  operating  temperature  is  based  entirely  on  known 
stability  at  this  temperature  from  their  use  in  precision  gyroscopes 
and  accelerometers. 

No  actual  data  is  available  on  the  anelastic  limit  but  it  is  in  the 
neighborhood  of  1,  000  psi. 

(2)  Beryllium; 

The  characteristics  of  beryllium  vary  widely  with  the 
method  of  manufacture  and  purity.  Table  IV  Refs.  (3)  and  (8)  provide 
many  details.  The  most  common  impurity  is  beryllium  oxide  and  its 
effect  is  not  especially  large.  The  material  selected  here  is  a  medium 
grade  material,  readily  available  from  several  sources.  Higher  strength 
and  higher  purity  are  readily  available.  Nuclear  purity  grades  are  not 
necessary. 

Beryllium  dust  is  highly  toxic  and  the  effect  is  cumulative. 

Proper  facilities  for  machining  -i.re  imperative.  Hand  "dressing", 
fitting  or  scraping  must  not  be  done.  However,  the  finished  parts  are 
perfectly  safe  to  handle  and  assemble. 

Actual  numbers  are  not  available  but  beryllium  is  known  to 
have  high  internal  damping.  It  is  not  useful  as  a  spring  material  in  a 
high  "Q"  system. 

The  recommended  heat  treatment  and  operating  temperatures 
are  based  on  the  known  characteristics  from  their  use  in  precision 
gyroscopes  and  accelerometers. 

The  thermal  coefficient  of  expansion  quoted  is  for  this  particular 
grade  and  process.  It  varies  rather  widely  with  other  purities  and 
processes. 
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(3) 


Beryllium- Copper 


Beryllium-copper  is  a  strong  stable  material  that  is 
slightly  diamagnetic  and  is  easily  worked.  The  percent  composition 
from  Table  IV  Ref.  (6)  and  QQ-C-5306  is: 

B&ryllium  '  1.8  -2.0 

Cobalt  0.  20  min. 

Nickel  plus  cobalt  plus  iron  0.  60  max. 

Copper  plus  beryllium  99.  50  min. 

plus  elements  named 

Note:  Analysis  shall  be  made  regularly  only  for  the 
elements  specifically  mentioned. 

Beryllium-copper  is  a  precipitation  hardened  alloy.  Alloy 
No.  172,  as  reived  in  the  quarter  hard  condition,  has  been  solution 
heat  treated,  sometimes  called  solution  annealed,  at  1450  ±  5°F  for 
several  hours,  quenched  and  then  cold  worked  to  produce  the  desired 
quarter  hard  temper.  After  finish  machining  it  must  be  further  precip¬ 
itation  hardened  to  attain  the  desire  ultimate  properties.  The  standard 
heat  treatment  to  attain  the  ultimate  strength  is  from  MIL-H-7199A, 

2  hours  at  600  ±  5°F.  However,  for  maximum  dimensional  stability, 
overaging  is  desired  and  the  specified  treeiment  is  4  hours  at 
650  ±  5  F.  This  overaging  will  cause  some  reduction  in  ultimate 
tensile  strength  and  some  reduction  in  hardness  but  these  are  not 
serious  drawbacks. 

The  magnetic  susceptibility  of  the  quarter  hard  precipitation 
hardened  alloy  is  quoted  in  Table  IV  Ref.  (6)  as  -0.  62  x  10"^  cgs  (emu) 
units.  The  magnetic  effects  of  impurities  on  copper  base  alloys  is 
discussed  in  considerable  detail  in  Table  IV  Ref.  (16).  This  reference 
indicates  that  as  long  as  the  iron  content  is  less  than  0.  15  percent 
the  susceptibility  will  be  no  greater  than  1.5x10"^  cgs  units.  Thus, 
for  critical  applications  each  batch  should  be  purchased  to  the  govern¬ 
ment  specifications  and  then  tested  for  magnetic  properties. 


181 


(4)  Elgiloy; 

Elgiloy  was  developed  to  provide  a  very  strong,  stable, 
nonmagnetic  material  for  precision  spring  applications.  It  is  manu¬ 
factured  and  distributed  by  the  Elgiloy  Company,  853  Dundee  Avenue, 

Egin,  Ill.  60120.  The  composition  and  characteristics  shown  in  Table  IV 
were  obtained  from  Specification  No.  54-71A  of  the  Elgiloy  Company. 

Elgiloy  has  one  characteristic  that  limits  its  use  in  RGG 
applications.  The  stock  sizes  of  the  rod  (wire)  are  limited  to  about 
0.130  inches  diameter,  maximum.  This  limit  is  due  to  the  severe 
cold  working  required  to  attain  the  desired  properties.  It  can  be 
obtained  in  strip  up  to  four  inches  wide  and  a  thickness  of  about 
0.130  inches  maximum. 

(  5)  310  Stainless  Steel; 

This  stainless  steel  can  be  considered  where  minimum 
magnetic  permeability  is  not  necessary.  Its  permeability  is  approxi¬ 
mately  p  =  1 .  002  and  in  different  specification  ranges  from  1.0012  to 
1 . 003.  The  permeability  varies  only  slightly  with  cold  working  up  to  50%. 

The  biggest  uncertainty  with  310  is  that  the  permeability  is 
not  guaranteed  by  any  specification  located  so  far.  One  possible  alter¬ 
nate  is  Carpenter  No.  10  which  has  a  guaranteed  maximum  p  =  1.005 
from  the  annealed  condition  up  to  50%  cold  reduction. 

(6)  Mallory  1000; 

This  material  has  been  developed  by  P.  R.  Mallory  and 
Company,  Inc.  Its  general  characteristics  are  high  density,  high 
strength,  nonmagnetic  and  corrosion  resistant. 

The  data  in  Table  IV  were  obtained  from  that  table's  Ref.  (14). 

It  is  seen  that  except  for  the  low  thermal  coefficient  of  expansion,  which 
does  not  match  that  of  aluminum,  it  is  generally  an  excellent  material 
for  end  masses . 

A  word  on  the  magnetic  characteristics  is  in  order.  P.R.  Mallory 
states  that  the  material  is  nonmagnetic  but  they  could  not  even  estimate 
its  permeability  or  susceptibility.  Table  IV  Ref.  (15)  discusses  this 
Cu(4)  -  Ni(6)  -  W(90)  alloy  in  some  detail,  but  fails  to  mention  its 
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susceptibility.  However,  P.  R.  Mallory  states  that  the  copper -nickel - 
tungsten  forms  a  ternary  that  is  Cu(4)  -  Ni(6)  -  W(18)  in  terms  of  the 
original  composition.  The  remaining  72  percent  tungsten  is  in  the  form 
of  rounded  grains  of  pure  tungsten  embedded  in  the  ternary.  There  is 
no  porosity  and  the  alloy  can  be  easily  and  accurately  controlled. 

Table  IV  Ref.  (16)  shows  that  a  copper  nickel  alloy  of  Cu(40)  -  Ni(6)  has 
a  susceptibility  of  approximately  2.8xl0"^cgs  and  iron  contamination 
up  to  0.  15  percent  has  little  effect.  There  is  little  reason  to  suspect 
that  the  addition  of  tungsten  to  this  binary  alloy  would  chan^ie  its 
susceptibility  or  its  ability  to  deactivate  small  amounts  of  iron.  Thus, 
we  could  expect  the  permeability  of  Mallory  1000  to  be  nearly  equal 
to  that  of  pure  tungsten. 

Qualitative  measurements  were  made  on  a  sample  of  Mallory  1000 
and  it  was  found  to  have  a  permeability  slightly  greater  than  that  of 
6061  aluminum.  Therefore,  the  permeability  of  Mallory  JOOO  was 
estimated  to  be  1. 00007,  slightly  greater  than  that  of  pure  tungsten. 
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2. 


Paramagnetic  and  Diamagnetic  Material  Tests 


a.  Introduction 

Obtaining  susceptibility  or  permeability  data  on 
paramagnetic  and  diamagnetic  structural  material  is  difficult  and  time 
consuming.  The  relatively  few  researchers  who  have  the  equipment  to 
make  such  tests  are  interested  in  the  theoretical  aspects  of  magnetism 
and  not  in  the  permeability  of  commercial  materials.  Furthermore 
even  if  data  is  obtained  on  a  commercially  available  material  no  one 
will  guarantee  the  permeability  of  another  similar  commercial  batch. 
However  a  number  of  commercial  materials  have  very  low  permeabili¬ 
ties  and  are  suitable  for  use  in  sensitive  instruments  if  they  are  not 
acc’  lently  co.itaminated.  A  simple  qualitative  test  for  accidental  fer¬ 
romagnetic  contamination  of  low  permeability  materials  was  desired. 

b.  Test  Apparatus 

A  large  horeshoe  permanent  was  fitted  with  pole  pieces  to 
provide  a  gap  approximately  0.  6  inches  long  and  a  cross  section  near 
the  center  of  about  0.  125  square  inches.  The  magnetic  field  in  this 
gap  was  measured  and  found  to  be  2,  000  oersteds  near  the  center  and 
3,  000  oersteds^near  the  tips  of  the  pole  pieces.  The  average  linear 
field  was  estimated  to  be  2,  500  oersteds  and  the  effective  gradient 
2,  000  oersteds  per  centimeter. 

Material  samples  were  made  0.  5  inches  long  and  0.  125  inches 
in  diameter.  These  samples  were  suspended  in  the  magnetic  field 
between  the  pole  pieces  by  a  very  thin  fiber,  about  10  inches  long, 
attached  at  the  center  of  the  sample  so  that  the  samples  acted  as  a 
"compass"  needle  in  the  established  field. 


C.  Material  Test  Results 

A  number  of  material  samples  were  prepared  and  tested 
in  the  magnet  gap.  The  speed  and  direction  of  the  sample  reaction  to 
the  magnetic  field  was  used  to  estimate  the  sample  susceptibility. 

The  following  equations  are  from  Subsection  B-4,  "RGG  Torques 
Due  to  Paramagnetic  Materials.  "  The  torque  equation  is  corrected 
to  account  for  the  homogeneous  sample. 


1^ 

Torque 


1  +  4TrXp 


X  p  Vol  sin  6  cos  6 


where: 

p,  =  permeability 

X  =  sa?ceptibility 

p  =  density 

_2 

Vol  =  volume  of  the  sample,  2.  5  x  10  cc 

e  =  the  angle  the  sample  makes  with  the  field 
H  =  magnetic  field  strength,  2500  oersteds 
dH/dX  =  gradient  of  the  field,  2000  oersteds/cm 

I  =  length  of  the  sample,  1.27  cm 

The  peak  torque  for  the  field,  gradient  and  sample  size 
previously  given  is: 

Torque  =  xp  x  0.  322  x  10^  peak  dyne  cm 

The  design  susceptibilities  and  qualitative  observations  for 
several  materials  follow. 
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Substituting  for  I  and  t,  k  is  found  to  be 

k  =  0.  02  d  cm/ rad 

It  should  be  possible  to  rotate  the  upper  end  of  the  suspension 

G  =  r  =  =  82  rad  =  13  revs, 

k  0.  02 

But  the  torques  and  the  revolutions  that  the  suspension  can  be  rotated 
are  not  nearly  this  large. 

d.  Conclusions 

The  apparatus  can  distinguish  paramagnetic,  diamagnetic 
and  slightly  paramagnetic.  This  is  all  that  is  required  to  identify 
accidental  contamination  of  material  batches. 
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Material 

X  X  10° 

P 

Observation 

Crown  Glass 

-0.  90 

2.47 

Slightly  diamagnetic 

6061  Alum. 

+0.  63 

2.  70 

Slightly  paramagnetic 

Be-Cu 

-0.  63 

8.  2 

^Very  slightly  diamagnetic 

Elgiloy 

+0.  39 

8.  3 

Paramagnetic 

310  Stainless 

+20 

7.  9 

Strongly  paramagnetic 

Malloy  1000 

+0.  33 

17.  0 

Paramagnetic 

^Se''  followint-  paragraphs 


The  theoretical  peak  torque  ‘‘or  beryllium  copper  3e-Cu)  is 
Torque  =  -0.63  x  8. 2x0. 322  =  -1.65  dyne  cm 

When  the  Be-Cu  sample  was  first  tested  it  would  have  been 
classified  as  "slightly  paramagnetic.  "  After  etching  to  remove  tool 
marks  it  was  very  slightly  diamagnetic.  It  is  not  as  strongly  diamag 
netic  as  the  crown  glass  so  it  is  probably  slighLly  contaminated. 

Some  idea  of  the  sensitivity  of  this  test  can  be  obtained  by 
calculating  the  stiffness  of  the  suspension  fiber.  The  period  of  a 
torsional  pendulum  is 

t  =  2TTyi7k 

and  the  measured  period  was  15  sec.  The  momei/  of  inertia  of  the 
cylindrical  sample  suspended  at  its  center  of  gravity  is 


=  0.  1 1  5  gram  cm^ 


Material  Test  Results 


c . 

A  number  of  material  samples  were  prepared  and  tested 
in  tlie  magnet  gap.  The  speed  and  direction  of  the  sample  reaction  to 
the  magnetic  field  was  used  to  estimate  the  sample  susceptibility. 

The  following  equations  are  from  Subsection  B-4,  "RGG  Torques 
Due  to  Paramagnetic  Materials.  "  The  torque  equation  is  corrected 
to  account  for  the  homogeneous  sample. 

HL  =  1  H  4irXp 

Torque  =  X  p  Vol  sin  0  cos  0  j  ^  ^ 


where : 


^l  =  permeability 
X  =  susceptibility 
p  ~  density 

Vol  =  volume  of  the  sample,  2.  5  x  10"^  cc 

e  =  the  angle  the  sample  makes  with  the  field 

H  =  magnetic  field  strength,  2500  oersteds 

dH/dX  =  gradient  of  the  field,  2000  oersteds/cm 

t  =  length  of  the  sample,  1.27  cm 

The  peak  torque  for  the  field,  gradient  and  sample  size 
previously  given  is: 

Torque  =  Xp  x  0.  322  x  10^  peak  dyne  cm 

The  design  susceptibiliues  and  qualitative  observations  for 
several  materials  follow. 
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